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NOTICE TO AUTHORS 


1. Communications._-Papers must be communicated to the Society by a Fellow. ‘They 
prcenenly sicse meget atc ampes. Mr agers scat. Sree conveying briefly the 
a , and drawing attention to importan one ieicaeiiae and to the main 
conclusions, should be intelligible in tm in itself, wi without reference to the paper, 
to a reader nH aon Seaedadee of the subject; it should not normally exceed 200 words 
in | Authors are requested to submit MSS. in duplicate. These should be 
mas ens eee margin of not ese than one incl onthe 

Corrections to the be made in the text and not in the 

Unless a paper reaches the Secretaries more than seven days before a Council 

cielia it will not normally be considered at that meeting. By Council decision, MSS. of 

accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 1g, 239, 1938. 

H. Jeffreys, Theory of Probability, and edn., section 5.45, p. 258, Oxford, 1948. 

Alternatively, the Harvard reference system may be used. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 


recommendations of Commission 3 of the International Astronomical Union (7vans. 
1.A.U., Vol. VI, P. 345, 1938). Council has decided to t the 1.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Sym Signs and Abbreviations 
— : Royal Society, 1951) except where these conflict with LA.U. practice. 


4 rams.—These should be designed et ete Sa 
about ees cee ceded bs cae ad eae for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alteration exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When pa pers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 
8. When submitti po authors are requested to indicate whether they will be 


willing and able to paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysica! Discussions or both should notify the praritheel Secretary. 
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MEETING OF 1956 MARCH 9 
Professor Sir Harold Jeffreys, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Vishnu Gianesh Bhide, Ph.D., Royal Holloway College, Englefield Green, 
Surrey (proposed by W. H. McCrea); 

Rodney Deane Davies, 18 Granville Road, Wilmslow, Cheshire (proposed by 
A. C. B. Lovell); 

Robert Allen Rach, Lick Observatory, Mount Hamilton, California, U.S.A. 
(proposed by W. P. Bidelman); 

Ragihubir Saran Sharma, Ph.D., Royal Holloway College, Englefield Green, 
Surrey (proposed by W. H. McCrea); and 

Stanislavs Vasilevskis, Lick Observatory, Mount Hamilton, 


California, 
U.S.A. (proposed by W. P. Bidelman). 


The election by the Council of the following Junior Members was duly 
confirmed : 


Brian Geoffrey Marsden, Fixby, 7 Perne Road, Cambridge (proposed by 
D. W. Dewhirst); 

Christopher John Millns, The Brow, ‘The Brecks, Rotherham, Yorkshire 
(proposed by Z. Kopal); and 

Brian Warner, 2 Diamond Cottages, Crawley Down, Sussex (proposed by 
P. Moore). 


Ninety-four presents were announced as having been received since the 
last Meeting, including : 


Sir Harold and Lady Jeffreys, Methods of Mathematical Physics, 4rd edition 
(presented by the authors); 

G. O. Jones, J. Rotblat and G. J. Whitrow, Atoms and the Universe (pre- 
sented by G. J. Whitrow); 

Sir John Herschel, Outlines of Astronomy, in Chinese (presented by Sir Eric 
Millar); and 

Le Livre du Centenaire de la Naissance de Henri Poincaré (presented by 
Le Comité d’Organisation du Centenaire de la Naissance de Henri 
Poincaré). 
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MEETING OF 1956 APRIL 13 
Professor Sir Harold Jeffreys, President, in the Chair 

The election by the Council of the following Fellows was duly confirmed : 

Adly Salama Asaad, B.Sc., Helwan Observatory, Egypt (proposed by 
C. W. Allen); 

Wilbur Norman Christiansen, D.Sc., C.S.1.R.O. Division of Radiophysics, 
University Grounds, Sydney, N.S.W., Australia (proposed by R. G. 
Giovanelli) ; 

Bernard Yarnton Mills, B.Sc., Radiophysics Laboratory, University Grounds, 
City Road, Chippendale, N.S.W., Australia (proposed by R. G. 
Giovanelli) ; 

Ali M. Naqvi, A.M., M.Sc., Ph.D., Professor of Physics and Astronomy, 
University of Sind, Hyderabad, W. Pakistan (proposed by M. K. Vainu 
Bappu); 

Kevin H. Prendergast, Yerkes Observatory, Williams Bay, Wisconsin, 
U.S.A. (proposed by 5. Chandrasekhar); 

Elizabeth Roemer, Ph.D., Leuschner Observatory, University of California, 
Berkeley 4, California, U.S.A. (proposed by W. P. Bidelman); 


Amos Hudson Ronan, M.C., Highover, Westerham, Kent (proposed by 
C. A. Ronan); 


Alexander George Stanford, ‘Thraves, Beech Close, 
(proposed by A. Hunter); and 


William Walton Wilson, 28 Norley Drive, Sale, Cheshire (proposed by 
Z. Kopal). 


Cobham, Surrey 


The President announced ‘that the Council had elected the following 
Associates of the Society : 

Julius Bartels, Professor of Geophysics at the University of Gottingen, 
Germany ; 

Mrs. Cecilia Payne Gaposchkin, Phillips Astronomer at Harvard Observatory, 
Cambridge, Massachusetts, U.S.A.; and 

Charles Donald Shane, Director of the Lick Observatory, Mount Hamilton, 
California, U.S.A. 


Ninety presents were announced as having been received since the last 
Meeting, including : 

C. 'l. Moss, The development of Spectroscopic Astronomy (presented by the 
author); 

R. H. Phillimore, Historical records of the Survey of India (presented by the 
Surveys-General of India); 

Liége University, Institute of Astrophysics, Atlas of representative cometary 
spectra (presented by the International Astronomical Union); 

E. F. Savarensky and D. P. Kirnes, Elements of Seismology and Seismometry ; 
Seismological Council, Academy of Sciencies, U.S.S.R., Guide to the 
performing and working of observations at seismic stations of the U.S.S.R.; 
Geological Institute of the Academy of Sciences, U.S.S.R., Questions of 
the Earth’s internal structure and development (these three, in Russian, 
presented by Colonel E, Tillotson); and 

Three volumes of astronomical pamphlets, presented by Sir Eric Millar. 





SIX VISUAL DOUBLE STAR ORBITS 
W. D. Heintz 
(Communicated by R. v. d. R. Woolley) 
(Received 1956 February 29)* 


Summary 


Orbital elements and physical data of the pairs ADS 5949, 8949, 9392, 
9425, 14412, 16057 are derived. 





A 2855 (ADS 5949. 7% 14™-2, +1° 4’) 
Elements: P 71%°+3 eo'16 i+ 135° 
(precession neglected)  7'1953°5 a o”*40 8 42°: 
Observations P, P, O-C 


1914°45 A 313°0 313°4 
20°93 A 2811 281°0 
27°97 VBs 253°2 252°6 
36°44 B 225°9 224'2 
38°54 VBs 214°9 216°7 
42°06 VBs 200°9 202°6 
46-00 VBs 182°4 181°5 
50°58 VBs 143°6 144°4 
54°23 Bz see wee eee round 


Ephemeris: 1956°0 93° 
55°0 79°: 
60°0 67°: 
62°0 7s 
64°0 45° 
66-0 40°5 


° 
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° 
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oyu nN OS 


This very close couple is within reach only of the largest apertures, and the 
observations are not numerous. ‘lhe predicted places diverge rapidly from those 
following from Baize’s elements (Circ. d’ Inf. 4). 

21757 (ADS 8949. 13" 31™-7 


Elements : P 334 € 0°24 
(1900) T 1723°0 a 2°35 


Normal places : P,(1900) } O-C 
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* Received in original form 1955 July 6. No changes in the resulting orbits have been made 
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Ephemeris 1956°0 101° 
(equinox of date): 60°0 103°: 
64°0 105" 
68-0 107° 
72°0 109°: 
76°0 111° 
$0°'0 114°3 


The point of maximum distance has been passed, and the companion closes in 
slowly. ‘The observations show marked deviations from a constant areal velocity, 
but so far there is no periodicity indicated by the residuals. Also, the run in the 
distance differences could be partly due to error of the interpolation curve as 
the older measurements are unusually discordant. Accordingly, the recent 
elements published by Florsch (Cire. d’Inf. 3) are based on a different areal 
constant. ‘lhe couple has a proper motion and parallax not much different from 
the 19’ distant star ¢ Virginis. ‘The collection of observations in a previous 
paper on © 1757 by Luplau Janssen (2) was used and completed by more than 200 
recent observations, especially the numerous measurements by Prof W. Rabe, 
which were kindly communicated to the writer, up to 1955:0. 


< 1883 (ADS 9392. 14! 46™-4, 46° 10’) 
Elements : P 10279 e062 t+115°°8 Ww 291 °5 
(1900) T 1941°9 a 1"*45 S4 88 6 ft 0°3505 


Normal places : P,(1g00) P, O-C Po 4 O-C 


1830°74 269° 269 °4° °° "£3 ' 0°04 
43°01 2606 267° ° 


‘10 ° 
56° 203 264° 03 ‘O° 93 
65°0° 260°: 261°75 ‘ O'S5 . (ots) 
74° 260°: 255°: ' ‘go 3: 07 
53°: 256: 255°1% ' ‘75 ‘75 ° 
95° 247° 245° ‘ . : 02 

1903 °62 242°: 241° ' s “ 05 
10° 230° 333° o2 
15°55 218° 218°5 “4 , Ol 

197° 196° ; y 2 02 
169° 872°4 ‘ >*2 2°25 ° 
153°5 153°3: "" “4 ; Ol 
139 14! 

125 

116 

114°7 


Ephem«e ris 
(equinox of date): 


660606000 
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Like & 2924 (below), this pair has been difficult during the last decades, but now 
the separation increases. Due to the fairly high inclinations, the orbital elements 
are less certain than would be expected from observed arcs of nearly 180°. The 
computation of these elements was dropped in 1953, after publication of orbits by 
Pensado (% 1883) and Arend (2 2924), cf. the ephemeris catalogue (1), but the 
work was resumed as recent observations seem to show residuals. 


OX 288 (ADS 9425. 14% 51-0, 
Elements : P 215¥-4 e 0°60 
(1900) T 1819°3 a "09 
Normal places : P,( 1900) P. O-C 


1844°385 224° 
43°63 216 
64°40 203° 
73°74 197° 
82°18 195 
85°49 193° 
96°34 19g! 


224°3! 
217°02 
202°41 
1938°0O7 
195°12 
193° 
1Q1° 
189° 
18701 


Seumbrerodnu 


1903 °42 als) 
12°51 137° 
26°23 155° 184°85 
30°54 182° 182°78 
43 1814 151°44 


50° 179°9 179°93 
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A751 (ADS 14412. 20!! 52-5, 
Elements : P 58¥-8 eo'si r w 109 
(precession neglected) T 1918-0 ao”'19g iv pm O12 


Observations Pe r. Po 


1904°53 ¢ 33° 32°7 ; o'16 
08°65 « 24° 5°: 17 
eyes 1 " é; . 20 
14°71 J 339° a 13 
15-16 too clos« 
20°63 « ‘ 15 
21°54 « 5 . 13 

2°79 : 14 

1 te 17 

3°72 2 17 
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“” 


Ephemeris: 1956°0 52°9 0°22 
58°0 47°6 22 
60'0 42°3 22 
62°0 36°7 22 
64°0 30°8 21 
66:0 24°3 o'19 


Also for this pair, Baize (Circ. d’ Inf. 4) gives a recent set of elements, very similar 
to those derived here. Some measurements do not satisfy the law of areas, as 
seen from the leap in the residuals about 1920. ‘The interferometer observations 
1923°74 (3n) by Maggini are apparently spurious and have been rejected. 


X 2924 (ADS 16057. 22! 318-6, 469° 39’) 
Elements : * P 225%6 € 0°27 i +790 w 6°°7 
(1900) T 1975°8 a o”'78 S& 82°-6  1°596 
Normal places : P,(1g00) P, O-C Po P, O-C 
1830°0 257'°5 257°29 
40 260° 259°54 
50 263°1 261°70 
60 264°2 263 °33 
70 265°1 265°20 
$o 266°3 2607°24 
go 268-6 269°69 
1goo" 271°6 272°93 
10 275°9 27794 
20 2886 287°82 
30 318°5 317°79 
40 318 33°33 
50 64°5 63°83 


0°85 0°05 
0°93 07 
0°98 06 
0°99 04 
0°96 ol 
0°90 t+ 03 
o'sI ol 
0°68 o2 
0°52 ° 
0°34 t+ O2 
ors t+ OF 
o'18 t 02 
0°34 + 0°02 


Orne One ne OTOH" 0 
©0000 00000°0°? 


Ephemeris 1956°0 
(equinox of date): 58'0 
60°0 

62°0 

64°0 

66-0 

63:0 

70'0 


Last elements by Arend (1953, 2 orbits). 

The following table gives the physical data of the six systems, total apparent 
magnitude m, magnitude difference Am, and spectrum, from which the dynamical 
parallaxes 7,,, and total masses #l, according to Russell and Moore (3), have 
been calculated. ‘The table also shows the radiation-energy parallax (strahlungs- 
energetische Parallaxe) 7,,, and component masses, computed from tables 
given by Brill (4) and Hopmann (5). (‘lhis method uses as arguments the data 


also employed for the dynamical parallaxes, semi-axis major, period, apparent 
magnitudes and spectra, but also takes into account an empirically corrected 
mass-luminosity relation, thus obtaining individual instead of merely statistical 
parallaxes.) From the latter values the other data have been derived : 


magnitudes M, radii R and mean densities D of the components (unit = Sun), and 
the semiaxis major in A.U, 


absolute 
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Commonwealth Observatory, 
Canberra, A.C.T.: 


1955 June 29. 
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MICROMETER OBSERVATIONS OF DOUBLE STARS 


W. D. Heintz 


(Communicated by R. v. d. R. Woolley) 


(Received 1956 March 10) 
Summary 
Measurements made by a filar and by a Muller double image micromete: 


are given and the errors of both series discussed. 


'. Filar Micrometer Observations.—By courtesy of the Commonwealth 
Observatory, Mount Stromlo, the Oddie refractor (aperture 9”, focal length 139"), 
equipped with a filar micrometer, was made available for double star work in 
the winter of 1955. Krom June 18 to September 10, g25 observations were 
obtained, 873 of which are given below; the remainder consists mainly of negative 
or highly doubtful results of difficult pairs obtained under insufficient atmospheric 
conditions. Also, a few of the very first observations had to be rejected. 
11 observations, also listed below, had been made in August 1954, in order to 
test the micrometer for other purposes. ‘The separation limit of the telescope is 
0”+52. 

The programme consisted of stars, RA 12-24", with known orbits or 
appreciable orbital motion, ‘couples prioritaires’’ from the 1.A.U. Double 
Star Circulars, and, for comparison, stars which are also regularly observed by 
Professor Rabe, Munich. A number of fairly close stars from the programme 
of Woolley, Gottlieb, and Simonow (1), 1945/47, at the same instrument (but 
by another micrometer) have been added. 190 double, 8 triple, 1 quadruple 
system, in total 209 objects, have been measured. 

The angular revolution value resulting from 65 determinations is 
29"°718+0"-027. All measures are made with field illumination of properly 
adjusted brightness, a 1/4” eyepiece, and with eyes parallel to the line joining 
the stars. ‘The mean error of one observation in angle increases from +0°-6 
for the wide couples to 5° for the closest (not separated) ones; in distance it is 
t 0-06 for the close pairs, and 0-11-0-13, practically constant, for all separations 
-1”. ‘The increase in accuracy for bright stars and for couples with equal 
components is small. 

The systematic differences of these observations minus the most recent 
values of 46 couples obtained by W. Rabe, Munich (11” aperture), distributed 
according to brightness of primary and separation, are, after the small epoch 
differences have been made allowance for in all cases where an orbit exists, by 
interpolation from the orbit ephemeris (2) : 

my» p-—-0o"-75 0”°75-2"-00 2”-00 
50 31 | o-oq (2) fOr2 0°02 (2) +o'4 orn (4) 
5'0-7'5 t4°o -+o'12 (7) +10 +004 (17) o'8 oor (3) 
7°5 o'2 o'13 (3) +o°6 0'06 (4) +or7 006 (4) 
The position angles given here are apparently somewhat too large, and this 
difference increases for close couples. In separation, only the close pairs show 
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appreciable differences, due to the fact that Rabe corrects these measures by 
comparison with large-aperture observations. ‘The agreement with his 
uncorrected values is satisfactory. Rabe’s personal equations are well determined ; 
he measures position angles, perhaps, slightly smaller than the mean of all observers 
does (3); his distance system shows no peculiarities. ‘lhe writer's thanks are 
due to Professor Rabe for making his unpublished results available for this 
comparison, which is of interest with respect to the double-image results given 
later. 

In the table of observations, only the means are given. A + sign with the 
number of nights indicates that the star has also been unsuccessfully observed 
on one or more nights. An asterisk refers to the remarks, For the comparison 
with orbits, the ephemeris catalogue (2) has been used, in some cases directly 
the elements in the orbit catalogue (4), where references to the orbits may be 
found. 

For negative observations, the limiting distances (e.g. <0":4) are based on the 
correctness of the ADS or SDS magnitude differences. It may be added that 
a number of these need revision, not only in the tar south. For example, the 
companion of p Cap, given as 1o"™-o, is certainly brighter than 8™ (Wirtz, 1907, 
also estimated 7™-7); similarly the companions of § 227 and A244 seem brighter 
than 10°5 and 9-4, respectively. On the other hand, € Pav b, given as 8"-6, ts 
certainly fainter than 1o™, about the same as 0 PsA B (10°5). 


ADS 


A 
p 


p 


8555 Bb 
8573 B 
h 
8630 & 
Ru 
8708 Ob 
8728 I 


3759 
3864 
5949 


Howe 


goig & 
p 


I lowe 


ys 


~ 


143 
605 
606 


225 
258 
4539 
1670 


4590 
395 
43 
1751 
343 
25 
1735 


y Cen 
y Vi 
p Mus 


x 


1goo 


55°49 
5° 


5° 


Pp 


1194 
155 + 
too difficult 
115°4 
too close 
3304 
62 
z10°9 
15‘! 
205°3 
round 


207°7 
259°4 
182°5 
239'3 
100°3 
252°6 
127°! 
136°6 
333°4 
O31 
round 
Sg°2 
5°6 


wm 
* 


~ 


~-Fe Hwwvyv 
> 


w 


of ee 
* = 


auwwe ws 


> wa w 





9343 
9344 


9392 


9425 
9453 


9579 
9586 


9659 


g7o!l 


9775 
9323 
9836 
g88o 


9909 


9913 
9951 


939 
181g 
4007 

225 


40387 


414 
1865 

345 

309 

1883 

236 

h = 4707 
OX 288 
p 239 


Cape 15 
h 4728 
p 618 
329 
228 
35° 


620 

30 
977 
393 
258 


1995 
204 


947 
120 


AB 
BC 


x Cen 


¢ Boo 


59 Ilya 


m7 Lup 
t Lib 


pe Lup AB 
AB/2 ( 


« Lup AB 
AC 
y Cir 


y Lup 
0 Ser 


«' Nor AB 
AB/2-C 


€ Sco AB 

AB/2-C 

AB 

Bb Sco AC 

v Sco AB 
AC 
CD 


p Oph 


a Sco 


W. D. Heintz 


a 


h 
1411 
13 
18 
22 
22 
32 
30 
39 
39 
39 
43 
40 
45 
5° 
51 
55 


‘ol 
o2 


10 
10 
12 


- goo 


55°52 
49 
51 
50 
50 
54 
51 
54 
49 
54 
45 
52 
55 
52 
52 
55 


5°53 
49 
54 
54 
52 
57 
49 
49 
55 
55 
53 


> 
a 


53 
49 
55 
56 
53 
58 
56 
56 
54 
53 
56 
58 
57 
50 


P 


6 


85°7: 


291°3 
140°2 
295"! 
93°! 
104°2 
379 
347°! 
305'9 
292°2 
round 
110°7 
116°9 
33°7 
177°3 
337°2 


44°3 
74°4 


330°6 
26°3 
141°9 
322°6 
125°7 
169° 
292° 
87° 
251° 
170°9 
459° 
3192 
282°6 
179°5 
3425 
159°7 
40°0 
349°0 
273°2 
276°3 
163°7 
234°9 
245°2 


320°6 
53°7 
153°1 
201°7 
3°2 
330°3 
50°6 
170°7 
390°4 
345°0 
194°4 
277°6 
195°0 


Ww 


e689 6588 55% 


tS 
t 


bt 


7 


ox 0 


~- 


Ow NK NS HW 


= WN 


yo 2 ees 


I 
°° 
o: 
oO: 
6° 


Con onoonee 


o 
~I 


Vol. 116 


* 


PW AW BWtFOUW 
* . * j 


= 
L 
- 


* 


NUS WwW 
+ 


* 


n+ Qwr FWOWHWHETHOUM NSF SHS MNWFEWHWFEYH 


> 
* 


> N 
* 


* 


* 


aust et Ute wFN 





10229 
10230 
10257 
102605 
10312 


1034! 


10355 
10374 
10417 
10465 
10564 
10595 
10657 
10696 


10828 
10858 


10912 


11005 
11046 


eeee! 
11186 


11324 
11353 
11354 
11378 


11483 
11484 


11530 
11622 


11643 
11950 
11971 


11989 


Micrometer observations of double stars 


21 Oph 


24 Oph 


7 Oph 
36 Oph 


AB 


Tt Oph 
70 Oph 


73 Oph 
& Pav 


59 Ser 


Stone 
Or 


Ox 
Ho 
A 
B 
Ru 
h 
Hd 


s 
~ 


a 


h 
16°28 
24 
36 
39 
45 
49 
52 


54 


1900 + 


55°57 
55 
57 
61 


P 


329°7 


166°5 

37°4 
196'5 
130°4 
352°3 
293°9 


177°3 


82°5 
175°5 
round 
343°! 
166°0 
254°5 
230°3 
263°8 
113°7 
290°7 
2388 + 
unmeasurable 
164°3 
195°7 
202°5 
287? 
154°0 
258'1 


269°7 
97° 
29°" 
29° 
94° ; 
154°5 
358°0 
3181 
253'2 
103°5 
135°4 
174°5 
173°0 


339°2 
44°4 
round 
270'0 
114 
242'5 


1061 
15°° 

221°9 
37°0 

242 
30°4 


="OQO, NNW OW OOK UN 
J 


on 
* 


~+-+ee++ + 


* + 


>= 


Ofte Feuseseeuuv SH SH 


uw > 


uw Ww 
eee & © 


wewet ene uve eu 
* * * 


a 


MB wet Wwnwrs 


as uu & 
** 


wn 





W. D. Heintz Vol. 116 


uv - 
MW 
aa © 
ON r 
.- & 


- Ww Ww WN 


13175 


- 


13265 


wuuUMuununm s 
oor OO KF WY 
wurk SUM WwW SW 


oo 


13493 
13506 


> 


4 
e 
oe 

- 


I 3800 
13357 
13961 
14056 


157°0 


round 


Nnawny 
2a Ui = 


round 
36°2 
114° 
205°3 
4 


355° 


6 
° 


14233 
14250 


14300 4 Aq 


aM 
- WN 


2 


3°5 
63:1 
14499 p3 7: « Eql AB 5 287 


9 
AC 57 607°9 


wn 
w 


ee oases 


14573 139°: 
139° 
194° 
193°S 


a) 


14592 


146045 

14066 527 236 
1626 153° 

14759 fp 270 18° 

14773 535 4 Eql 23°. 
5258 # Ind 274° 
766 6 Mic 270° 

15007 , 2799 

15176 1212 24 Aqr 
§278 A Oct 

15313 ) 2825 

19 

15371 4600 

15447 / 75 

15530 fp 276 


mwwt 


aut UN 
o 4 


99°: 


A Wat O 


nN 


~@OGunenwen ev eee. a 


= 
— 


15596 ; 2854 
5319 
15731 ¢ 2292 
15753 5 339 «641 Aar 
134 
1sgo2z pi 172 51 Aqr 
15934 345 53 Aqr 
15971 , 2909 _«¢ Aq 


uuwtat fe eu 


* 


* 


15975 474 AB 
15985 2912 37 Peg 


lt+tww wre ee 
- One kt OOM ON SN 


was] = 





No. 3, 
ADS 


16131 
16145 
16208 


16215 
16270 


163605 
16425 


16613 
16643 


16649 
16650 
10791 
160350 


10979 


17107 


golg 


gis2 


9909 


9913 
998 2 


h 
p 
p 


Od 


Je 
p 
Hu 
p 
Hu 


479 
2928 
2935 


289 
2944 


340 
5307 
178 
772 
453 


20 
714 
399 

79 
400 


Micrometer observations of double stars 


AC 
y PsA 


Oo PsA 


52 Peg 


? Gru 


_ 


1900 + 


55°06 
65 
4°55 
55°93 
605 
66 

65 

07 

63 
64 

67 


/ 


. 
67 


07 
67 
67 


P 


Ww 


were wustet - HU 


~ 


Stone 59 
Hu 599 
A 492 
I 693 
Rmk 27 
Sh 356 
Sir 14 
‘ 


“ 3046 
I 1477 


C Boo 


p 309 


p> 1883 
h 4707 
OL 288 
y Lup 


u'. Nor 
€— Sco 
A 264 
fh Sco 


>» 2026 


? Phe 


NN ON 


wn uw vv 


67 


www 
> +> 


66 


107 Aqr 65 


—_ 
Ww 


. 


Jl 
30 « 


67 


a 
_ 


67 


-+iuwuteteeeuustuut wv 
* - 


57 67 


Remarks: Comparisons with orbits. 
1 *5, —0°°04 (van den Bos 1936). 
Obs. irreconcilable, evidently spurious 
The orbit Arend (1949) gives 328 °7, 1 
o'or (van den Bos 1936) 
o'os (Strand 1937) 


22, Baize (1954) 282°1, 0°92. 
ob, 


009 (Heintz 1955), —1°0, 
0:06 (van den Bos 1935) 
0:07 (Rabe 1951, Elem. 1), » 0°2, § 0°16 (Ziller 
0°13 (Arend). 

0°03 (Hopmann 1945) 


O'l4 (Florsch 1954 


1954, 


0°40 (Finsen 1926) 

The single observations are : 
55°401 36° 

‘450 37'5 
"453 

“494 
5°s, 


5°93 ‘50 37°3 66 
‘ Or re) 35°2 41 
37°60 5°49 505 z51 41 
35°2 63 617 39°3 5°25 
0°05 (van den Bos 1937), —0°9, | 0°06 (Wierzbinski 1954). 
‘The separation according to van den Bos’ orbit (1950) is 0”:56, but 
the recent elements of Finsen give O'17. 
o'r, 0°00 (Heintz 1955) 
1*4, 0:03 (Woolley/Mason 1945) 
0°7, | 0°13 (Heintz 1955) 
The orbit Dawson 1919 is astray 
work. 
o's, 


A new determination is under 


0°08 (Woolley/Mason 1945) 
3°9, 0°14 (Baize 1942). 
5°7, —0°02 (Woolley 1948). 

Close companion not seen. 
2°3, + 0°06 (Woolley/Symms 1937). 
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10057. A Oph 3°6 0°06 (Rabe 1948). 

10374 4 Oph 4°38, — 0°14 (Baize 1950). 

brs 13 0°7, ~ 0°13 (Wieth-Knudsen 1954). 

py 2173 5°7, | 0°15 (Duncombe/Ashbrook 1952) 
3125 Probably spurious, companion too faint. 


10595 
10855 
11005 7 Oph 1°0, + 0°13 (Hertsprung 1942) 
11046 70 Oph 0°8, o-oo (Strand 1952). 


~ 
5 
~ 


h 5014 0°4, ~ 0°14 (Dawson 1926). 
73 Oph 9°7, ~ 0°13 (van den Bos 1937). 
11186 & 2 o'6, — 0°03 (Wilson 1935) 
11483 >» 1954: 1°5, —~0°16 
1955: 04, ~ 0°03 (Heintz 1955). 
11484 OL 357 Separation 0”:25 according to Florsch 1955. 
11950 C Ser 16, + 0°05 (Finsen 1937). 
11989 HN 1%, - o-os5 (Gottlieb 1948). 
2°0, ~ 0°16 (Dawson 1924) 
8-4, ~ 0°08 (Finsen 1935). 
13493 One obs. : 1955°50, 107°'9, 0” 6, but single at four subsequent good 
nights. 
143600 1954: 50, ~ 0°07 
1955: ~§ 2°7, | 0°08 (Baize 1943). 
14499 0’9, 0°00 (van den Bos 1933) 
14773 , | 0°04 (Luyten/Ebbighausen 1934) 
15176 ‘5, » 0°03 (Danjon 1942). 
1§971 ; 3°38, o'ls 
2°7, 0°20 (Rabe 1954) 
15988 37 Peg , 0°03 (van Biesbroeck 1927). 
16850 A 492 »-7, ~ 0°05 (Woolley/Mason 1948). 
Shr 14 , ~ 0°05 (Woolley 1948). 


Titil 


2. Double Image Micrometer Observations.-.On invitation by Professor 
Woolley, some observations were made in September 1955 with a newly 
acquired Muller type double image micrometer at the same telescope. Muller’s 
formula for the angular screw value, evaluated by Mr K. Gottlieb, yielded 
o”-971 +":002; from 5 determinations using circumpolar stars, 0”:963 + ":027, 
was obtained and adopted. 

23 couples were selected from the previously observed programme in such 
a way as to give the widest possible range in brightness, magnitude difference, 
and separation. Usually, observations were made in two positions: having 
the four images in a straight line and in the vertices of a square, the latter was 
viewed from two sides (by turning the head go’). ‘lhe close couples 6 Eql and 
> 2825, were measurable only in linear position, 6 PsA (very wide, Am=6) only 
in the square. The linear arrangement generally seemed more easily to be 
judged, except for distance measures of close pairs. From 19 couples, the mean 
difference square minus line, is +0 +4, +0”:06, which is not significant. 

‘The differences against the filar micrometer observations for various brightness 
and separation are as follows: 

Square minus Filar Linear minus Filar 


: Ol! (2) ‘ 000 (2) 
(10) q 0°03 (11) 
(6) , 0°05 (4) 
(2) . fohrers) (4) 
(5) : 0°06 (5) 
(7) 0°02 (38) 
(9) , 0'07 (3) 
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Position angles and separations are both measured smaller (the opposite result 
for bright stars being caused by ( Sgr only). ‘Thus, according to the foregoing, 
the double image micrometer decreased the personal equations. 

‘The accidental errors of both series of measurements cannot be fairly compared. 
The mean error of one double-image observation is throughout slightly larger 
than for the filar micrometer, but this is certainly due to the fact that the 115 
observations made are insufficient to gain much experience with the instrument, 
and that they are made with a magnifying power of about 300 (against 560 before) 
and under almost continuously unfavourable seeing conditions. ‘Thus it is 
preferable to list the double-image results separately : 
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It may be concluded that the Muller double image micrometer should not 
replace the filar micrometer but be employed together with it. ‘lhe double image 
instrument yields less observations in a given time than the filar, but gives higher 
accuracy, especially at distances of, say, 1 to 2°5 times the separation limit of the 
telescope, over which range the filar measurements are frequently subject to 
errors. ‘lhe loss of light by the division in two images is almost compensated 
by the fact that no field illumination is required. 


Mumech Observatory : 
1955 December 14 
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MEASURES OF DOUBLE STARS AND THEIR 
SYSTEMATIC ERRORS 


I’. Schmeidler 
(Communicated by R. v. d. R. Woolley 


(Received 1956 March 10) 


Summary 


Measures of double stars made at Mt Stromlo by a filar micrometer and 
by a double image micrometer of the Muller type are discussed and published 
A systematic difference of 1 or 2” in position angle between both micrometers 
is found. A comparison of filar micrometer measures made by the author 
at Munich in 1938/39 with modern orbits shows that the position angles 
are systematically too great for pairs with a separation smaller than 2”. It is 
concluded that double image micrometer observations are less subject to 

ystematic errors than filar micrometer observations. 


1. Purpose and programme of the observations.—During the author’s stay 
at the Commonwealth Observatory on Mt Stromlo, Canberra, measures of visual 
double stars were made with the g-inch Oddie refractor in September 1955. 
‘The main purpose of these observations was to gain experiences with the newly 
constructed double image micrometer of the Muller type rather than to make 
a contribution to the number of existing series of double star observations. 

‘The double image micrometer was constructed by the workshop of the 
Commonwealth Observatory, the prism having been sent to Canberra from 
France. Unfortunately, the only eyepiece available for the micrometer had 
a comparatively low power of about 300 times. ‘lhe construction of an eyepiece 
with higher power was not completed before the author left Mt Stromlo and 
therefore the present double image micrometer measures are not as accurate 
as they could be. 

A number of double stars within the greatest possible range as to distances, 
magnitudes and magnitude differences were selected as programme. Altogether 
23 pairs were observed for which 45 measures were made by a filar micrometer 
and 63 measures by the double image micrometer. Observations by filar micro- 
meter were included in the programme, as it was intended to test the results for 
systematic deviations. 

‘The filar micrometer was the same as that used by W. D. Heintz (1); its 
screw value 29"°718 as determined by Heintz was adopted. For the double 
image micrometer the theoretical screw value 0":g71 was checked by two indepen- 
dent methods. ‘The first method consisted in estimating the time interval between 
the passages of both images of any star of high southern declination through the 
same wire in the field. ‘lhe second method compared distances of wide pairs 


as measured by the double image micrometer with measures of the same pairs 
by the filar micrometer. ‘lhe following three results were found for the screw 
value of the double image micrometer 
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R=0"-971 + 0"-002 as theoretical value 
1”-010 + 0"-017 (m.e.) found from time intervals 
0-938 + 0"-019 (m.e.) found by comparison with filar micrometer 
measures of wide distances 

No plausible explanation for the discrepancy between the two empirical values 
could be found, although the difference is greater than the mean errors would 
allow. As the mean of the two empirical values agrees very well with the 
theoretical value, the latter was definitely used. 

The measures by the filar micrometer were made with the eyes parallel or 
perpendicular to the line joining both stars, whichever position was more 
convenient. ‘I'wo different methods were used with the double image micrometer, 
the four images being arranged either in a square or in a straight line. In most 
cases, the same pairs were observed in both ways at the same night successively. 
‘The author believes that measures with the four images in a straight line are a little 
more convenient and more accurate, but the number of observations is too small 
for a definite judgment of this question. 

2. Systematic differences.—F¥or all stars for which double image micrometer 
measures were made in both ways, the difference between the results obtained 
with a square position of the four images (Dq) and those with a linear position 
(Dl) was formed. ‘The mean difference in the sense Dg minus D1 is 

Ap= —0'+55 +0 *37 (m.e.) 

Ap = +0"-02 + 0-04 (m.e.) (13 stars) 
As no difference significantly exceeding the mean error was found, all observations 
made by the double image micrometer by either method were combined in mean 
values (D) and compared with the mean values of the filar micrometer results (F’) 
for the same stars. ‘The mean difference F minus D is 

Ap = +2°°37+0-70 (m.e.) 

Ap = —0"-06 + 0-03 (m.e.) (15 stars) 


The systematic difference in p is small enough to be neglected, but there is a 
remarkable amount Ap in the sense that all position angles measured by the 
filar micrometer are about 2° too great. 

The question arises which of the two micrometers gives the wrong results. 
A comparison with the results of other observers strongly supported the assump- 
tion that the present author measures positon angles systematically too great 
with the filar micrometer. Some double stars of this programme were also 
observed by Rabe at Munich and by Heintz at Mt Stromlo(1)in 1955; the average 
differences between the present filar micrometer measures and the results of 
Rabe and Heintz are 


Schmeidler—Rabe Ap = +2°'0+0°62(m.e.) Ap 0”-02 + 002 (m.e.) 
Schmeidler- Heintz + 1 6+0 -g2 (m.e.) 0”-05 + 0:02 (m.e.) 


which again indicate a systematic error in position angle. 


In order to examine the question more carefully, the observations made by 
the author in 1938/39 at Munich (2) with a filar micrometer were compared with 
modern orbits. For 48 of the 224 couples contained in the above-mentioned 
publication orbits were found which agreed satisfactorily with all observations 
until at least 1940. Any residuals of individual observers against these orbits 
can therefore be considered as differences between the corresponding observers 


'g 
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and the “‘ average’’ observer, at least for observations before 1940. The mean 


values for the 48 pairs were 


Ap=+0-g0 +075 (m.e.) 
Ap= +0”-028 + 0”-013 (m.e.) 
in the sense Schmeidler minus average observer. ‘The value of Ap is neglegible, 
whereas Ap is again positive, although not much greater than its mean error. 
If the material is subdivided according to distances, the following results are found 

pat + 1°O Ap= +0"-10 18 stars 

p=1-2" +12 0”-00 1g stars 

2-4" o°'5 -0"-06 5 stars 

p>4" +08 0”:05 6 stars 
which are a clear indication that the position angles of pairs with a separation 
smaller than 2” are measured 1 or 2° too great, whereas the observations of the 
wider pairs are systematically right. A similar result holds for the distances 
which are generally too great for pairs closer than 1” and practically correct for 

wider pairs. 

An additional confirmation that the measures made by the author in 1938/39 
at Munich are systematically correct for pairs with a separation greater than 
2” can be found by comparison with photographic observations made by 
Hertzsprung (3) and Strand (4). After corrections for orbital motion had been 
applied wherever necessary, the following differences were found 

Schmeidler- Hertzsprung Ap=+0''2 Ap o”-O1 


22 pairs 
Schmeidler— Strand + O'2 


0°02 26 pairs 
which are entirely neglegible. ‘he distances of these pairs are, with very few 
exceptions, greater than 2”. 

3. Conclusions,—\t seems that the systematic difference of about + 2° between 
the author’s filar micrometer and double image micrometer measures of position 
angles is almost exclusively due to systematic errors of the filar micrometer 
observations. As the residuals of the Munich observations against modern orbits 
show a systematic error in position angle which is not much greater than 1°, 
it is possible that the double image micrometer measures also are subject to a small 
(negative) systematic error in position angle. ‘lhe author is, however, inclined 
to believe that the mentioned discrepancy is due to the uncertainty of the used data 
and to consider the position angles as measured by the double image micrometer 
as systematically right. 


It appears that double image micrometer observations are systematically 


more reliable than filar micrometer measures, specially for close pairs. ‘The 
present discussion shows that it is still interesting to study the systematic errors 
of double star observations, a problem which has been rather neglected in the 
last decades. 

4. The observations.—All individual measures are published. For each pair 
the usual designation is given and, if existing, also the ADS number. In the 
next line, the coordinates for 1950 and the magnitudes are given. For each 
observation, the date, the measured position angle and the measured distance are 
given; / means observations by the filar micrometer, Dg and D/ observations by 
the double image micrometer with the four images in a square or in a straight line 
respectively. For the F observations, . 


. means eyes parallel and : perpendicular 
to the stars, 
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Equal weight was given to all observations in combining the mean values. 
No systematic correction has been applied to the filar micrometer position angles, 
since the preceding discussion gives no clear result about the exact value of the 
error. As for most pairs the double image micrometer observations outnumber 
the filar micrometer measures, the mean of all observations should not be affected 
much. 

‘The author wants to thank Professor Woolley for allowing him to work with 
the Oddie refractor. Acknowledgments are also made to Professor Rabe and _ to 
Dr Baize for communication of residuals against unpublished or unavailable orbits, 
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Residuals against existing orbits were formed wherever possible. ‘The 
theoretical values of position angle and distance were taken from Muller’s 
catalogue of ephemerides (5). ‘The following residuals in the sense O-C were 
found : 


& 
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Finsen 1938 

Baize 1943 

van den Bos 1933 
Luyten 1934 

Rabe 
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wu 
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©20900000, 


van den Bos 1934 


Sternwarte Manchen, 
Munich, Germany: 
1956 February 28 
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THE ORBIT OF ‘THE SPECTROSCOPIC BINARY HD 170523 


William Buscombe 


(Received 1956 February 11) 


Summary) 


Radial velocities have been measured from twenty-four new spectrograms of 


the giant L-type star 0* Telescopii. ‘The orbit has been derived by combining 


these with eleven observations made in Chile during 1908-17 and three in 


South Africa during 1952 
The mean velocity of the system is —7°6 km/sec and the velocity amplitude 


of the primary is 34°9 km/sec, with a minimum semi-major axis of 10°2 * 10° 


km. ‘I'he period is 219°70564, eccentricity 0°218 and mass-function o’o8g1 
solar masses. 


Introduction and observations.—A three-prism spectrograph has been in use 
since early in 1954 at the Cassegrain focus of the 30-inch (Reynolds) reflector 
at the Commonwealth Observatory. ‘lhe first principal programme is the 
determination of radial velocities of stars in the Scorpio-Centaurus cluster. 
At the outset, the observing time has been concentrated on stars brighter than 
visual magnitude 6-0 with uncertain velocities. 

The star HD 170523 (6?'Tel), R.A. 18" 28-3, Dec. —45° 48’ (1950), vis. 
mag. 5°3, Was announced as a spectroscopic binary from two plates obtained 
in 1908 by the D. O. Mills Expedition of the Lick Observatory at Santiago, 
Chile (1). Additional spectrograms were secured there between 1914 and 1917 
(2), at Radcliffe Observatory in 1952 (3), and at Mount Stromlo in 1954, but 
the period of velocity variation remained unknown. During 1955 a determined 
effort yielded 22 additional spectra from which an orbit could be computed. 

The observed radial velocities are listed in ‘lable I. Column (1) shows the 
date and Universal ‘Time. In column 2 is the phase in days, calculated on the 
basis of the finally adopted period and an arbitrary origin JD 2435 200-0. ‘The 
initials in column 3 indicate the observer and spectrograph used, for which a key 
is given in ‘lable Il. ‘The numeral shows the number of prisms. Column 4 
shows the mean plate velocity, reduced to the Sun, and column 5 the residual 
(observed — computed) from the final orbit in km/sec. 

The design and operation of the spectrograph on Mount Stromlo will be 
described subsequently. ‘The straight slit employed is 60 microns wide, corre- 
sponding to 15 » projected on the plate. ‘lhe star image is transmitted through a 
rocking plane-parallel glass plate to produce spectra of uniform density across a 
width of o-6 mm at the slit. Eastman Kodak emulsion 103a-O was used 
throughout, developed for 5 minutes in D-1g. Experiments on the distortion of 
emulsions (4) led to the adoption of Kodak ‘Tropical Developer after 1955 
August 15. Also at this time the exposure time could be reduced from 12 
minutes to 6, due to the higher reflectance of freshly aluminized mirrors. 
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TABLE I 


Radial velocity observations of Hb 170523 


Date (U.T. Phase Observer Velocity Residual Remarks 
d 
July 5°353 » L 31 
Sept. ‘10 $22 : 12 1! 
July 31° 003 ‘ +28 
July . ‘078 
7065 


: ‘0738 
Aug 4° 302 


302 
‘ 262 ’ 
June 5° ‘293 : 
Sept. 564 
May ‘ 8173 MF, 
ok86 Mi 
Sept. "7 315 MF, 
Aug. , ‘O99 AP, 
2°770 WH, St: 2 ») Weak comparison 
Sept. , 0°290 MIF, 
333 MF, 
May 523 HiG, 
June 742 WB, 
QQ! WB, 
‘000 WB 
5 OO! HG, 
‘244 WR, 
‘230 WB 
182 WB, | ‘ underexposed 
456 AP, ! fogged 
2449 WB, !{ : ; underexposed 
Sil WH, 5 15 underexposed 
5 °406 GH, St3 10 , underexposed 
‘417 AP, St3 35 
ye SY, WBE, St3 28 
693 WB, St3 32 4 overexposed 
‘O10 GH, St3 34 
983 WB, St3 38 
003 HG, St3 36 2 weak comparison 
‘425 ‘936 GH, St3 43 
10°408 $969 WH, St3 24 
11°390 9°95! WB, St3 21: 4 underexposed 
29°429 6°284 WB, St3 35 
*'These additional observations, communicated privately through the kindness of 
Dr A. D. Thackeray, are not included in Fig. 1 


TABLE Il 


Observers and Spectrographs 

Obs. Name Inst. Location Dispersion at Hy 
WB W. Buscombe Li Lick Obs., Santiago, Chil 37 A/mm. 
MF M. W. Feast L2 Lick Obs., Santiago, Chile 

HG H. Gollnow 

GH (;. Hagemann ni 
WH W. D. Heintz R2« Radcliffe Obs., Pretoria, S. Africa 


AP A. Przybylski St3 Mount Stromlo, Canberra, Australia 


Rab Radcliffe Obs., Pretoria, 5S. Africa 
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The plates were measured on a Hilger long-screw engine with a microscope 
of 30-power. ‘The reductions are based on the wave-length system developed by 
Petrie (5), with corrections for curvature of the iron-arc comparison lines. 
Standard-velocity stars are observed nightly; the velocities on 20 plates taken 
during the period, covered by the observations of the binary star, never departed 
by more than 2 km/sec from the corresponding values adopted by the International 
Astronomical Union. On each spectrogram of 6?'Tel, from seven to twelve 
absorption lines were measured, depending on quality. ‘The internal probable 
error of the plate mean velocity is in each case between 3 and 4km/sec, which 
corresponds to the scatter from the final orbit. However for the weak plates, 
for which the velocity is followed by a colon (:) in ‘Table I, the probable errors, 
due to the scatter of measures from diffuse lines, reach 6 to g km/sec. 

The derivation of the elements.—'rom the observations obtained in the winter 
of 1955, the period was found to be approximately 21 days. Since a number of 
observations had already been secured 40 years previously, it was possible to 
determine a provisional period of 21°-707 graphically by trial and error. However 
it is seen from Fig. 1 that those early observations fall almost entirely on the 
rising branch of the velocity curve, and do not delineate the extrema very 
precisely. ‘lhe trial period was first adopted, and a graphical solution made 
by the method of Lehmann-Filhés. ‘These elements were then adjusted by a 


= oe 
Velox ity HD 170523 | 


e '908-17 
> + 40 km/sec + 1952-55 


Weak plates 
> +20 


> Oo 

+ 

> -2 ~e 
(+) 

> - 40 


Oo 
- \ 











The velocity curve of HD 170523, drawn according to the final elements. The origin of 
abscissae is taken arlntarily at J.D. 2,435,200°0. 


first least-squares solution based only on the individual observations made in 
1955. 

‘To facilitate the computation with a desk electric calculator, the equations 
of condition were expressed in terms of the Thiele-Innes rectangular 
coordinates 


y+LdX/dM+NdyY/dM, 
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and 
AV =Ay+ALdX/dM+ANdY/dM + Pde + OnAT + R,An, 


0 {dX Ody 
"os (3a) ad (su ) 


0 (dX 0 {dX 
=~ hem (su) aad) (sm) 
=-(t-T)Q,; 


L=1-92 = 10“asinisin w, 
N= 1°92 x 10 4a sini cos w. 


where 


‘The period was held unaltered. 

The theory is given by van den Bos (6), who has also published tables of 
and Y for each degree of the mean anomaly and each step of o-o1 in the 
eccentricity (7). Values of the derivatives (dX/dM and dY/dM) for any 
eccentricity can be interpolated easily, and the second partial derivatives with 
respect to eccentricity and anomaly obtained by differencing. 

As the early Lick observations showed a mean residual of +6km/sec from 
this orbit, it became clear that a slight decrease of the period was necessary. 
‘The amount was derived from the slope of the ascending branch of the velocity 
curve. For computing the final elements, eight normal points were adopted, 
representing all the observations in groups of from two to five within close 
intervals of phase. ‘lhe unweighted normal equations were solved again by 
least-squares to yield the final elements shown in ‘Table II]. The individual 


residuals are shown in ‘lable I. On many of the plates the density of the 
continuous background is not sufficient to permit a measurement of the velocity 
from the K line of Cat, but from the measures shown it appears likely that the 
absorption feature is mainly stellar. 


Taste Lil 
Orintal elements 


Preliminary Adjusted Final, with Probable Error 
P 21°707 days 21°707 21°70564 + 000005 
y 10°0 km/sec 13°3 7°6 o'ls 
kK, 40°0 km/sec 34'0 34°9 0°3 
L + 238°4 km/sec 340°2 429°9 8-7 
N 2171°0 km/sec 1756°6 1915°6 30°2 
T-J.D.2435200 16°60 days 17°0958 16-6688 0°055 
e 0°*302 0°401 o°218 0°023 
w 6°26 10°96 12°65 O's 
a, sini 11°4 * 10° km 9°3 - 10°2 00% 

(0°068 a.u.) 

f(m) o’o%gI1(-) 0'0021 


Discussion.—By comparison with spectrograms exposed with the same 
equipment of standard stars in the MK system (8), the spectral class of 
HD 170523 is approximately B6 11-111. The hydrogen lines are always more 
diffuse at the phase six days, i.e. at maximum velocity of approach. Blaauw (9) 
lists it as a doubtful member of the Scorpio-Centaurus cluster, and in fact at 
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the probable distance of 1400 parsecs, this star would be far beyond the utmost 
extent of the cluster, whose known members are almost exclusively main-sequence 
stars. 

The absolute (visual) magnitude of — 3-5: corresponds with that derived 
statistically from the assumption that the star’s proper motion is purely a 
reflection of the solar motion. ‘The corresponding mass of the primary from 
the mass-luminosity relation would be about 16 times the Sun’s. Since no 
trace of the companion can be seen visually or spectroscopically, it must be at 
least a magnitude fainter, with a maximum mass of 8©. If in the equation 
f(m)=0-08g1” we insert the minimum value of m,/m,=2, we find that the 
least possible inclination is i= 29 ', so that the maximum length of the semi-major 
axis 18 4;=21 ~10°km (o-1ga.u.). A smaller value of the inclination could 
only occur if both components depart from the mass-luminosity relation. 

On the other hand, any larger value of i 29° could be made consistent with 
the primary’s mass and luminosity if we are willing to accept a larger value of 
the mass-ratio (i.e. a greater difference in magnitude of the two stars). ‘The 
maximum value of the mass-ratio, in the case where 1=go", is m,/m,=5. If the 
inclination is high, and the companion an A-type giant or main-sequence star, 
the chance of detecting it even from infra-red spectra would be very slight. 

Detailed photometric studies have not been made, which would exclude the 
possibility of a shallow eclipse. However, from the absence of a deep eclipse, 
it is reasonable to infer that the inclination is not near go’, as the radius of the 
primary component forms a large fraction of the radius of the relative orbit. 
Hence it appears probable that the mass of the secondary component is between 
four and eight times the Sun’s. 

Icknowledgments.—I\t is a pleasure to thank Professor Woolley and Dr Eggen 
for advice on the Thiele-Innes method, of which no previous application to 
exclusively spectrographic observations seems to have been published. I am 
also grateful for the helpful interest of Dr Gollnow throughout the work and 
for the share taken in the observing by Drs Hagemann, Heintz (of Munich) 
and Przybylski. 


Commonwealth Observatory, 
Canberra: 
1956 January 12. 
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STELLAR PARALLAXES DETERMINED AT MILL HILL 
(THIRD LIST) 


University of London Observatory 
(Communicated by the Director) 


(Received 1955 December 28) 


The following is the third list of relative parallaxes determined with the 
Radclitfe 24-inch photographic refractor; the earlier lists were published in (1), 
(2). In the second list there was a mis-identification pointed out to us by the 
Naval Observatory at San Fernando, viz. the star No. 15 on that list was A.C. 

6 3939-98 and not A.C’. —6 3939-101. 

‘There has been a noticeable deterioration in the quality of images on plates 
taken in recent years which is probably due to a small relative movement of the two 
components of the Radcliffe 24-inchlens. [tis thought that this has had an adverse 
effect on the accuracy of the more recent parallax determinations and for this reason 
the measurements are being terminated with the list published in this paper. 
‘The mounting and performance of the lens is now under investigation, 

All stars in the present list except No. 15 were taken from our first parallax 
programme. ‘The results for two of the stars (Nos. 15 and 17),were obtained by 
measuring the plates with a recording microphotometer, as described by Kiang in 


(3). No. 17 was also measured with our usual measuring microscope, and the 


result, 7 0"'003 + o”"-022 was obtained. Ihe very large probable errors 
obtained in both measurements were due, at least in part, to the fact that the 
comparison stars were very unsymmetrically placed; they were situated at three 
corners of a square surrounding the star. 

All stars were photographed for at least six epochs. More specifically Nos. 11, 
16, 17 were photographed for six epochs, Nos. 2, 3, 4, 5, 6, 8, 10, 12, 14 for seven, 
Nos. 1, 7, 13, 1§ for eight, and No. 9 for nine epochs. ‘lwo exposures were made 
on each plate and an average of 4°5 plates per epoch was obtained. Most of the 
plates were taken between 1946 and 1951 inclusive. Nos, 1, 13, 16 each had an 
early epoch in 1939, and Nos. 15, 17 each had two late epochs in 1953 and 1954. 

Details of exposures, types of plates, and reference plates were the same as 
before (2). 

‘The photography was carried out by Mr C. C. L. Gregory, Dr E. M. Burbidge, 
Miss A. C. Robinson, Dr G. R. Burbidge, Mr R. W. Pring, Mr 5S. K. Wang, M1 
'T. Kiang and Professor C. W. Allen. ‘Twelve stars were measured by Mr 5. K. 
Wang and Mr ‘I. Kiang, the remaining five by Mrs G. de Vaucouleurs, Miss 
A. C. Robinson, and Mr C. C. L. Gregory. 

In the list, the second column gives the designation of each star. Radcliffe 
numbers give the number of the selected area followed by the Durchmusterung 
number, and are taken from the list of stars with proper motions exceeding 0”: 100 
in the Radcliffe Catalogue of Proper Motions, p. xl; Yale numbers refer to the 
1935 General Catalogue of Stellar Parallaxes. ‘The third column gives the 
approximate visual and photographic magnitudes. ‘hese, and also the spectral 
types, are taken from various sources and should be regarded as provisional only. 
‘l'wo values of the annual proper motion in R.A. are given, the first being a catalogue 
value and absolute (as are also the proper motions in declination ;) the second value 
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is derived from the parallax plates and is relative; the average difference between 
this and the catalogue value is +0":02z0. The probable errors of unit exposure 
were divided by V 2 to give the probable errors of unit plate, which are shown in 


the eighth column. In the last column the upper value refers to the parallax and 
the lower to the proper motion in R.A. 


Unit — 0” oo! 


-_* 


Designation Mag. R.A. Proper P.E. 


and Vis. Spec. Dee. Motion Unit Parallax P.E. 
B.D.or A.C. No. Phot. Type — 1goo’o 





a a | 
R.A. Dec. 
hm —°5s 

10 §4 156 

15° 43 


39 57 1 113 
44 574 45 11 


Radel. 68°354 G3d 
15 27 


Radcl. 22°752 G 


3 06 30 
A.C. +79° 1584 . 79° 34 
Camb. IIIc, 67 3 25 14 
43° 744 43 20 
Yale 1364 


19 55 
15 621 


15° 43° 
54 3° 
39° 49’ 
Radel. 75 12 57 
15° 1535 g° 12 
McC, 120°25 


Yale 1567 
39° 1138 


42 54 

A.C. + 29° 26176 ‘ 28° 53° + 78 
43 10 + 2 

44 1910 44°46 — 1 

Yale 3526 

68° 652 

SW Dra 


Il 32 12 10 
67° 53 3 
12 12 48 40 
+7o°4' +S 
13 03 52 196 
390° 237! +29" 55° 265 
Yale 4495 15 03 06 826 
25° 2874 +25° 18’ -891 
Yale 5071 16 55 29 +247 
65° 1157 5° +65° 17’ 220 
Yale 6323 20 28 00 )=—+- 300 
32° 3853 +32° 40° +301 
Camb. VIIa, 7 K2 21 58 56 +150 
11°8 59° 06’ +151 
HD 5223 8:8 R3 00 48 54 +135 ° 
23° 123 +23°32' + 93 
* No previous determination of parallax. 
| Measured with a microphotometer 
Umiversity of London Observatory, 
Mill Hill Park: 
1955 December 21. 


Radel. §7°192 
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[IONIZED AND NEUTRAL CALCIUM IN THE ATMOSPHERES 
OF THE CEPHEIDS » AQUILAE AND ¢€¢ GEMINORUM 


W. N. Abbott 
(Communicated by F. J. M. Stratton) 
(Received 1955 November 7)* 


Summary 

The equivalent widths of the ionized calcium H and K doublet and of the 
neutral calcium singlet at A 4227 in the spectra of these two Cepheids have 
been measured in relation to the phase of the light variation. 


‘To render this 
possible, certain simplifying assumptions were necessary. 


Abundances of 
Cat and Ca I as well as ionization temperatures in the atmospheres of these 
variables in relation to phase are derived. Gas-pressures permit, through the 
application of the perfect-gas formula, the computation of the values of the 
specific volume of the calcium atmosphere which is indicative of the general 
behaviour of the atmospheres of these stars in the course of their cycle. 
Results show a difference of phase existing between the general, photospheric, 
pulsation and the volume variations of the atmospheres of the Cepheids. 


1. Material.—The spectrograms for this investigation of the variations of 
the equivalent widths, W,, of the ionized calcium lines H and K, and of the neutral 
calcium line 44227 (W,), in the spectra of 7 Aql and ¢Gem were as follows. 
In the case of 7 Aql spectrograms were obtained at Ann Arbor with the 37-inch 
reflector, equipped with the Cassegrain two-prism spectrograph (24-inch camera), 
of the University of Michigan giving a dispersion of 12°8 A/mm at 4000 A. ‘The 
plates used were Ila-O Kodak, and exposures depended on phase, the region of the 
spectrum, and seeing. ‘The slit was open at o-o4mm, and the plates were 
calibrated with the calibration spectrograph of the Observatory recently adjusted 
and equipped with a double-wedge slit. ‘he spectrograms were very carefully 
developed, tracings were made with the MacMath—Hulbert deflection micro- 
photometer, and results were, in general, satisfactory, although occasional poor 
seeing due to the position of the star, early evening hours and local conditions 
constituted a handicap. ‘The spectrograms of ¢Gem, on the other hand, are 
Mt Wilson plates of a higher quality in all respects, taken with the coudé grating- 
spectrograph of the 100-inch at a dispersion of 2-9 A/mm. We took tracings of 
these with the Mt Wilson deflection microphotometer at Pasadena; after careful 
examination a number of plates of both groups were discarded as unsuitable, 
and those used, eight for 7 Aql and ten for ¢ Gem, are given in ‘able I. 

The equivalent widths W, (the sum of the equivalent widths of the ionized 
calcium H and K lines) and W,, (of the neutral calcium line A 4227) were measured 
on these plates with the object of acquiring the maximum possible information 
on the respective abundances of these calcium ions and atoms, as well as on other 


physical data, such as ionization temperatures, gas pressures, and specific volumes 
in the calcium atmospheres in relation to phase. 


* Received in original form 1955 May ig 
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‘TABLE | 
7 Aql ¢ Gem 
Pl. No. Date* Ph. : Pl. No. Date* Ph. 
R Ce 
19686 


19675 


Oct. 10 0°06 2301 Apr. 016 
2 Sept. "19 2263 Feb. *30 
2 Sept 33 2722 Jan. 3¢ "41 
2 Oct. !$ ‘64 i4) 2724 42 Jan. : SI 
Sept. 2 76 2764 Apr. 2 ‘51 
Oct 78 2742 42 Mar 66 
Sept. 2 “gl 2977 43 Mar. 2: ‘07 
2 Oct 3 o'92 cF4p 2980 43 Mar. a 
2707 42 Jan. ‘85 

2988 42 Mar. 28 og! cF gop 


mw Wa 


"- 


196 77 / 


19690 
19659 4 
19696 / 
19664 


~~ we ee oe 
wueuwuuwuuu Vv 


19033 


* The decimal fraction of the date is about o'8 for mid-exposure 


2. The 4227 line.—This resonance singlet is the first member of the principal 
series of the singlet system of the Ca 1 atom and its study presents no difficulty. 
Although small, the variations of W,, were easily measurable. 

3. The Catt doublet.— This resonance doublet, being the first member of the 
principal series, arises in absorption through the transition of this potassium-like 
ion from ground level, 4°5 to 42P. ‘The wave-number interval between the two 
components is 223cem™'. ‘The duplicity of the P level gives these two lines, 
45,0 +47 Poo and 4°5,5->47P,, the shortward one 4?5,)->4°P,,. being more 
intense, and in the doublet energy levels of atomic systems containing but one 
valence electron the fine-structure level, j=1-— 4 generally lies deeper than the 
}. Between levels 475 and 4*P there lies in the calcium 
ion a metastable level, 32D, from which allowed transitions to 47D give, in 


corresponding level) = 1 


absorption, the infra-red triplet A 8662, very pronounced in the Sun and probably 
so in the classical Cepheids if we could explore that region; unfortunately we are 


restricted to the H, K doublet, and information on the Ca 1 1on and its abundance 
may be gained only through measurements of the W values of these lines. 
\ formidable obstacle to such measurements, however, especially to precise 
absolute measurements, constitutes the very heavy blending of this region in these 
spectral classes. Hundreds of metallic lines of atoms and ions crowd this area 
and in assigning the contours of H and K we have made due allowance for the W 
of the blends. ‘lhe systematic errors, important in the case of research on the 
absolute values of W, are in our case unimportant, as we are interested in the 
variations of W as a function of phase and not in the absolute W values themselves. 

4. Assumptions.—Short of abandoning the problem altogether we have 
made the following working assumptions: 

(a) ‘he contours of H and K are assumed to be exact provided they do not 
cut the line of zero intensity. If they do so they are discarded as useless; no 
correction is made on the reduced contours even if they merely touch the zero line. 
‘The effect of instrumental scattering is therefore neglected. In the case of the 
Ann Arbor plates, instrumental scattering cannot but be negligible; the 
Mt Wilson plates are more affected by it and the contours of the ¢Gem lines are 
markedly less deep than the contours of the » Aq] lines. 

(b) As the blends are fine and belong both to atoms and ions, the resulting 
variation of their total W, or XW within the region of H and K as a function of 
phase may be assumed {dW =o for all practical purposes. ‘The dW are of the 
order of + 0-01 A while the variation of H and K is of the order of 10 A. 
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(c) The non-systematic errors made under the above assumptions (a) and 
(+) are no greater than, but are of the order of, errors resulting from measurements 
on curves, the use of a slide rule, plotting, tracing, reducing, and, in general, 
errors of approximation. 

5. The He line.—TVhis blend comes under assumption (4) as the variations of 
its equivalent width cannot exceed +0°03A; and these are probably balanced 
by the contrary variations of other blends. 

6. Continuum.— Placing the continuum may, especially in cases like this, 
be a highly controversial undertaking. Wherever we place it, however, it is 
important to keep it there all along. ‘The highest-peak method gives a continuum 
which is too low, as the actual continuum certainly lies a good deal above the peaks 
ot the microphotometer tracing, its actual location being to all practical purposes 
unknown. ‘The highest-peak method offers nevertheless the advantage of 
supplying a conventional continuum which is solidly anchored. As the spectral 
classes of these stars are G and F we assume that we may advantageously, if not 
decisively, check the continuum thus obtained by comparing it with the continuum 
of the Utrecht Atlas. Line too of this Atlas was taken as the continuous back- 
ground; if A, is our continuum intensity, A the intensity of a point on our 
spectrum, and U, and U are, respectively, the Utrecht continuum intensity and 
our intensity at the same point, we have: 


: | od Idd 


si ‘fue 


The agreement for the two sets of spectrograms is satisfactory. 


1 ("Ad= 2 ('UdA; A, =U, * 


A.J, l 


" H 











Contours of H_ and K in » Aquilae at different phases 
(The vertical lines are metallic lines of reference.) 
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7. Profiles.—The profiles of H and K for the two stars resulting from the 
reduction of the tracings by means of the characteristic curves are shown in 
Figs. 1 and 2. ‘The variation of their contours as a function of phase is very 
similar for both stars and will be discussed more at length in a subsequent 
paragraph. 











Contours of H and K in © Geminorum at different phases. 
(The vertical lines are metallic lines of reference.) 


8. W-values.—TVhese were derived in the usual way: areas of the contours 
were measured with a planimeter and W=S/A, values for H and K and for 
A 4227 are tabulated in ‘lable II and plotted in Fig. 3; the scattering of the points 
is comparatively small and mean curves illustrating the mode of variation of 
W, and W,, may be easily traced. ‘The two different values of W, of ¢Gem at 
phase 0-51 may give an idea of the mean error, provided the differences are not 
real and due to different conditions in two different cycles. ‘The mean error is 
for » Aql + 1°3A and for (Gem+o-8A. In Fig. 3 the light curve of 7» Aq) is 
by Wylie (1) and of €Gem by Harris (2), whom we thank for his collaboration 


on this occasion. ‘lhe radial velocity curves are both by Jacobsen (3). 


g. Abundances.—'\'o derive the abundances of Catt from W, values we used 
the curve of growth for the Schuster-Schwarzschild model and pure radiation 
damping for the line K (3933 A) which Aller (4) has computed, with pure radiation 
damping absorption coefficient a=0'0287 and oscillator strength value 
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alues of A 4227 and of H (dotted), and K (full) with phase for » Aquila 
and A Geminorum (higher curves). 
Middle curves are light curves. 
Lower curves are radial velocity curves. 


f=2/3x 1-19. ‘To simplify things we used the same curve of growth for the 
A 4227 line and the values of log N, and log N,, thus derived are tabulated in 
Table Il (columns 4 and 5). ‘The lower W, values and the resulting lower N, 
values for € Gem are probably real. 


Tasie Il 


Ph. log W; log W, log N; log N, log P, log P, T 5040 V 


deg. K 7 10’ cm" 
7 Aql 
0°06 "565 0'087 19°95 
‘19 624 ‘042 
"33 ‘708 134 
"04 644 ‘170 
‘76 653 ‘179 
‘78 ‘634 ‘170 
‘9! 545 "046 
0'92 ‘519 "152 
¢ Gem 
"16 403 0°409 19°75 
0°30 634 ‘252 20°08 
‘41 544 143 19°88 
‘51 "506 194 19°33 
51 "519 ‘104 19°35 
‘66 ‘532 ‘060 19°90 
‘68 478 ‘222 19°77 
‘TF 449 ‘137 19°73 
85 381 284 19°60 
o'"9l 362 0°244 19°56 
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10. Electron pressures.—-M. and B. Schwarzschild and W. S. Adams have 
given values of the electron pressures in the reversing layer of 7 Aql (§) derived 
by an ingenious if a rather indirect method, by comparing metal abundances 
in the Sun to metal abundances in 7 Aq] and using rather freely the method of 
least squares. Kelative values are given by them of log P,/P, ©. We used their 
relation and by taking P, 1°48 dynecm * we derived P, values for the calcium 
atmosphere which to all practical purposes are the same for the reversing layer. 
‘The same values of P, were assumed for the calcium atmosphere of (Gem. 
‘The calcium lines are probably formed higher in the atmosphere than the neutral 
Ke lines but we may assume that the increase in degree of ionization is approxi- 
mately compensated by the diminution of the density, at least within the practical 
limits of error of this research. Under this assumption the values of log P, are 
given in ‘lable Il, column 6. 

11. Temperatures.—'Vhe Saha formula was used to derive the temperatures 
from the abundances and the electron pressures; thus for each phase of the 
cycle of light variation we have: 


S040 nn ) 

log N, —log N,, + log P, z 7 J + 2:5 log 7 —0-48 + log 2B/B, (J = 6:09 ev) 
‘Lhe ratio of the partition functions was assumed B, & 1°0, 

and we have for the temperatures : 


log N, — log N,, + log P, + 0°28 12250 


a— =log 7. 
2°5 l 
The temperatures were derived by a graphical method and their values as a 
function of phase are tabulated in ‘lable II, column 8. ‘Thus we see that €Gem 
has higher ionization temperatures all along the cycle than 7» Aql. Maximum 
temperature for both stars, in the calcium atmosphere which is formed of the 
mean and higher layers, is at phase 0-15 to 0-20 instead of o-oo, and minimum 
ionization temperature is at phase 0°65 ; this is indicative of a certain independence 
of ionization temperature from surface temperature, as minimum light or 
minimum surface temperature of 7 Aql is at phase 0-7 while for (Gem it is at 
phase o°5, the light curve of this star being almost perfectly symmetrical. ‘The 
temperature curves have a steeper descending branch than the ascending one 
and the respective gradients are : 

yAql AT/At 200 K for At=o-10 phase 

(Gem AT/At= + 226K 
on the descending branch, and: 

» Aql AT/At 145 kK 

€Gem AT/At= +233 


on the ascending. 


As the use of the Saha formula is based on the assumption of thermodynamical 
equilibrium, the ionization temperatures derived therefrom may legitimately 
be considered as being the actual thermodynamical temperatures of the atmosphere 
at these levels, such temperatures being defined by the Maxwellian distribution 
relation : 


m.. T. 


m 
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12. Gas pressures.—By using Aller’s relation (6) between electron pressure, 
temperature, and gas pressure we derive for each phase the gas pressure, P,, 
at the calcium level. Values of log P, in dynecm™ are tabulated in Table II, 
column 7. 

13. Specific volume.—As we have assumed the existence of thermodynamical 
equilibrium we may consistently associate T and P, through the perfect-gas 
formula; and we have for the specific volume of the calcium atmosphere which 
we assume to be representative of the medium-to-high atmosphere of these 
supergiants : 


(A =8-31 x 10° erg mole deg). 


‘These values of the specific volume are contained in column 10 of Table II and 
are plotted in Fig. 4. ‘The full curve shows the variation of the specific volume 
of the calcium atmosphere with phase for 7 Aql, and the broken curve for { Gem. 


5 Gem 





Y Aq! 





a Se a eee 4 
02 )4 0 0.6 0.2 


hic. 4.—Variation of the specific volume of the calcium atmosphere with the phase for 4 Aquilae 
and € Geminorum. 

The vertical dashes locate the phase of maximum volume of the pulsating photo- 
sphere resulting from the radial velocities as deduced from the metallic lines, 
Apart from cycle peculiarities there is a very marked difference of phase between 
maximum star-volume and maximum volume of the calcium atmosphere, the 
dynamics of which may definitely be assumed to be the dynamics of the higher 
atmosphere. ‘This assumption together with the mode of variation V 1s, if not 
conclusive, at least indicative that this difference in phase of maximum volumes, 
photospheric and atmospheric (which is about 7/3-65 for 7 Aql, and 7/4 for 
€ Gem), is real, the maximum atmospheric volume following by these quantiti 
the maximum photospheric volume. 

Working with the metallic lines on the lower atmosphere the Schwarzschilds 
(5) find that, while the pulsation of the photosphere has the characteristics of 
a wave which is neither standing nor progressive but intermediate between the 
two, the pulsation of the higher atmosphere has the characteristics of a progressive 
wave. ‘hese results are interesting in the light of our measurements, 
notwithstanding the doubts which mav exist as to the validity of some of the 
Schwarzschilds’ assumptions, 


20* 
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14. Asymmetries of the lines.—Figs. 1 and 2 show that the profiles of H and K 
present asymmetries slightly after maximum equivalent width, at least in the 
case of 7 Aql. ‘The theory of the influence of a moving atmosphere on the line 
profiles developed by Miss Underhill (7) is not applicable to H and K as this 
would imply velocities of the order of 1,000kmsec™ which is impossible. The 
condition of existence of the star fixes an upper limit to velocities in radial motion 
of the material. ‘lhe velocities cannot exceed for any length of time the critical 
velocity, otherwise gravitational disintegration of the star would be unavoidable. 
There is therefore a relation 

U,<\ 2GM/R 
which links radial motion to the potential; a pulsating star is less stable at 
maximum volume than at minimum volume, since, the potential being lower, 
the condition of disintegration is less remote. 

15. Fluctuations in the position of the H, K lines.—Such fluctuations in the 
position of H and K as compared to the metallic lines as are suggested by Fig. 1 
and 2 where the lines marked M-M’ are metallic lines of reference, are not very 
convincing and may be due to the variation of the depths and intensities of the lines. 

16. Conclusion.—'Vhe variation of the equivalent widths of the Cat and Cat 
lines reveals the variation of the abundances of ionized and neutral calcium in the 
atmospheres of 7 Aql and Gem, two Cepheids of different optical characteristics 
but physically very much alike. By using the Saha formula and assuming 
thermodynamical equilibrium we derive the temperatures in relation to the phase, 
and their variation reveals a certain independence of surface temperature. By 
associating these temperatures with the gas pressures, derived from available 
electron pressures by means of the perfect gas relation, we obtain the specific 
volumes of the calcium atmosphere and its variations related to phase. 
A difference in phase amounting to about 7/3°65 for 7 Aql and to 7/4 for €Gem 
between the volume variation of the atmosphere and the pulsation of the photo- 
sphere and reversing layer is rendered very obvious. ‘Che asymmetries of the 
contours of H and K at maximum W, and certain apparent oscillations of the 
centre of the lines compared to the metallic lines, would be, if real, subjects open 
to further investigation. 

Acknowledgments.—-1 wish to express my thanks to Dr L. Goldberg for 
making available to me the facilities of the Observatory of the University of 
Vichigan, to Dr D. McLaughlin for his assistance, and to Dr I. Bowen and 
Dr W. S. Adams for making available to me the facilities of the Mt Wilson and 
Palomar Observatories in Pasadena. I! likewise thank the Committee of the 
Fulbright-Smith—-Mundt Fund for their grant. 


Athens: 
1955 April 20. 
References 
(1) C. C, Wylie, Ap. 7., §6, 217, 1923. 
(2) D. Harris III, Ap. 7., 118, 346, 1953. 
(3) T. S. Jacobsen, L.O.B., 12, 138, 1926 
(4) L. Aller, Astrophysics, 1953, Ronald Press, p. 289. 
(5) M. and B. Schwarzschild et a/., Ap. 7., 108, 207, 1948 
6) L. Aller, Astrophysics, 1953, Ronald Press, pp. 91, 223. 
(7) A. B. Underhill, Ap. 7., 106, 128, 1947. 





EQUIVALENT WIDTHS OF THE K LINE 
Antoinette de V aucouleurs 
Communicated by the Commonwealth Astronome 
(Received 1955 December 30 


Summar\ 
Determinations of the equivaicnt widths of the A line were made for 
46 stars of spectral types ranging from Ao to M3. ‘The observations are 
compared with the theoretical gradient-—-equivalent width relation for 
radiative models computed by Przybylski. Comparisons are also made 
with the results of other observers and it appears that the equivalent width 
of the K line ts fairly constant for stars later than type Ko 


Determinations of the equivalent widths of the A line were made for 46 stars 
of spectral types ranging from Ao to M3. ‘The plates were those taken with an 
objective-prism telescope of four inches aperture for the monochromatic magni- 
tudes programme of Woolley, Gascoigne and de Vaucouleurs (1). ‘The 
dispersion was 212 A/mm at K. 

Additional material, obtained with the same equipment, was provided 
recently by Dr Heintz; from this were selected late-type stars (6 to M3) 
taken with a longer exposure on 103a-E plates and for which the density in the 
K line region is more adequate. 


The spectra were recorded partly with the Leeds and Northrup micro- 


photometer of the National Standards Laboratory (Sydney) and partly with 
the new Hilger microphotometer now in use at the Observatory. ‘The enlarge- 
ments of the tracings were x 42 and x 61-5 respectively. For each plate at least 


four spectra of a star taken with sectors varying by steps of o-5 magnitude were 
recorded for calibration. 

Most of the densities used for the contours were on the straight part of the 
characteristic curves but for some late-type stars the density at the centre of 
the A line was too small for accurate measurements. ‘l’o this source of error 
is added the usual difficulty of drawing the continuous background and the 
wings of the lines. 

‘The contours were drawn for values of ///, at intervals of 0-03 or 0-04, where 
/ is the intensity of any point along the line and /, the intensity of the continuum. 
‘The two wings were integrated separately at ///, =0°95, 0°35, 0°75, etc., along 
a smooth curve through the points. 

‘lable I gives the name of the star, its MKK spectral type and luminosity 
(or that given by Miss Woods in her catalogue of southern stars (2)), the 
Stromlo or Greenwich gradients (when in parentheses the interpolated value 
corresponding to the spectral type was used), the equivalent width in angstroms 
with its mean error, and logEW. ‘The weights are the number of spectra 
measured (on two plates or more). 

On Fig. 1 the logarithms of the equivalent widths were plotted against the 
relative gradients. Open circles represent supergiant stars and though their 
scatter, especially for F and G type stars, is rather high, they were included in 
the mean curve. On this figure is also shown the theoretical gradient-equivalent 
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‘TABLE | 


Spectral E.W. Mean 
Type adient in A error 


Ao Ill 03 + 0°06 
\o III "02 “67 11 
Ao V "O04 

\o V Eu ‘o8 . ‘08 
W21V 13 31 


Aa V 
A3 met 
A3 \ 
A3 V 
As Ill 


‘15 
15 
‘il 
13 
13 


As met 
A7 V 
Fo la 
ko Ib 
ko I-Il 


ko V 
Ko V Sr 
F2 V 
Fs5 IV 
6 II 


FS la 
(70 Ib 
Go | 

Go IV 
Gi lv 


G2v 
(s2 Il 
G5 Il 
G5 Il 
Gs5 IV 


G8 IIT 
G8 Il 
Ko 

Ko III 
Ko III 


Kko III-1IV 
Ko III 
Ko III 

Ko 


K2 III 


K2 III 
Ka III 
Ke Ill 
K3 Il 
M3 II 
M3 Il 


log EW 


‘100 
‘212 
"155 
‘240 
‘441 


0°633 
"438 
"403 
"292 
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width relation for radiative models computed by Przybylski (3), with a correction 
of —o-15 applied to the logarithm of the computed equivalent width and a zero 
adjustment of + 1-10 to the observed relative gradients. ‘Che agreement between 
the observed and the theoretical curves is very good and the zero point for the 
relative gradients agrees exactly with the value determined by Greaves (4). 
Our observed curve indicates a maximum value of log EW for gradients approxi- 
mately equal to 2-0 (corresponding to the spectral type Ko) and the values of 
the A line for later-type stars seem to stay fairly constant, within the limit of 
the precision of the measurements. 


log EW my . T T 


1.50 


i - A 
”) 2.00 
AELATIVE GRADIENTS 


bic. 1.-Obser 1 gradient equivaicnl dth relation (Vit fy + determinations) 


A list of contours of the K line for late-type stars was published by Hogg (5). 
‘These contours were integrated to obtain equivalent widths. ‘The 15 stars 
ranging from K2 to M7 (2:1 to more than 4:0 for relative (interpolated) gradients) 
give a mean value for log EW of 1-18 with an average deviation of only 0-06, 
thus showing a fairly constant value of the A line over this range. 

Our observed curve can be compared too with the curve given by Unsold (6) 
and based on Williams’ values (7). ‘The constant intensity of the A line, 
reached for type G1 approximately, is log EW = 1-45, to compare with the value 
log EW = 1-37 on our curve. 

A comprehensive list of A line intensities has been published by Barbier, 
Chalonge and Morguleff (8). ‘Their values for A, determined by the wide-slit 
method, are expressed in logarithms of the total intensity of the line, i.e. 
log ///, for the minimum of the line. A direct comparison of our equivalent 
widths with their logarithms of the intensities has shown that the relation between 
these two quantities is not linear for late-type stars. But if we plot V EW 
against log //J, the points fit closely to a straight line. Fig. 2 shows this relation 
for the stars in common with Barbier, Chalonge and Morguleff. ‘The most 
discordant point is for 6CMaj, a supergiant, which has a low weight in both 
lists. Even in this case the error is only 25 per cent of the value in A. It has 
been shown that such a difference between observers is quite possible; thus we 
felt justified in transforming the AK values of B.C.M. into equivalent widths, to 
allow a direct comparison with our observed curve. 
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Fig. 3 gives for their data the relation between log EW and the relative 
gradients (Stromlo or Greenwich gradients, or the interpolated value corre- 
sponding to the spectral type). The scatter is more pronounced than for our 
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bic. 2.—-Comparison EW Mt Stromlo with K values B.C.M. 


values and the mean curve includes the supergiants as well. When a correction 
of —0°26 1s applied to the logarithm of the computed equivalent width and, of 
course, the same zero point to the observed gradients, the theoretical relation 
computed by Przybylski fits exactly the observed mean curve. 

‘The mean of the two corrections to the logarithm of the computed equivalent 
width, i.e. —o-15 (Stromlo) and —o-26 (B.C.M.), equals —o0-20,, in agreement 





yRADIENT 


Observed gradient—equivalent width relation (Barbier, Chalonge, Morguleff values). 
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with the correction derived by Przybylski from the mean observed curve 
obtained with other data. 

Przybylski’s results do not extend to relative gradients greater than 1-2. 
On the other hand, curves published by Unsdld (6), Russell (g) and Pannekoek 
(10) indicate that the logarithm of the number of Call atoms, or the intensity 
of the K line should decrease for stars of types later than Ko. Of this there 
is no sign in our observations and it seems likely that any correction for difficulties 
in drawing the continuum would tend to increase rather than decrease the 
observed equivalent widths. 

Acknowledgments.—The present work was carried out at the suggestion of 
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THE RADIO SOURCE 
IN ‘THE DIRECTION OF THE GALACTIC CENTRE 


k. G. Smith, P. A. O'Brien and fj. E. Baldwin 
(iteceived 1956 March y)” 


Summary 


Measurements have been made over u range of frequencies from 38 to 
500 Mc/s of the flux density of radio waves from the intense source in 
Sagittarius \erials with high resolving power in the East-West direction 
were used, giving a strip distribution of brightness across the source. It was 
found that the flux density increases with decreasing frequency in this range of 
frequencies, and that the brightness temperatures exceed those expected from 
Hit regions. ‘The suggestion is now put forward that there is in this direction 
a non-thermal source of radio waves lying behind ionized hydrogen, and it is 
shown that this might account for the unusual spectrum indicated by the 
available measurements of flux density, including those at high frequencies 
‘The source itself may be similar to other non-thermal galactic radio sources, 
and its position suggests that it might lie in the galactic nucleus. 


Introduction.— Although the radio source LAU 1752A in Sagittarius, close to 
the direction of the galactic centre, is one of the most intense in the sky, it has 
proved difficult to observe owing to tts peculiar location in the intense band of 
radiation from the galactic plane. ‘his band, and the source itself, are not 
easy to resolve with the aerial beam-widths used, and observations are therefore 
particularly difficult at the longer wave-lengths. In this paper we describe 
measurements using aerials with high resolving power in the East-West direction, 
at frequencies from 38 to 500 Mc/s. ‘the measurements at 38 Mc/s and at 
210 Mc/s were made with interferometers with various spacings up to 60, 
sufficient to resolve the source completely. At the latitude of Cambridge the 
source reaches an elevation of only 10° above the horizon, so that resolving 
power in declination is necessarily very restricted. In right ascension the 
resolving powers of the various aerials were of the order of 1. 

1. Observations.—\'wo ditlerent kinds of observation were made. On 
38 Me/s, 210 Mc/s and 500 Mc/s two small aerials were used as phase-switched 
interferometers of variable spacing in the way described by O’Brien (1953) in 
his measurements of the distribution of radio brightness across the Sun. On 
81-5 Mc/s and 160 Me/s individual aerials of the large Cambridge interferometer 
were used; these were used both in a conventional manner with a receiver 
which measured total power, and also in an arrangement in which two parts 
of one aerial were used with a phase-switching receiver. ‘lhis latter arrangement 


provides better discrimination between radiation from the discrete source and 
that from the intense belt of radiation along the galactic plane. 

The axes of all the interferometers were aligned nearly East-West. ‘The 
receivers for all frequencies were calibrated against thermal sources with an 
accuracy which we believe to be about 10 per cent. 


* Received in original form 1955 August 25. 
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(a) 35 Mc,s.—The interferometer aerials consisted of pairs chosen from 
seven tixed-corner reflectors, each containing one capacity-loaded dipole one 
wave-length long. ‘The aerials were directed at = +7°, and although the 
directivity in the declination plane was very small the sensitivity at 5= — 29 
is uncertain owing to the effects of ground reflection. ‘The flux density may 
therefore be in error by a factor as great as 2. 

Fig. 1 shows the recorded amplitude as a function of interferometer spacing 
(a curve referred to as the visibility curve), and the strip brightness distribution 
resulting from a Fourier transformation of this curve. (The latter curve will 
represent the actual East-West strip distribution provided it is symmetrical 
in R.A.) The width to half brightness is 1°°6, and the total flux density is 
1:0 x 10°"? wm? (c/s)-!. 














Fic. 1. Amp litude-spacing curve at 38 Me/s and the corre sponding trtp brightness distribution 


(b) 81-5 Mc/s.Observations were made at spacings of 12, 15 and 18A, 
using various sections of one aerial of the large interferometer (Ryle and Hewish, 
1955). Each section had a beam of either 4° x 15 or 8 x 15, and in conjunction 
with the resolving power of the interferometer it was then possible to distinguish 
between the discrete source and the general background of radiation near the 
galactic centre. ‘lhe source itself was only partially resolved in the East-West 
direction. Fig. 2 shows the record obtained with a 12A spacing, together with 
the shape of the beam with which the source was scanned in right ascension, 

At the 1955 General Assembly of the I1.A.U. some preliminary results were 
presented of a survey of the galactic centre region made at nearly the same 
frequency with a large cross aerial. ‘This survey, by B. Y. Mills, shows the 
source to be somewhat extended along the plane, so that it appears to have an 
extent of about 3° North-South and 1° kast- West. ‘The East-West distribution 
found in our observations is in good agreement with Mills’ observations.* 

Our measurement of the total flux density at 81-5 Mc/s is 

12x 10°" wm-*(c/s)!. 

(c) 160 Mc/s.—At this frequency the beam of one aerial of the large 
interferometer is 1 x74 . Observations were made of the total power received 
by one aerial, and also by using two halves of one aerial as an interferometer. 
In this arrangement the main lobe of the interference pattern was 1°-2 wide 
between zeros. A record obtained with this system, and the shape of the aerial 
beam, are shown in Fig. 2. On this record the rise at 1} in the following 


* We are indebted to Mills for allowing us to make measurements on a copy of his map prior to 


publication. 
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direction cannot be due to the main source, and must be due to a separate object. 
This separate source is also seen on the record of total power. ‘The asymmetry 
in the 81-5 Mc/s record is possibly to be attributed to this object. 

‘The East West distribution observed with these arrangements was very 
similar to that found at 81-5 Mc/s, and to the distribution of brightness found 
by Mills at 85 Mc/s. The source was nearly symmetrical, and extended with 
fairly uniform brightness over about 14°. ‘The total flux density was 

og x 10°-** wm~* (c/s) 











Records obtain d at SI's Vie s and al 100 VW ) hull line: recorder de flecti nly broken line b 
aerial beam. 


(d) 210 Mc/s. wo movable aerials were used: each one was a broadside 
array of 20 half-wave dipoles in front of a reflecting sheet. ‘lhe receptivity 
pattern was sufficiently narrow in the vertical plane to avoid effects from ground 
reflection. Spacings up to 60A were used, giving the visibility curve in Fig. 3. 
The Fourier transforms here show the source superimposed upon part of the 
background radiation. ‘his background derives in the Fourier transformation 
from the low order components of the visibility curve, which were not included 
in the observations at 38 Mc/s. ‘The width to half brightness in this strip 
distribution is 1°-o, and the apparent flux density is 1-0 x 10°* wm~?(c/s) 
The subsidiary source observed on 160 Mc/s with about one-quarter the flux 
density will have been confused with the main source, giving a flux about 
25 per cent too high, and an uncertainty in the diameter of about o °2. 














a 
-y 
7 


40 


Interferometer Spacing (Wavelengths) 


Amplitude-spacing curve at 210 Mc/s and the corresponding strip brightness distribution. 





No. 3, 1956 The radio source in the direction of the galactic centre 285 


(e) 500 Mc/s.—"l'wo movable broadside arrays were used, each comprising 
f 


22 full-wave dipoles in front of a reflecting sheet. No detailed visibility curve 
could be obtained since the signal-to-noise ratio was insufficient at the larger 
spacings, but it was found to be similar to that obtained at 210 Mc/s. ‘The flux 
density was found to be 0-45 x 10 * wm *(c/s)"!. 

2. Positional observations... Accurate measurements were made of the R.A. 
of the centre of the source at 81-5 Mc/s and 160 Me/s, and these both gave 
R.A. 17" 44™. At 81-5 Mc/s the declination was found to be — 28° + 3°. 

‘These values agree well with optical determinations by Oort (1952) of the 
position of the centre of the Galaxy. Kraus and Ko (1955) have summarized 
the results of other radio determinations, particularly at higher frequencies, 
and these all give similar positions. ‘lhe close coincidence of positions over the 
whole range of frequencies may, however, be accidental, since the source may 
not be one single object. 

3. Discussion of results.-The present observations show the existence of a 
source of radiation virtually unchanged in size over the frequency range 
38-500 Mc/s. ‘lhe measurements with high resolving power at 38 Mc/s and 
at 210 Mc/s show that there is no appreciable fine structure in the source. 

‘The experimental values of flux density are given in ‘lable 1, together with 
the corresponding values of brightness temperature. ‘hese latter values are 
based on our experimental determinations of the East West distribution com- 
bined with the North South distribution given by Mills, and represent an 
average brightness over the source. 

bre querne y 
Vic 
Flux Density - 2 ) 0°45 
10 wim * (c,s) 
Brightness 150 000 ( ’ 300 deg. K 


gs” mperature 


‘The experimental values of flux density are known to an accuracy of about 
20 per cent, except at 38 Mc/s where the value is less certain. 

We note first that our values of flux density are considerably higher than 
those for higher frequencies, which have been discussed by Davies and Williams 
(1955) and by Priester (1955). ‘These authors suggest that the flux density is 
independent of frequency, a relationship which is expected in thermal radiation 
from ionized hydrogen. ‘here is, however, only incomplete evidence on the 
size of the source at higher frequencies; it has been suggested that at gjoo Mc/s 
the source appears to be considerably smaller (Haddock and McCullough, 1955). 

It, therefore, appears that the simple behaviour which we have noted at the 
lower frequencies does not extend over the whole range. We note, however, 
that some experimental values are in doubt: for example McGee, Slee and 


' at 400 Mc/s, lower 


Stanley (1955) find a flux density o-16 « 10 * wm *(c/s) 
than our value at 500 Mc/s, whereas the preliminary results of Denisse, Leroux 
and Steinberg (1955) at goo Mc/s indicate a flux as high as 0°3 x 10 ** wm *(c/s)"!. 
Nevertheless there seems little evidence against our conclusion that, except at 
the highest frequencies, the source can be treated as one simple object, and we 
shall propose a model which will explain the apparent departures from this 
simple picture at high frequencies 
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If we bear in mind that the flux densities at 38 Mc/s and at 81-5 Mc/s 
must be affected by absorption in ionized hydrogen in the galactic plane to an 
extent depending on the distance of the source, it is clear that at the lower 
frequencies the flux density corrected for absorption cannot be independent 
of frequency, but could be represented by a spectrum in which flux density 
decreased with frequency. ‘lhe brightness temperatures may similarly have 
been reduced by absorption, but even without allowance being made for this 
the experimental values greatly exceed 10000 deg. K, the temperature usually 
assigned to Hit regions. ‘hese two conclusions immediately remove the 
possibility that we are here observing thermal radiation from an H11 region. 

4. The identification of the source. —he close coincidence of the position of 
the source with the direction of the centre of the Galaxy has naturally suggested 
that the source is, in fact, a nucleus at the galactic centre. On the other hand a 
remarkable coincidence with the position of a compact group of H11 regions 
has also been pointed out (Davies and Williams, 1955; Priester, 1955), and it 
has been shown that the flux density observed at the higher frequencies is not 
incompatible with such an origin. At 1420 Mc/s, from studies of the absorption 
in the hydrogen line by Davies and Williams and by McClain (1955), there is 
evidence that the source lies at a distance of about 3 kpc, the distance of the 
H 11 region. 

‘The possibility now arises that radiation is received from two sources, one 
of which is the group of H11 regions and the other a non-thermal source which 
may reasonably be considered to lie in the galactic nucleus. ‘The flux density 
of the non-thermal source may be assumed to decrease with increasing frequency. 
Some evidence that the non-thermal source lies behind the ionized hydrogen 

provided by the apparent broadening of the source at 38 Mc/s, which is 
consistent with the effects of absorption by ionized hydrogen lying between us 
and the galactic centre and concentrated within +100pe from the galactic 
plane. 

Under these circumstances the spectrum of the source would show the 
following characteristics : 

(i) At frequencies where the optical depth of the hydrogen is small and 
where the brightness temperature of the non-thermal source exceeds that of the 
hydrogen, the non-thermal source will be seen relatively unaffected by the ionized 
hydrogen. 

(ii) At lower frequencies, say below 100 Mc/s, appreciable absorption will 


occur and the spectrum will become flatter. 


(iii) At higher frequencies, say above 1000 Mc/s, the brightness of the source 


may become less than that of the ionized hydrogen, so that the source may not 


be clearly distinguishable from irregularities in the distribution of the hydrogen. 

s. Conclustons.-We have shown that the radiation observed from this source 
at frequencies from 38 Me/s to 500 Mc/s cannot be thermal radiation from 
ionized hydrogen, although the group of Hit regions which may be mainly 
responsible for the radiation observed at higher frequencies may, however, produce 
appreciable absorption at 38 Mc/s and at 81-5 Mc/s. Allowing for this absorption, 
the brightness temperature at 38 Mc/s would probably reach 400000 deg. K, 
and the flux density vary as (frequency), a spectrum comparable with those 
from other non-thermal sources. ‘lhe absorption varies across the source, so 
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that a unique value of emission measure of the ionized hydrogen regions cannot 
be found without more detailed observations. Values between 10° and 10° are 
consistent with our hypothesis. 

It seems reasonable to suggest that this source is situated in the nucleus of 
the Galaxy. At this distance its size would be about 150 parsecs by 400 parsecs, 
the longer dimension lying in the plane of the Galaxy. Nuclei of dimensions 
as small as 20 pe have been observed in Sc galaxies (Seyfert, 1954) and in these 
a high degree of turbulence is apparent. No nuclei of such small dimensions 


as 150 pe have been reported from optical observations of Sb galaxies, but Kwee, 


Muller and Westerhout (1954) have shown that in our own Galaxy large 
turbulent velocities exist in neutral hydrogen over a region several kiloparsecs 
across. It seems possible that the discrete source represents a central concen- 
tration in this region. 


( ave ndish Lahoratoi 
C'ambridee 
1956 March 8 
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STUDIES IN ‘THE EQUILIBRIUM OF GLOBULAR CLUSTERS (TI) 


R. vw. d. R. Woolley and Denise A. Robertson 
(Received 1956 January 3) 


Summary 


‘Times of relaxation are calculated for radial and circular motion in an 
isothermal gas sphere.’’ ‘The times of relaxation increase outwards with the 
radius. A model is developed in which relaxation is complete only up to a 
limited distance from the centre. Calculations with this model give projected 
densities which agree much better with observation than do those from the 
simple isothermal case 


In a previous paper, M.N. 114, 191, 1954 (Paper 1), attention was given to 
the distribution of density in globular clusters, with special reference to (1) the 
fact that the ‘isothermal gas sphere”’ is not a finite object and (2) that it is only 
a correct solution of the differential equation concerned (Poisson’s equation) 
if the stars all have the same mass. 

The first point has been dealt with by many authors by introducing an 
arbitrary cut-off in the solution. In Paper I this appeared as an arbitrary upper 
limit to the stellar velocities. ‘lhe second point was examined by introducing 
a spread of stellar masses —a mass function. In the present paper we attempt 
to enquire more closely into the nature of the cut-off, in other words to suggest 
a physical justification for one. ‘Throughout the present paper we deal only 
with the case where the stars all have the same mass: firstly, to avoid complication, 
and secondly to avoid producing results which only have meaning in relation to 
an arbitrary mass function. (All stellar masses equal is also an arbitrary mass 
function, but it is the simplest.) Secondly, observations of globular clusters 
continue very obstinately to give no support to the theoretical expectation that 
massive stars are concentrated towards the centres of clusters.* Perhaps this 


only means that the mass luminosity law for nearby stars is quite inapplicable to 


stars in clusters, but it does suggest some caution in assigning a mass function 
and a mass luminosity law to cluster stars, both derived from observations of 
“main sequence” stars. 

1. By “globular cluster’? we mean an assembly of stars possessing spherical 
symmetry. ‘The relation between the gravitational potential ¢ and the radius r 


I d dd 
24 _r 
ah (, 7, ) 471 'p. (1.1) 


If the cluster is in complete equilibrium, Liouville’s theorem can be applied 
to it, and hence so can Jeans’ theorem. ‘The only satisfactoryt solution for the 


° , 
is then Poisson’s equation 


* But the galactic cluster M37 as analysed by von Zeipel shows a concentration of luminous stars 
towards the centre, and so do the Hyades, according to van Bueren 
ft i.e. stable against collisions 
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velocity distribution is Maxwell’s distribution, and if v,,(v, 7) dm dv is the number 
of stars per unit volume at r with velocities between v and v+dv and masses 
between m and m+ dm, 


V»(v, 1) = Nof(m) + (Bm)*?. exp {2Bm(d¢ — dy)} . 0? exp | — Bmv*} (1.2) 
\ 7 


where N, is the number of stars per unit volume at the centre of the cluster, 
{(m) the mass function at the centre of the cluster, and f is a constant. * 

Then if v,,(r)dm is the number of stars at r with masses between m and 
m+ dm, integrated over all velocities, 


Vm (1) = Nof(m) exp \2Bm(¢ : do) § 
leading to an expression for the density, p(r)=Jmv,,(r)dm, which depends on 
the central mass function f(m). However, throughout this paper we confine 
ourselves to the case where all the masses are the same, and then 


v,(r) = No exp {2B8m(¢ — dy) }. 
We write Bm=j? and 2j*%(¢~— 9) y. Then 
P=Poe ' 
Further we introduce a dimensionless radius z defined by 
z =r (87l'p,)*)'* 


and then equation (1.1) becomes 


d (=? dip a2 a ¥ » 
dz\* 5) iil (1-5) 


Equation (1.5) is the well known equation of the isothermal gas sphere and its 
solution has been tabulated. 

In passing we may note that if in (1.4) z=r//, where / 1s a length, then if p, 
is in solar masses per cubic parsee and 7 in (km/sec) |, 

/ = 3°05(py/*) '* parsecs. 

2. We now turn to the calculation of the time of relaxation in the cluster. 
It was remarked in Paper I that “once a high velocity star has got a velocity 
anywhere near the velocity of escape it spends almost all its time well away from 
the centre of the cluster and experiences very few collisions relative to the number 
experienced by a low velocity star which spends its time near the centre.” In 
this section we follow up this remark. 

A formula, number 2.355, given by Chandrasekhar in his Principles of Stella 
Dynamics asserts that UAE? the sum of the squares of the exchanges of energy 
experienced in time dt by a star of mass m, and velocity v, moving through a 
field of N stars per unit volume, all of mass m,, and having a Maxwellian distri- 
bution of velocities with a parameter j, is given by 

VLAE?* = 82 N12m,?2m,20,G(x_) ln (qv,*) dt 
whee G(x) = (2x2) ! 4 erf x - vexp(— v*)> 


‘13 


) . {EP (m, +m.) \! 


Paper I, section 4 
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Here m,;=m,=m. If we take* as our definition of the time of relaxation dt = 7 
when LAF? = E*/e, where E = 4m/4)", then 

‘7° 128 } rT: "4 4 ¥ - . 4 

T >) = re Nj*1*m*v,G( jv) In (qv,”). (2.2) 
lo calculate this for a star moving in a circular orbit is simple, as v, is constant 
throughout the history of the star, but for any other motion it is necessary to 
follow the history of the velocity of the star as it describes its orbit. 

Consider first a star which has the circular velocity at some point in the cluster. 
Since gdr db we have 2j*gldz=dj, where 1=(87I'p,j?)'*. But v2 =gr 
so that 
dif; 
2j*v,? = 2j*glz=2z =. 
JU. = 28 . 
In any motion v*—2=constant = 2E/m where E is the energy in the motion: 
hence /*v* + y= constant = 2j"E/m. Hence if an object moves radially with the 
speed of the circular velocity at a point where y=, the apocentre of the move- 
ment occurs where i) =, such that 


dips 
Wo w,+4(2 o% 
n= +4 (250 ) (2.3) 


and from tabulated solutions of equation (1.5) giving # as a function of z, the value 
of , can always be found for any yy. 
‘The circular velocity at ys, is given by 


. dip\ 2 
J. = (ay)? = 4 i(2 ‘) p (2.4) 
\ le J 
Substitution in equation (2.2) now gives 7), the time of relaxation for a star 
moving with the circular velocity. We have in fact 
Toh = = rep 8 ?m(ih, — b,)'*G (yb, — ,)'*} ln tgp (hb. — y)} (2.5) 


where 


) 13 
q (=) (21m) '{ pyexp(—,)) °°. 


‘Then 
11 my 


en exp (— bs) (dh. — b)!?@G (bh, — hb) * En {qy {by — b)} 


T 
using / * = 8ap,ly*. 

On the other hand if a star is moving radially so that its apocentre, occurs 
where = yy, at any other point /2z2 =y,— 4 so that the time dt required to {pass 
through an element of dimensionless radius dz is 
ldz jlds 

z (yb, —y,)'? > 


Hence the time to complete an oscillation ts 


dt (2.6) 


a : aie " difs 
4b=4jl} (2-H) ( 7 ) ' dip 
y= ™ 
and ‘T can be found from the tabular solution. Now by (2.1) and (2.6) the sum 
of the squares of the energy exchanges experienced in moving through dz is 
LAE? = 8r Ny exp(— ph) F2m'1G (yb, — b)"* | In [gp * (he — Pb) } dz 


*1e. e | of the quantity defined by Chandrasekhar as 7 y. 
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Hence in performing a quarter oscillation (moving from ¢ =o to # =) in time T, 
LAE" ’ 1207 V2m2)4] ure” Gs 2) In fgj-2 , de 
Ge = Nol tmyj'l | “exp(— 4) G (a #'} In (gj * a W)} 
Accordingly LAE*/E?=1/e intime. T7,=n'Tl if 
Ie dz 
I= ne po V¥mjl | exp(— PG —¥)"*} In {gj * he W)} So 
Itence 7,, the time of relaxation for radial oscillation is given by 


Va 


shoTed J exp (— ys) G {(uig — bh)" *} In fq7 2 (bg — oe, (d z/ dibs) ds 
9 . ¥ (yb, i,) '(dz/ dis) di 
0 

Krom formulae (2.2) and (2.8) we can calculate the times of relaxation for the 
extreme cases of circular and radial motion for any total energy E = \mj~*.. 
Both times of relaxation are inversely proportional to ;//*, characteristic of the 
particular cluster, and also inversely proportional to m, the common mass of th« 
individual stars. 

Numerical results are given in Table 1. For the purposes of this table we have 
taken m= the solar mass, j~' = 1 km/sec, and /=1 parsec. 


. (2.8) 


TABLE I 


Times of relaxation for circular and radial motions, etc. isothermal gas sphere 


In z (where z is 


radius) 


| 
dimensionless 2 


log 7 years 
log ‘ years 
log 7’, years 8 , . 
log T/T | g2 | S| 

| y 7 ; 

| 


i 








le 1 P Po 
| 





In the example shown in ‘Table 1 complete equipartition cannot be established 
beyond £ = 3 or 4 ina life of 10" years : and if the life of the cluster is substantially 
less than that the full quota of circular velocities can hardly extend much beyond 
&=3(orz=e%=20). ‘The effect of this defect on the distribution of density in the 
cluster is discussed in the next section. 

3. Weare led to the idea that at some considerable time after the formation of 
a cluster, relaxation is substantially complete at and near the centre, but that in any 
finite time after formation there are some great distances from the centre at which 
there has not yet been time to relax the stellar velocities. If the cluster started out 
in a restricted volume, it would push out some stars beyond this initial volume, at 
first on orbits initially present in the distribution, and later on orbits arising as a 
result of relaxation of velocities at and near the centre. 

We now investigate a model of a cluster in which relaxation is complete inside 
a certain radius R, but stars only occur outside this radius if they travel in orbits 
part of which lie within R (or touch RX). ‘The stars outside R are thrown out, as 
it were, by the equilibrium at and inside R. 

21° 
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Consider an orbit such that the velocity is v, and the direction makes an angle «, 
with the radius vector at R. Then if the corresponding quantities are v, and a at a 
distance x«R,from the centre of the cluster (where x> 1), we have, since the force 
is always central, 

VU, SIN 4 = XVg SIN Ay. 
Again, 

j? (v,? — v4?) = p. — fy. 
Accordingly if V,, V, are the velocities in the orbit which touches both r= R and 
r= xR (i.e. such that sin #, = sin a= 1), 

V,=xV, 
PPV? = (Yq — fy) 07/(x®— 1), J?V 22 = (Yo — Yy)/(** — 1). 
‘Three cases arise, as follows, 
Case |. If j*v,? <b, — 4, the orbit cannot get to xR for any value of a. 
Case I]. If fb, —, <j*v,? <j? V,,’, the orbit intersects r = xR for some values of «, 
(see Fig. 1). 














Fic. 1 —Examples of orbits intersecting r xh. 
(a) % he 
(b) oOo a 
(c) a2,=a,—0 


Case II. If 2y) 2 2V 3 the orbit intersects 7 rR for all Va lues of Ty. 
1 1 1 
‘These cases arise from 


‘ ° ty P P 4 4 , ¢ 
SiN a = SiNa,.—— =sina,.v, {x*v,2—(x?—1) V,? 
2 — 


to 


' 
j 


Hence, sin «> 1, which is impossible, unless «, <A, where 


sin? = x? — (x? — 1) V,?/e,?. 
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When J2v,? = eg — ph, = {(x? — 1)/x?} . 7?V,,?, sin® A =o. 
When v,?= V,3, sin? A=1. 
Accordingly in Case II, o~ sin? A <1, and A can always be found. All (v,, «) 
stars for which o<a<A gettoxR. In Case III, all values of «, give a possible 
value of a, (but not all values of «, are possible.) 

These rules enable us to calculate the number of stars at xR from the number at 
R. If v, is the number of stars per unit volume at R, the number with velocities 
between v, and v, + dv,,moving in directions between «, and «, + dx, is 


vy + j®v,? exp ( —j*v,") dv, | sin x, | dx. 
\ 77 


‘he number moving inwards and outwards over the sphere of radius R in unit time 


Is 
vy 3. pv, exp ( —7*v,") dv, | sin x, | cos «, da, . 47R* 
\ 7 


integrated over all values of «, from o to 7/2 and v, from oto #. ‘The number 
crossing the sphere at xR is the same integrated over all permissible values of o, 
and v,. Each class (v,, «,) contributes to v,, the number per unit volume at xR, 
inversely as 7, COS %, so that 

4 fen4 pars v2 exp(— j*v,*)dv,.sina,cosa,da, R* 


vo = Vy). ee 
2 — - ° ne 
. V TJ a,=9/ jv,2-y —y, Vy COS Ze vk? 


Since v, sin @, = xv. sin a, we have xv. Cos « = (x*v,? — v,* sin, «,)"? and, carryin 
i 1 2 2 2 2 2 1 2% 
out the « integral, we get 


4g [Pm ea , v,*sin?A\ v2 
TrAD sae te. o, UEERP(— S01") a1 4 Mr («" 3 +) is 


"y 0," "Wea 
When pb, — fh, <j?v,? <j?V 4? we have sin *4 = x*V,?/V,* and when j?0,? >j?V,? we 
have sin*A=1. Accordingly 


UV, exp ( —]* ") dv, 


ye 


v,exp(—s*v,*) dv, 


U,% 2 EXP | me 1") de 1 


| o(« s_ ) exp j?v,*) dv, 
ve ~ 


and remembering that /?7, = )?v,* » ~,), vdv, =v, dr, 
Zs )*2 ” 1 
exp { —(y.—y,)} <1 bt y : 
t 1 \ 7 y x 
in which y =/*v," and y, =j?V,? — (yy — hy) = (Wg — yy) /(%? — 1). 
\gain 


9 


“exp ( y)dy, 


-(4 Vy), 


so that 
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The differential equation to be solved is 


5 (# +t) = 2*exp(—y) 


from z=0 to z= R/l (say z= 2z,) subject to =o and a 0 at z=0. (3.2) 


dis a1" | .”( y— py 
z(* z) = 2*exp(— nfo ( Met)" exp J = (i ae | 


4 


for all z>z,, subject to # and s being continuous at z=2,. 4, is the value of 
# at z= 2). 

A solution of equations (3.2) was made, with z, = 30, and carried from z = 30 to 
z=100. It is shown in ‘Table 2, with some comparisons, including projections. 
It was found that if we denote p/p, by 7, then 72° was approximately constant (there 
is a small term linear with z) from z = 60 to z = 200, and in making the projection it 
was supposed that this relation held good outwards for such a distance as could 
sensibly affect the projection. Notice that this gives a projected density roughly 
proportional to r-*, as compared with r ! for the isothermal case. 


‘TABLE 2 


Solutions and projections of isothermal and modified isothermal 








log projected intensity 


log 1) ———_. —— 


mod 80, mod isothermal modification 


0'00 1°52 ae 5 
1°76 1°36 
B°as o'96 1d 
2°35 0°54 ‘75 
3°65 15 2°34 
00 
‘$8 
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To obtain a good comparison of these results with Gascoigne’s observed 
curves,* the zero of the log of projected light was chosen at 
z = 30 for isothermal and modified isothermal solutions 
r= 3.5 minutes for 47 ‘Tuc. 
and r= 7°23 minutes for w Cen 


These four points were plotted as one in Figure 2, thus determining the scale of the 
abcissae. 


* VM.N., 116 (in press). 
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Fic. 2.—Comparison of theoretical models with observations 


The modification certainly agrees better with observation than does the 


unmodified isothermal solution, but much more work is needed, both theoretical 
and observational, to decide these outstanding questions : 
What is the mass function in clusters ? 
W hat is the mass-luminosity law ? 
How far out in radius do clusters extend and what is the nature of the cut-off ? 
Is equipartition of energy confined to the centres, or even absent ? 


Commonwealth Observatory, 
Mount Stromlo, 
Canberra, Australia: 


1955 December. 





ON ‘THE EQUILIBRIUM OF CLUSTERS 
lil. THE PLEIADES 


R. v. d. R. Woolley 


(Received 1956 January 4) 


Summary 


Data concerning the Pleiades, principally trom Hertzsprung’s catalogue, 
are compared with a model (a truncated negative polytrope of index 3) with 
the following constants : 

central density py 15°6 solar masses per cubic parse« 
1 0-60 km/sec, 
truncation 0°72 « velocity of escape at centre, 

‘The projected density, the relative distribution of the brighter stars 
and the observed internal motion are reasonably well represented. ‘The time 
of relaxation of the inner parts is between 107 and 10% vears 


velocity parameter 7) 


‘Lhe Pleiades form a star cluster which has few members in comparison with 
the large numbers seen in great clusters such as w Centauri and 47 ‘lucanae and 
has perhaps been regarded as unsuitable for comparison with equilibrium theory 
on this account. On the other hand the relative closeness of the Pleiades makes 
it comparatively easy to sort out members of the cluster from background stars 
by using proper motions, and also enables us to identify members with low 
absolute luminosities. ‘hese advantages appear to make an attempt to represent 
the Pleiades by a theoretical model worth while. 

Hertzsprung* in ‘lable 4 of his Catalogue de 3259 étoiles dans les Plétades 
lists 2g1 stars which are “ plus ou moins suspectes de participer au mouvement 
propre d’Alcyone.”’ Only one star in eleven in the region is a member: a fact 
which illustrates the great difficulty in eliminating non-members from clusters 
(when stars of low luminosity are taken into account.) ‘The area measured by 
Hertzsprung was 160’ 160’, and according to ‘Trumplert a number of stars 
well outside the area considered by Hertzsprung have the same proper motion 
as Alcyone. Again, Hertzsprung’s limiting magnitude is 16"-8, while 
van Maanen], using the 100-inch telescope at Mt Wilson, found some stars 


fainter than the 17th magnitude which he considered to show the proper motion 


of the cluster. Hertzsprung’s list of 291 stars is therefore almost certainly 
incomplete in extension (distance from Alcyone) and faintness, but it is never- 
theless the most complete catalogue of membership that we possess for any 
cluster. 

‘Lhe stars listed by Hertzsprung are catalogued with photographic magnitudes 
determined by Binnendijk. If the stars are grouped according to apparent 
photographic magnitude, and according to distance from Alcyone, as in ‘lable |, 
it can be seen at once that there is a tendency for the brighter stars to be concentrated 
towards the centre. 


* Ann, Sterrewacht te Leiden, 19 1, A, 1947. 
t L.O.B., 10, 110, 1921 
t Ap. 7., 94, 399, to4! 
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‘TABLE | 
Lstribution of 290 stars less than 85° from Alycone in Hertzsprung’s catalogue, 
arranged in annul of equal area, centred on Aleyone and according to apparent photo- 
eraphic magiitude (One star is 88' from Aleyone.) 





6" to 8" 8" to 0" 10" to12” 2” to a" 14" 016" 
AREA 
O te 20: me: mz: ma: i: 
‘to a2 - a: BB: | i | 
‘tos2 | iT i> mi: i: 
seo}, o = iE: 
61 to67" a | 2 a: 
68 to73 rE B: a. 
74 10 84° | & a> 











We now attempt to estimate the masses of these groups of stars, in terms of 
the solar mass. According to Eggen* the modulus of the Pleiades 1s 5™-68. 
On account of absorption the corresponding parallax is 0":0077 (and not 0°-0073 
as would be deduced from this modulus if there were no absorption). From the 
modulus we can convert apparent photographic magnitudes to absolute photo- 
graphic magnitudes, and then make an estimate of the mass, following Eggen| 
for stars of mass and Kuiper] for stars of mass . Masses are adopted 
as follows: 


Apparent photograph Number of stars Adopted mean mass 
magnitude (Hertzsprung) ( 1) 


We can now form a table showing the distribution of mass in various annuli, 
These annuli are of course as seen projected on the plane of the sky, and the 
annular densities could be transformed into space densities in the manner used 

* A.J. (in press). 


t A.J. (in press). 
t Ap. 7., 88, 429, 193%. 
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by von Zeipel, namely, if p(r) is the space density and o(a) the projected density 
at annular radius a, then 
a'(a)da 


AN= | Geir’ 


’ 


his method was used by von Zeipel himself in his analysis of the open cluster 
M 37, and by van Bueren in his analysis of the Hyades. It has disadvantages if 
the annular function o(a) as observed is so rough that its differential o’(a) is 
unreliable, and also if the observed o(a) is incomplete through lack of knowledge 
of faint stars or of stars with large values of a. 
An alternative method is to start with a model of the cluster and project it 
by using 
‘).rdr 


which has advantages for numerical work as the values of p(r) are theoretical 
and therefore smooth, It is a disadvantage that one has to select a particular 
model of a cluster and make a solution of the differential equation before a 
comparison can be made, and that one is in some danger of finding out more 
about the model than about the cluster. However, in the inner parts of a cluster 
the potential as a function of the radius is not at all sensitive to the choice of model, 
nor is the projected density critically sensitive*. In this paper we have chosen 
a particular model and projected the density (and also the distribution of the 
bright stars) for comparison with the annular results for the Pleiades. 

The model chosen is similar to some described in the writer’s first papert 
but it is more sharply restricted (cut off) than any yet given. — It 1s, in fact, similar 
to models C, D and E in the first paper but the cut off is at ¢ = 2°5 (or k= 4°5). 

‘The notation used is m =m where m is the mass of any star and m a suitably 
defined mean mass. ‘lhe velocity distribution is such that v,(v, r)dudvdr is the 
number of stars per unit volume with masses between pw and y+ dy, velocities 
between v and v + dv and distances between r and r+dr from the centre of the 
cluster, It is supposed that 


v(o,r=C verre” for jet<k—yp 
1@) for je >k— ib 
where k is a constant and ¢/ is related to the gravitational potential ¢ by 
p= 2)°(by— >) 


and ¢, is the value at the centre of the cluster. Here ¢ is reckoned decreasing 
outwards, so that 4 increases outwards from zero at the centre of the cluster. 
‘This arrangement of velocities satisfies Jeans’ relation. 


Integrating over all velocities we have 


v,{r) + e (yy?) *?C,, erty (mk —p)j"* 


where 


ert, x | ye V*dy. 
V7™Jo 





* A. W. Rodgers and D. M. Myers, M.N., 114, 620, 1954. 
ft VLN., 114, 191, 1954 
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At the centre, in particular, 


v,(0) = + (yj2)#2C, erfy (uk)! 
\ 7T 


so that 
v(r)/v,(0) =e ert, y {u(k—)j/erf, y (ue). 
It is necessary to adopt a particular case of C,, as a function of «4, and the models 
have 
C,=const. x p> *e"% 
which gives 
v(0) = const. x we erf, y (pk) 
and 
v(r) =const. x pe? erf, v/ |u(k — Wh)! 


Poisson's equation, assuming spherical symmetry, is 


1d ( , . r 
r dr P dr} er 


and if we put r=/z where /=(87I"p,j?)-'? and p = yp, the equation is 


d (=* dip ns 
dz dz 


subject to =o and p=oatz=o0. From the expression for + v,(r) we get 
( | 3 (R-¥\3? [ (126+ 12k) + 15s") 
I -+ : — “I+ ——_———_—_-—____—- C 
enn #13) k 315 + 84k + 8k? f 


The solution tor k= 4:5, called model F, is given in ‘Table 2. 


‘Taser tl 
Vlodel F 
1) 2 t 1] 
1°000 ; 2°O35 
0'°S39 f 1 ? 0°023 
0°555 i "O54 
O°34!1 H 0°O0O0Qg 
o°210 4 o 2°005 
0'132 i 2 002 
oO O84 I d ool 
oO'OS5 000 
14°54 
Truncation of velocity at centre 
‘Oo 
Ug" 


where i, is the velocity of escape at the centre of the cluster and )*7 


0 


‘The projected densities o(a)/o(o) in model F are as follows 

° I 2 3 } 5 9 2) 10 

1°00 o'39 0°67 o'4 0°33 0°25 160 0°07 0°03 
Here the projected radius a is in the same dimensionless units as are used for z, 
so that the annular radius in real units is a/. Before comparing these projections 
with densities from Hertzsprung’s catalogue we consider briefly the results given 
by ‘lrumpler. For stars brighter than 10™-o photographic, his “ very probable 
members”’ and “‘ probable members”’, distant more than 1 from Alcyone, and 
not given by Hertzsprung, are as follows :— 
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Distance from 6! 9 to 6™-g 7-9 to 7M -g 8/19 to SIG go to of -g 

Alcyone 

rr tors fe) 

16 to 2° 2 
2°10 2°5 I 
2°6 to 4°0 ° 


‘Trumpler’s additions indicate an extension to at least 3° from Alcyone and also 
(for stars brighter than 10"™-o) a twenty-five per cent increase over the mass derived 
from Hertzsprung’s catalogue, which gives a total mass of 337©. We therefore 
suppose that the total mass of the Pleiades is about 500() (probably an 
underestimate) and that a= 14-84, the limit of model F, corresponds to 180’. 
These two identifications fix the scale of a and the mass, and we now proceed to 
compare the projection of model F with the observed values. ‘The quantity 
chosen for comparison is the mass inside annular radius a, and the result is shown 
it Fig. 1. ‘The break in the observed curve at about 60’ from Alcyone is taken 
to indicate that the data are incomplete beyond this radius. 


Minutes of arc from Alcyone 


60 80 100 120 140 160 
T T T T i 











1 l 1 ] 1 
6 7 S.-@ WH 8 
Projected Radius 


Projected ma inside annular radius a from model F (full line) compared wit) 
values derived from Hert prun catalogue and adopti 1 masses (crosses) 


Since in model f 


‘ i 7 
Aé a 2ma da = 4n( 27s), = 331, 


a(O) 


the mass of model F is 331»,/*, in solar masses if py is in solar masses per cubic 
parsec and / is in parsecs. Now since 180’ at a distance of 130 parsecs is 
6-79 parsecs, we have 14°84/=6-79 parsecs or /=0°46 parsecs, and then from 
33 1pol® = 500 we have 


Py = 15°6 solar masses per cubic parsec. 
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Again 1 = (871 'pyj?)"1? so that j~! =1 . (871'py)'*, from which j-! = 0-60 km/sec. 

In a full Maxwellian distribution, the average projected velocity of stars 
relative mass p is }7*(yj*)-"?, but if the distribution is truncated to j*v*< L* 
the average projected velocity is 


V = (uj?) Merf, L/{1 —exp(— L’)}. 


Using k= 4°5 and the empirical mass function derived in the next section for 

the Pleiades the average of the reduction factor 
wo" erfy v/ (u(k—p)}/{1 — exp (— wk —p))} 

is about o-g. Accordingly, the average projected velocity with ; -' = 0-60 km/sec 
is 0=0'54km/sec. ‘This corresponds to a proper motion of 0”-g per 1000 years 
at Eggen’s distance of the Pleiades. Discussing the observed proper motions 
Hertzsprung says “le mouvement interne des Pleiades semble d’étre environ 
+o0"-7/1000a"’. ‘Titus* gives 0”:8/1000 years for the average internal motion, 
and Smith and Struvet concluded from spectroscopic observations that “‘ there 
may exist internal motions of the order of 1 km/sec.” Most weight should be 
given to Hertzsprung’s value of 0”-7/1000 years. ‘The theoretical value of 
0”-g/1000 years is, however, in quite satisfactory agreement with the observations. 

We now compare the adopted mass function of model F with the observed 
distribution of stellar masses in Hertzsprung’s catalogue. Before we can do 
this it is necessary to convert the central mass function of model F 

v,{0) =const. we “ erf, y/ (ky) 

to a mass function for the cluster as a whole, which can only be done by calculating 
the quantities 


A,={v,(o)}-' . 47 | v,(2)22d2 


7 


the mass function of the cluster being /, = const. A,v,(0).  ‘Vhe following values 
were found (constants having been subtracted from the natural logarithms). 


pe 2 I 5 0°50 0°25 


In (v,{0) ) 2°97 1°47 fore) 0°70 1°29 

In (f,) 4°38 2°21 7 2°15 2°75 

log (f,, 4) 1°90 o0'96 00 0'50 o'9 I 
Bw 3 


7 

In An 1°46 o'79 oo 0°42 o's 1°16 
] 
1 


oo I 


From p=0'5 to p=1°5 we have d/dulog(f,/4) = 2-0, nearly. 
On the other hand the mass function of the cluster from m-=o-6o\ 
m= 1-7() is quite well represented by 
d 
—log(f,/m)=1°6 
dm B(Sn/m) 
and by comparison m= 1°25 ( 


The projected annular distributions of stars with = 2 and «= 1 in model F, 
aa) and o,(a), are now given: 


4 


a ° I 
ofa) oO 76 
aa) o'sS oO OS 
a 7 
ofa) 0008 0004 


aa) 0'076 0045 ). 0°004 


* A.J., 47, 25, 1938 
t Ap. 7., 100, 360, 1944 
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With m= 1-250, w=2 gives m=2-5() or about 7™-5 apparent photographic 
magnitude. We can compare the g stars between 7™ and 8" with o,(a) using 
a=12':1, and we find 
Annulus Expected from o,(a) Hertzsprung's catalogu 
‘ ‘ 
o to 12 1°3 I 
13 to 24 2°6 
25 to 36 2°2 
37 to 45 is 
49 to 60 oO" 
61 to 73 ol 
74 O° 
an agreement fully as good as could be expected. 
Stars with w=1 have m=1:25© or about 9-8 apparent photographic 
magnitude. We compare the theoretical distribution o,(a) with 37 stars from 
g™ to 11 and with g1 stars from 8™ to 12™. 


\nnulus 37 stars from 9g" -o to 11-0 gi stars from 3-09 to 12™-o 
Hertzsprung’s 
catalogue 


fertzsprung’s 
Expected o,(a) 
cataloguc 


| 


Expected o,(a) | 
| 
} 








bright stars, well away from Alcyone, do not fit into the scheme: to this 
question we will briefly return later. 

The time of relaxation of the system can now be mentioned. According to 
calculations by Woolley and Robertson*, based it is true on another model 
(the “isothermal gas sphere’), we should expect the following results with 
m=1':25(, l=o-46 parsecs, } '=o-60km/sec and z =e? =7°39 (1.¢. 3°4 parsecs 
from Aleyone, or go’ apparent distance): 

log 7, 7°96 log : 7°47 
Here 7), is the time of relaxation, in years, of stars having the circular velocity 
7°39 and 7’, is the time of relaxation, in years, of stars describing radial 
orbits and having the speed of the circular velocity at a=7-39. ‘These figures 
suggest that if the age of the Pleiades as a cluster is 10‘ or 10° years (or more) 
the stars within 3 parsecs of Aleyone should be substantially in a state of equi- 


ata 


partition of energy. ‘The present analysis appears to support this view. 

On the other hand the bright stars (<8"-o) at distances greater than 1°°5 
from Aleyone seem to be too far out to be in equilibrium with the centre of the 
clustert. ‘There is nothing inherently unreasonable about this, and it may 
indicate that the cluster is not more than 1o* years old. More light would 


VON., 116, 288, 1956 
t There is a suggestion in Fig. 1 that equilibrium may fail at 60° from Alcyone and outwards 
but one would have to be very certain that the data were complete before pursuing this 





| 


he agreement is av good as could reasonably be expected. ‘lrumpler’s 
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thrown on this question, as well as on the appropriateness of model F as a represen- 
tation of the Pleiades, by the collection of more data about membership of the 
Pleiades concerning stars too distant from Alcyone or too faint to be covered 
by the work of Hertzsprung or of ‘Trumpler. 

Acknowledgments.—'Vhe writer is indebted to Olin J. Eggen for many valuabk 
discussions about the Pleiades, and to Denise Robertson for her conscientious 
checking of this paper. 
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Note added, 1956 May 

‘The question has been raised whether the Pleiades can be considered as a 
self-gravitating object, in the presence of the galactic gravitational field. ‘The 
paper examines this question whether the Pleiades behaves like a model cluster, 
by appealing to the observations, and suggests that it does: but the theoretical 
point may be pursued. Clearly a uniform gravitational field would affect the 
Pleiades as a whole: and the question is whether differences in the galactic field 
over the diameter of the cluster are comparable with (particularly greater than) 
the internal field. Stability of a cluster in the presence of the galactic field has 
been investigated by Bok (H.C.O. Circular No. 384). Bok finds that the cluster 
is least stable along that axis which is radial to the centre of this galaxy, and that 
the cluster is stable at a distance r from its centre if the density A(r) satisfies the 
inequality 
ain > 304 
mwla 
where © is the circular velocity in the galactic field, A is Oort’s constant and 
a the distance to the galactic centre. With © = 216km/sec, A = 19°5 (km/sec)/kpe 
and a=8-2kpe we have A(r)> 0-114) pe ® 


If the mass of the Pleiades is 500() and we take the distance r as 6°79 pe. the 


limiting radius in the paper A(r) = 0-353 (pe, which satisties the Bok criterion 
with something to spare. 

The cluster is therefore stable in Bok’s sense. ‘lhe etlect of this galactic 
field is to destroy gravitational spherical symmetry, but not presumably to an 
extent which would easily be observed. It is tempting to suppose that since 
Bok’s limit for a mass of 500() is reached at 10-2 pc, not very far outside the 
adopted limit of the Pleiades, it is the galactic gravitational field actually which 
limits the cluster (causes the cut off)—a point of view which has been put to me 
privately on many occasions by S. C. B. Gascoigne, while discussing cluster: 
generally. 





FINITE ATMOSPHERES WITH ISOTROPIC SCATTERING 
I], INCREASE OF LINE STRENGTH TO ‘THE LIMB 


1. W. Busbridge 


(Received 1955 December 29) 


Summary 


The following problem in the theory of the formation of absorption 
lines is solved by the method developed in Paper I: An atmosphere of 
finite thickness lies above a photosphere of infinite extent, both being stratified 
in plane parallel layers. ‘Vhe absorption lines are formed by coherent isotropic 
scattering in the atmosphere, and there is continuous emission according to 
a linear Planck function throughout both photosphere and atmosphere 
It is required to find the emergent intensity at any point of the line profile 
Subject to certain assumptions, the exact solution is found in terms of the 
X- and Y-functions and their moments. Some numerical results are 
obtained and a theory is put forward to account for the variations in the 
intensities of faint and strong lines across the solar disk 


1. Introduction..-In recent years, most of the theoretical papers on the 
formation of absorption lines have dealt with the Milne-Eddington model, in 
which lines are formed throughout an atmosphere of infinite extent. One 
phenomenon which cannot be accounted for in this way is the strengthening of 
faint lines to the limb of the sun. ‘This was first observed by Dr Adam (1). 
She found that strong lines, on the other hand, decrease in strength towards 
the limb.* 


In some circumstances, a model which comes closer to physical reality is 


that in which the absorption lines are formed in a layer near to the surface and 
such a model will be considered in this paper. For the sake of clarity, the layer 
in which the absorption lines are formed will be called the atmosphere, and the 
semi-infinite layer below this will be called the photosphere, though there will, 
in fact, be no difference between the two outside the line frequencies. 

With this nomenclature, the following problem will be considered: A plane 
parallel atmosphere of thickness x, lies above a plane parallel photosphere ot 
infinite extent. ‘Throughout both atmosphere and photosphere there is continuous 
emission according to a linear Planck function. Absorption lines are formed 
by coherent isotropic scattering in the atmosphere. It is required to find the 
emergent intensity at any point of the line profile considered. 

In (3), to which I shall refer as Paper I, I have developed a method for finite 
atmospheres, which was originally given by Ambartsumian (2) for semi-infinite 
atmospheres. With the aid of this, it is a comparatively simple matter to find the 


* Houtygast (6) found that the wings of some strong lines are strengthened An account of all 
such investigations is given by Unsdéld in (7) 

t ‘his model is almost the same as that considered by Houtgast in (6). In obtaining the emer- 
gent intensity, he used the Eddington approximations 
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exact solution of the stated problem in terms of the X- and Y-functions and 
their moments. ‘This is done in Sections 2~—4. 

Extensive tables of X- and Y-functions (§) have been compiled, with the 
aid of an electronic computer, by Chandrasekhar, Breen and Franklin. For 
reasons which will appear in Section 5, only a few numerical results can be 
calculated with the aid of these tables. ‘Uhese results are, however, of consider- 
able interest, and with these as a basis it is possible to suggest an explanation of 
the phenomenon referred to above. 

2. The Milne integral equation._-We shall solve the problem outlined above 
subject to certain conditions. ‘Throughout both the photosphere and the 
atmosphere we assume that 

(i) the variations of the continuous absorption coefficient « over the range of 
frequencies considered may be neglected ; 

(ii) « is independent of the depth x below the surtace of the atmosphere ; 

(iii) the Planck function B(v, 7) can be represented in the form 


B(v, T) = by + by7, (2.1) 
where 


I 


| Kp dx, (2.2) 


and where the variations of the coefficients 6, and 6, over the range of frequencies 
considered may be neglected. 

Throughout the atmosphere we assume, in addition, that 

(iv) the ratio 7, of the line scattering coefhicient to « is independent of 4 

(v) the scattering in the line frequencies is coherent, 1.e. the frequency 
re-emitted within the line depends only on the radiation absorbed at the same 
frequency ; 

(vi) the coefficient «, introduced to allow for thermal emission within the 
line, is independent of both v and +. 

Let r=7, when x=x,. In the photosphere (7 »7,) and outside the line 
frequencies in the atmosphere (o<-7<7,) the intensity /(7,~) at the optical 
depth 7 and at an angle cos! with the outwards drawn normal satisties the 
transfer equation 


dI(r, 4) 
Yds 


I(r, 2) — (bg + b,7). 


‘The solution of this, subject to the boundary condition 
Ko, ~p)=0 =(O<p<1) 
together with the condition that /(7, 4) is at most O(7) for large 7, is given by 
I(r, +) =by +b(e +7), (2. 
I(r, — 2) = by + by(7 — p) — (by — bye” (2.6) 


where 0: <1. Hence the intensity emerging from the surface of the photo- 
sphere is 

I(7,, t 11) = by + bir, t byp, 
and that emerging in the continuous spectrum is 


Ilo, + L)= by + bp. 
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Now consider the radiation in the atmosphere in the line frequencies. Here 
the intensity /,(7, 1) satisfies the transfer equation 
dl fr, ! . , 
p—a H) =(1+9,)1 (7,4) — 41 —€)n, | Lf(7,p') du 
J] 
- (1 t 7), (bg + br), 


and the boundary conditions 


[(o, —p)~ 0 (O<p<1), 


V7, +4) by +by7, + byp (O<p<t). 
‘The latter comes from equation (2.7). 
If we write 
t, (1 \ Ny)T, t =(1 7 HT 
I + €7), I 


n 
, v , 
I 4 uP | 1, 


and if we replace /,(7,) by /,(t,,) (as can be done without ambiguity), then 
2.9) and (2.10) become 
d1 (t,, ps ;, , , 
é M <1(t,u)—Ma-A)| Ltn’) de 
1 


, dt, 
A,(b, + 6,n,t,), 


1 (0, p= o, 
(ty, + pt) by t byn,t, t bp, 


where Opt. 


Let 
1 


3(t,) = $1 —A,)] Akt eae + A(by + b,n,t,). 
Then (see Paper I, Section 2) 9,(t,) satisfies the Milne equation 
3,(t,)=(1—A,)A, {8,(’)} + BYt,), 
where 
Bt) =A (by + byn,t) + 40 d)| (by + byn,t, + by’ )eO" dy’ 
0 
A (by + byn,t) + 41 —A,)[ (by + bn, t, E(t, — t) + 6, E(t, — t)). 


‘The operator A in (2.16) is defined by [cf. Paper I, equation (2.10)| 


Aft’)} =3{ fe Ele 4h) ae’, 


x 1x 
: ez. 


ete —, 
i 


E,(t)=| 


‘The exact solution of (2.16), when B,(t) is given by (2.17), is a sum of 


multiples of the solutions of (2.16) when B,(t) has the values 1, ¢, F(t, —1), 
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E(t, —t). ‘The emergent intensities corresponding to the first two were found in 
Paper |; those corresponding to the last two are found in the next section. 


3. The solution of (2.16) when B,(t)= E(t, —t).—In this section I shall revert 


to the notation of Paper I in order to be able to quote formulae from that paper 
without change of notation. ‘These will be quoted in the form I(1.1) ete. 
We have to solve the integral equation 
B"(r) =wA{3"(t)} + BE, (7, -7) (Oo<m<1) (3.1) 
for n=2,3,...,08"(r) being the unknown function. We shall solve for the 


functions j,"(4~') and 4,"(u~"), where 


1,"(s)=s O"(tye “dt, (3.2) 


i"(s)=s| B"(r, the" dt. (3.3) 


As in I(2.17), we have the relation 
jut) = 4" — 1). (3.4) 
Since also, by [(4.26), 
Y(p)=e""X(—~p), (3-5) 


the value of J,"(u°') can be written down from that of 1)"(m ') if this is known 


in terms of X- and Y-functions. 
When n= 2, the solution of (3.1) can be found from Section 9 of Paper I in 


which the equation 
(P?(r) =A, {*(t)} + bwE,(7r) (o-m<1) (3.6) 
was solved. In (3.1) the surfaces r=0 and r=7, are interchanged. Hence 
jo (eu )=20 he '), (ur!) = 207" *("*), (3-7) 
and therefore, by I(g.15) and I(g.16),* 
by*(e 1) =2m0 If VoX (2) (1 Xo) ¥ (4) —e heat (3.8) 
(ut) =20 "ti —-(1 ~ Xy)X (1) — Vo ¥ (4) }- (3.9) 
Now let n>2. Then on using the differentiation formula I(3.3) and the 


relation 
d 


- E,, (r)= —E,,_, (7). (3.10) 


We obtain from (3.1) 


v7) 


d O"(r)=wA 
dr 


fd 
< oy" ry } oo” y T 
fy Bp + beB"(O)E (7) 


- ha 3"(7,)E\(7, T)+E,, (71-7), 
and hence 


d n n— = fd n — {r— 
§_ 8'()— Br) =A, } 7, BM) — BU) 


+ 4a3"(0)E,(7,) — wo 3"(7, )E\ (7, —7). 





1 1 
* xy jar [' wrx) du, vr a Y(u) dp, where r ©, 1, By woe » 
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On subtracting from (3.11) 4w3"(o) times 1(4.3) and adding }wS3"(7,) times 
I(4.4), we get the homogeneous Milne equation 


; Zs d x x 
: O"(r) — Br) — ba D"(0)| JS (7, x) = t bar3"(7,) | J (7, -17, no 
ar ! 1 ‘ 


d 


5 O"(t) a” '(t) ha" (0) | Htx) 2 + bar3"(7,)| Ir,-t,x)% s 
1 1 J 


a \ 


r| d 
and, ignoring solutions which increase at exponential rate, the solution of this is 
d ‘ dx = dx 
Fl) -— (7) — ho (0)| I(x, x) = + bod (r1) I(r, —7,*) — =0.(3-12). 
Multiply (3.12) by e*’ and integrate with respect to 7 over (0, 7,) integrating 
the first term by parts. ‘Then, by (3.2) and by 1(4.7) and I(4.10), we have 
CO B"(7,) — B"(0) + Ig"(s) — Ny” (5) 


Z da: “og - 
jo 3"(0)| R(s, x) = ha B" (7, )e omy R(—s, a . 
J1 x J x 
By 1(4.12) and I(4.13), this becomes 
iy"(s) — 8 Ny” (s) = B"(0)S (0, s) — B"(7,)A (74, 5). 
Finally, on putting s=' and using 1(4.21), we have 
Wy"(w*) — pw" (ue *) = B"(0)X (x) — B"(7,) Y(t). 


This holds for n= 3, 4,.... 
The solution of (3.14) in terms of J9?(y~') is 


iM(u?) = X(n){B"(0) + wB" (0) +... + w"™8B4(0)} + p"-44g4(u-) 
Y(u){B"(7,) + pO" (7,) +,-- +p" *B%(7q)}, (3-15) 
and, by (3 4) and (3.5), 
H™(m *) = X(p){B"(7,) — pO" (74) +... + (— pw)” FB(74)5 + (— BHA) 
— Y(p){B"(o) — pB"™ "(0)+...+(—p)™ 2B%(0)}. (3.16) 


The functions j,"(u~") and j,*(u~!) are given by (3.8) and (3.9). 
We now have to determine the (2m — 4) constants 5’(o), 3’(7,) (r = 3, 4,...,m). 
On putting + =0 and 7 =7, in (3.1) and proceeding as in 1(7.12) and I(7.13), weget 


Br(0)= her | sar) du + En), (3.17) 


2. 
3"(7,) =e | HMur) du —. (3.18) 


When the values of J,"(-") and J,"(~) given by (3.15) and (3.16) are substituted 
into these equations we get the following equations for the constants: 


B"(0)(1 — xy) — B""'(0)x, —... — B*(0)x,,_ 3 + 3" (7 yg + B'(7)y, +... 
“| 
+O%(r)ns= or] wom) du + Ear), (3.19) 
B"(7,)(1 =X) + BHT ay. + — 1)" 7 B(7, )X,_-3 + B"(O)yg — BD” “(O)y, +... 


se “ ‘7 - he I 
+ (= 1) 'B40)y_-a=(— 1)" | pW )dut ——. (3.20) 
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From (3.8) and (3.9), 
1 
l 


\ar| phy (un!) du = 2a" | VoX, 2 + (1 — Xo)V, 2} E(t), (3-21) 
0 


' 


4 ; 
bor | p4,2(ur) dum —— = 20 (1 — Xy)X,, 9 + Yona}: (3-22) 
0 n I 


‘The equations (3.19)-(3.22) enable us to determine successively the coefhcients 
B'(0), B(7,) with r=3, 4,.... 
When n = 3, (3.19) and (3.21) give 
B%(0)(1 _ Xo) = (7, Vo =20 1 Yo%) 7 (1 —X)) Ly (3.23) 
Also (3.20) and (3.22) give 
34(0)Vo + B(7,)(1 — Xy) = 27! (1 — Xp )X, + Vey (3-24) 
and hence we have 
3*(0) =2m~'y,, B4(7,)= 2m! x). (3.25) 
On substituting these values into (3.15) and (3.16) (with m= 3) and using (3.8) 
and (3.9), we obtain 
Wye") = 20" Ly, + Yor |X(H) = [x — (1 = Xe] V(H) — ee", (3.26) 
1 °(u ty = 20 M(x, +(1 — Xo) |X (2) aie [yy Yor] (4) — pj. (3-27) 
Equations (3.26) and (3.8) are the ones needed to complete the solution of the 
problem in Section 2. 
4. The emergent intensity. —We now return to the notation of Section 2. 
By I(2.5), the emergent intensity in the problem given by (2.13) and (2.14) 1s 
To, +) =(by + bynt, + byje "+ 5o(e-*), (4.1) 
wher« 
f dt 
tol 4 ‘\= (Le tl , . (4.2) 


0 


This is the sum of 1(7.21) (with a,=A,by, ay=A,hn,), (3.8) multiplied by 
}(1 —A,)(b) + 4n,t,) and (3.26) multiplied by }(1—A,)b,. In each case 7, is 
replaced by ¢t,. Since w=1—A,, we obtain the following expression for the 
emergent intensity : 


1 (0, +p) = X(p){bg(1 — xy) +4,[n,x, + (1 — 1, )y,| + bye[n 1 — xX) + (1 — 1,)¥o] } 


+ ¥(u)fbuvo— Olt =) +m) + Old —M,\r—a%y) +m Vo]}- (4-3) 


In (5), Chandrasekhar and Elbert have given the moments «,, B, (r=0, 1), 
where 


1 1 
4,=| X(p)u" du, B, Y(j2)e" dp, (4.4) 


so that, in our case, 


x,=h(1—A,)x,, —¥,= 3(1—A,)B,. (4-5) 
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Hence for calculations (4.3) must be written in the form 
Lo, +p) = (1 —A,)X(u){bo[2(1 —A,)* — 49] + by [may + (1 — ,) By] 
+ byp[n(2[t —A, J! ~ a9) + (1 ~m,)Bo)} + 401 —A,) ¥(#) 
x {boBy — b,[(1 — n,)a, +n,28,] + bp[(1 —n,)(2[1 —A,}-! — a) +n, Bol}. (4-6) 


As a check on this formula we can let t,-> ©. Then X(y)-> H(y) and 
Y(u)-—> 0, and we have 
10, +p) - s( —A,)H(p){bo[2(1 -A,y* — ty] + bn, 
+ b,pn,[2(1 —A,)-! — ay} }, 
where 
! 
z,=| H(p)y" dp. 
0 
By (3), Chap. V, Theorem 1, %, satisfies the equation 
b(1 —A,)ay = 1 —A,"", 
On using this and the values for A 
equation (66) of (3), Chap. XII. 
5. Numerical results.—'Vhe residual intensity in the line at frequency v 
(expressed as a percentage) is 


and n, given by (2.12), (4.7) reduces to 


v 


r,(u)= 1001 (0, + )/I(o, +p), (5.1) 


where /,(o, +) is given by (4.6) and J[(o, +) by (2.8). In the following 
calculations we have taken b,/b) = 1-5, «=o and (for comparison) € = 1. 

If «=o, then n,=A,. On writing w,=1-—A,, we get the following formula 
for r,(p): 
rp) or,X (42) 


\2m Le + S(Alay +a, B,) + iplA (2m, 


* 34 


Ay) + w,Po| } 


‘plo (20, . ty) + A,Bol}- 


The X- and Y-functions are tabulated in (§) as functions of 7, w aud px, where 
T=0°05, O10, O15, 0°20, 0°25, 0°50, 1°00, 
m7 = 0°50, O80, 0-gO, 0°95, 1°00, 

and o<p<1. In our case, 


r=t,=(1+7,)1 o=m,=7,/(1 +7,)- 


v 


Thus we get the following possible values of y,: 


y= 1s 45 9 195 


and corresponding to these, the largest possible values of 7, are 


respectively. 
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The following table gives the values of r,(1) and 7,(0) whenever these can be 
calculated from the tables. ‘lhe values in italics were estimated graphically. 
They could, of course, have been obtained by interpolation, but | doubt whether 
greater accuracy would have been achieved. ‘The last line of the table gives the 
Milne—Eddington case (rt, = *). 


TABLE | 


r(1) and ro) 3 


96°93 





g0°00 


88°50 





80°43 


Re 26 





‘The value of r,(o) has been given in ‘lable | because it enables us to see at 
a glance which lines will be strengthened towards the limb for the values of », 
shown. Moreover, by considering the difference 7,(1) —7,(0), it is seen that the 
change in residual intensity is almost independent of », for a givenr,. Hencein 
studying the centre to limb variations of r,() it is enough to consider the values 
of », and 7, corresponding to the extreme right-hand member in each row. ‘The 
corresponding values of r,(j) for 0-2 <j <1 are given in ‘Table I]. When7t,= « 
I have taken »,=1 because this gives about the same intensity when »=1 as in 
three of the other cases. 


"? 


"TABLE II 


Centre to limb variations in ¥ <4 
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From ‘lable If the following conclusions can be drawn: 

(i) ‘he shallower the atmosphere, the more are lines (for the given values 
of »,) strengthened towards the limb. For 7,=0-°05, 4,=19, the greatest 
strengthening, viz. max [r,(1)—r,()], is about 7 per cent of the intensity at the 
centre. For 7,=0°5, 7,=1, it is only about 1-4 per cent. 

(ii) ‘The greatest strengthening occurs at about pp =0-4 for r,=0-05 and the 
position moves towards the centre with increasing 7, until, for sufficiently large 
7,, the line is strongest in the centre. 

When «=1 then [by (2.12)] A,=1 and w,=0. As wm, +o, X(w)->1 and 
Y(u) -e &".* Hence on letting A, 1 in (4.6), we have 


[(0, +p) = by + bn, + b(t ~n, pe ’” (5.6) 
and, with b,/b, = 3/2, 
r (pu) = 100{1 — 3u(1 —n, (1 —e 8") /(2 + 3n)}, (5-7) 


a formula which can be verified independently. 
‘The following table for «= 1 corresponds to ‘lable II. 


"Tasie Il 


Centre to limb variations in r,(p) when €= 1 


69°66 $104 
69°17 80°73 
68°87 80°54 
OS S4 80°53 


69°26 sO 





79 
70°36 81°47 
72°46 82°73 
70 ob SS O4 y. 48 


K-66 SRS 4 ®! “46 

















Krom ‘lable III it is seen that the greatest strengthening is in every case less 
than 1°6 per cent of the intensity at the centre; the position of greatest strengthen- 
ing is near the centre (u=0-7); the lines fade away rapidly towards the limb. 
(Actually r,(j) + 100 as pO.) 

6. Discussion and conclusions.—'Vhe values of r,(j) found above for 7, <0°5 
are comparable with the residual intensities in faint lines or in the wings of strong 
ones. Unfortunately, the results given by Dr Adam in (1) are mainly expressed 
in terms of equivalent widths and are not comparable with the results found here. 
It is evident from ‘Tables I] and III that the observed strengthening can be 
accounted for theoretically only by taking « =o or, perhaps, small, and the follow- 
ing discussion will therefore be concerned with the results when «=o (‘Tables I 
and II). 

First consider the values of 7,(1) in ‘Table I. For a given value of », the 
difference between the yalues of r,(1) for two different values of 7, gives a measure+t 

* This follows at once on letting a-+o in equations (4.22) and (4.23) of Paper I. 


+ ‘his is not strictly correct because the values are calculated on the assumption that there is no 
scattering in the layers below them 
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of the contribution of the corresponding layer to the total absorption in the 
Milne-Eddington model. ‘Thus for »,=4, the total absorption 1—r,(1) ts 
57°19 per cent, of which about 12 per cent arises from the layer 0<7<0°05, 
8 per cent from the layer 0-05~ + 0-1, and so on. ‘The larger the value of »,, 
the larger is the contribution from each successive layer, the increase getting 
bigger as the layer approaches the surface. In particular, when »,=19 a 
contribution of about 38 per cent towards the total absorption of 77°22 per cent 
comes from the topmost layer. ‘his is what is meant by saying that strong 
lines are formed high in the atmosphere and weak ones low down. 

Nevertheless, it is also true that the larger the value of r,, the stronger will be the 
total absorption over (0, 7,) for a given »,. Moreover, the same residual intensity 
can be produced by taking different pairs of values of 7, and y,. For instance, 
r,(1) is approximately 66 per cent for 4, =9, 7, =0°1, for y,=4, 7, = 0°2, and for 
MN =1, 7, = © (cf. Table Hl, line 1). Only a more detailed knowledge of the 
conditions in the atmosphere of the Sun (including the variations of 7, with depth) 
than we at present possess would enable us to decide between these pairs of values 
in the case of a particular line in the solar spectrum. 


However, ‘lable II shows that, even though the values of 7,(1) may be about 


the same, the behaviour of r,(u) from centre to limb is very different. ‘The 
smaller the value of 7,, the more is the line strengthened towards the limb. ‘This 
suggests that faint lines which are strengthened towards the limb must be formed 
in layers of optical depth 7, less than o-5 (for which the strengthening is negli- 
gible); lines which are not strengthened must be formed in layers of optical 
depth 7, greater than 0-5. Since, as has been emphasised above, the total absorp- 
tion over (0,7,) for a given 7, increases with 7,, the former would probably be 
weaker lines than the latter. 

In order to account physically for the strengthening of any particular faint 
line according to this theory, one would have to show that, for the line in question, 
n, varied with depth in such a way that it could be well represented by a constant 
value for (say) 0° 70-2 and by zero for r 0-2. 

Acknowledgment..-1 am indebted to Dr Adam, Dr Stibbs and Professor 
Plaskett for criticisms from the physical point of view and tor some helpful 
suggestions. 


St. Hugh's Colles 
Oxford: 
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SOLAR ATMOSPHERIC HENTING BY HYDROMAGNETIC WAVES 
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Summary 


Previous estimates of solar atmosphere heating by hydromagnetic waves 


(due to Joule loss and viscosity) showed negligible heating \ new mecha- 


nism is described, involving ne al atoms, which leads to absorption (and 

heating) many orders of magnitude greater than previously believed. If 
hydromaygnetic waves are present they must cause strong heating effects 

‘The formation of such waves from granules is described, their absorption 

it low levels, and partial reflection being discussed in some detail A very 

rough order-of-maynitude est mate of cnert flow »btained The region in 

hich heating effects are important is restricted to the transition zone between 


lightly sonized and fully tonized gas The energy absorbed in this region is 


estimated as ™10° erg cm™*. sec™! Numerous particles with energies 
~™100 EV are released in this region and these may account for the quiet (non- 
pot) corona 


The mechanism is likely to be important in connection with sunspot 


heating and flares and perhaps also t el! hydromagnetic waves and 


cosmic ray veneratior 


Introduction. Of the ny theories advanced to explain the heating of 
the solar corona and chromosphere the most promising are those depending on 
energy transfer by wave motion upwards from sub-photospheric sources. ‘The 
types of wave which have been considered are hydromagnetx (Alfvén (1, 2)) 
and pressure (Biermann (3) and others) 

Recent observations of solar magnetic fields (4) indicate that these are 
widespread, covering polar and other extensive areas as well as spot regions. 
Furthermore the greatest heating ctlects are generally observed in regions 
permeated by magnetic fields Hence if wave are effective in transporting 
energy they travel through tomized gas permeated by a magnetic field. Four 
and only four types of wave may propagate in such a medium (5) and at the 
frequencies concerned (presumably much below radio-frequencies) one of these, 
the space-charge or electron pressure wave, is excluded. ‘The remaining three 
are all hydromagnetic waves although one has pressure (sound) waves as a special 
limiting case. ‘Thus, in general, any waves which heat the solar atmosphere 
must be hydromagnetic waves. 

It is desirable here to separate hydromagnetic waves into two categories (see, 
tor example, ref. (6)). “lwo of the waves have velocities of order V = Ho(47p)"' ?# 
where H/, is the magnitude of the steady magnetic field and p the gas density. 
‘These waves are largely transverse, they are called the ordinary (O) and extra- 
ordinary (I) waves according to their polarizations and have as limiting cases, 


when heavy ion motion is negligible, the corresponding radio O and E waves of 


the well-known magneto-ionic theory. ‘The remaining (S) wave has a velocity 
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of the order of sound and is mainly longitudinal ; it may be regarded as a moditied 
sound wave. ‘The present paper is primarily concerned with the O and E waves. 

Alfvén’s theory concerns the O and E waves; he considered that. the 
associated electric fields and currents might cause Joule loss sufficient to account 
for the heating. Cowling (7) investigated the Joule loss and decided it was much 
* A more detailed investigation (8, 6) 
gave an estimated chromospheric heating rate about 10° or more times lower 


too small to cause significant heating. 


than that required to explain the observed effects. ‘lhe possibility of ettective 
heating by the O and E waves appeared to be remote. Incidentally, similar 
considerations (and others) dispose of the possibility of effective heating by 
more or less steady electric currents, a mechanism which has been suggested by 
a number of authors. 

There is no theory of heating by the 5 wave but such would be a logical 
development from existing theories invoking sound waves. 

It is the object of the present paper to re-discuss the absorption rates of the 
© and E waves in the solar atmosphere. A new factor, involving neutral atoms, 
is revealed which may enormously increase absorption; in some regions factors 
like 10!” are involved. 

\lthough most solar atmospheric heating appears to take place near sunspots 
the new heating mechanism is not discussed here in that connection. ‘he 
example chosen for more or less detailed analysis concerns O and E waves caused 
by solar granules in regions of weak but extensive magnetic fields, such as neat 
the poles. ‘This simpler example allows a discussion of the formation and partial 
reflection of the waves, as well as their absorption. However, the results given 
and some further calculations indicate that the mechanism may be very important 
in connection with heating near spots, flare formation and interstellar phenomena. 

2. Heating by hydromagnetic waves. he Joule loss of hydromagnetic waves 
has been found (6) from conductivity transport equationst and in the solar 
itmosphere is generally negligible. Another effect which causes absorption 
and hence heating, is viscosity which has been considered by van de Hulst (9). 
Lhe ratio of viscosity loss to Joule loss is 40 where a is the electrical conductivity 

nd the kinematic viscosity (c.g.s. electromagnetic units used throughout). 
‘The ratio may be considerably greater than unity in parts of the chromosphere 
but the heating effect due to viscosity is still negligible. 


When neutral atoms are present in numbers comparable with the heavy ton: 


they have negligible effect on the electron collision frequency the reason 1 


their very much smaller collision cross-section. However, for varying electric field 
hydromagnetic waves), neutral atoms do change the electrodynami 

rtics of the gas since the mass of the gas can no longer be considered situated 

ions (as In a binary gas) but to some extent in the neutral atoms. ‘he two 

yases tend to move independently, thus causing a frictional loss. Consider a 
rather extreme case when the collision period of a neutral atom is of the same 
order as the wave pertod. Since the electromagnetic forces necessary to maintain 
hydromagnetic oscillations act only on the ions, the ion plasma (electrons and 


lo 


* ver is red mate of Joule le 
ave (Cowling’s Lo i ould be identified with the ivelength in the corona (6) and not with 
! 


the size of the vranule If the tead magnetic field 1 100 wal ! t na further reduction 


appears much too large hie dimension " of the 


of absorption by a factor of 10° or more 


t As opposed to Maxwell's momentum transfer equations wh ire also port equation 


! 


combination of the two ts found necessary below 
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heavy ions) tends to perform oscillations independently of the neutral atom gas 
while collisions between the gases tend to promote forced oscillations of the 
latter as discussed previously (ref. (10), p. 660). At times there will be con- 
siderable relative velocity between the two gases so that a collision between a 
heavy ion and a neutral atom will cause a substantial dissipation of wave energy, 
which reappears as thermal, excitation or ionization energy. On the other hand 
all ions share the common hydromagnetic oscillations so that collisions between 
them are not significant in this way. 

‘The new dissipative process is not an ordinary Joule loss; nor is it an ordinary 
viscosity effect since it occurs (for finite wave amplitude) when the gas particles 
have no thermal motions. It may, if desired, be regarded as Joule loss due to 
Hall drift in the ternary gas. ‘This drift may be large but as the plasma moves 
nearly as a whole the Hall current may be small or zero. 

When the atom collision period is much less than the wave period the two 
vases move nearly in unison. However, the loss may be far greater than the 
ordinary Joule loss even under these conditions. 

3. Hydromagnetic waves in a partially ionized gas.--Vhe properties of hydro- 
magnetic waves in a partially ionized gas may be studied in either of two ways. 
They may be found from a set of field equations already given (ref. (8), 
equations (14), (16) and their companion), combined with a dynamic equation 
describing the gas velocity in terms of the various forces acting on it. ‘The 
conductivity components are those of a partially ionized gas; they comprise a, 
parallel to the magnetic field and o, and o, perpendicular to the magnetic field, 
in and perpendicular to the plane of the electric field respectively. ‘lhe 
composite or “ effective’’ conductivity o, is given by 0, o, + 4,2/o, and when 
small is associated with strong damping of the wave.* ‘The gas velocity is the 
weighted mean velocity of the ion plasma and neutral atoms. 

‘The alternative approach, permissible under conditions described below, is 
to regard the ion plasma (heavy ions and electrons only) as an entirely separat« 
yas. ‘The relevant field equation is that for a fully ionized gas (ref. (8), equation 
(17)): 


H as 
V*H = 470, < curl(v » A) ? ~ curlH, (1) 


Oo, 


where H is the magnetic field, having a steady component H, along Oz. ‘Th 
conductivity Components and gas velocity v now refer to a fully ionized gas; 
this means that o, and o, are large when the collision frequency is low. 

‘The effect of the neutral atoms is introduced through their frictional drag 
on the ion plasma and its effect on v. ‘Vhus, instead of a single gas momentum 
equation, there are two. For simplicity the terms in gas pressure are omitted 
from these equations. ‘lhe relative magnitude of these terms in the region of 
main interest (mid-chromosphere and above) is small and they may be neglected. 
At lower levels they are large but it has been shown (6) that even under these 
circumstances the velocity and absorption of the O and E waves in a fully ionized 


* The small value of o, in a partially ionized gas when collision frequencies are low has been dis- 
cussed (10) for the case when only collisions with neutral atoms are important. Since the com- 
munication of the present paper T'.G.Cowling has kindly sent me a copy of his paper The dissipation 
of magnetic energy in an ionized gas (M.N., 116, 114, 1956) in which he completes the discussion and 
shows that the effect may lead to unexpectedly rapid decay of magnetic energy in certain cases ot 
istronomical interest. 
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gas are unaffected because the waves are mainly transverse. It seems safe at this 
stage to neglect the pressure terms, in which case the momentum equations are 
ov nV 
Ei lated i H,? 


n 


(H,, x curlH) =o, (2) 


ov’ = 

at bas ie. : 
where v’ and 7, are the velocity of the neutral atom gas and the collision period* 
of an atom and 7 is the ratio of mass densities of neutral atoms and ions. 

Strictly speaking, either method of investigation is applicable only to weak 
waves (10), although the results may be suggestive of the effects of strong waves. 
No further limitation in this respect is necessitated by the use of separate gas 
velocities since (v—v’) and also its effect on collision periods must both be small. 
‘The second method of investigation imposes a further restriction: the collision 
rate of electrons with heavy ions must be substantially greater than with neutral 
atoms; this is assumed in calculating o, and o,. Because of the larger collision 
cross-section of heavy ions the requirement is satisfied even though neutral 
atoms are about 1000 times more plentiful. 

An advantage of the second method is that it discloses more about the relative 
velocity of the two gases; this is particularly important when the wave period 
is comparable with the neutral atom collision period. ‘lhe second method may 
also reveal more about the underlying physical processes in the wave and is 
adopted here. 

Plane wave solutions of equations (1)-(3) are sought of the form expi (wt — ké) 
in a right-handed system x, y, z where O€ lies in the xz plane at an angle / to Oz, 
the direction of the steady magnetic field. Now let H,, H,, H, be the maximum 
values of the components of the perturbation field (H — Hy) and v,, v,, v, those 
of the plasma velocity v. In all cases the perturbation field lies in the phase 
planes so that H, tanwH,. Substitute in the component equations : 


« 


lw, = ik sin ys, 
x 
ik cos Wh, 


O - :I + ne t lwT ,) 1} - 


and 
V =Hy(4np) *, 


p being the mass density of the ion plasma alone. We find 


(hk? + iew*)H, + ie, cos kw , 4 2 cos* pk? H_, = 0, (4) 


(k? 4 iew*)H, t eH, coskwv, io k*H =o, (5) 
a; 


* ‘The mean period which elapses before atoms loose their ordered motion relative to the wn 
plasma. In equations of this nature, the mean velocity of the whole gas mixture is usually 
introduced and the collision period defined as the mean period which elapses before the atom 
loose their ordered motion relative to the gas’ mixture (that is their peculiar velocity). ‘The 
simplification results in a different value of collision period, the factor involved being small 
enough to neglect under the conditions assumed here. 
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H,,cosO leyz a” V*kH , = Oo, (6) 
H,O we, + V* cosphH, =0, (7) 
which are identical with those of a binary gas (6) except for the term O. 

4. Solution of the equations..-Vhe equations are solved as before (6); for 
convenience terms (eo,V*)/o, are replaced by Q/w where @ is the heavy-ion 
gyromagnetic frequency.* ‘The polarization of the waves (R=(H,,cos#)/H,) 
is found to be: 


R 


1 {OQ sin? y ; ts: .; 3 (8) 


2\ w cosy w* cos* 


Since O 1s a complex quantity, R has both real and imaginary parts; the waves 
are partly plane, partly elliptically polarized. In the absence of neutral atoms 
they are elliptically polarized. 
The dispersion equation is found to be: 
ro (v4 ! 
w*/k? Ol 282 + - (g) 
€ 


where 


S a (10) 


sin? i on w"* cos" is 12 
‘\ 4 020? 


> 


! and as we are dealing 


with much lower wave frequencies the last term of equation (10) may be 
neglected so that S=cosy for the O wave; S=1 for the E wave, just as for 


In a magnetic field of only one gauss £2~10' rad sec 


waves in the absence of neutral atoms. ‘lhe corresponding wave velocities are 
O'?V cosy and O'*V. ‘The quantity O has real and imaginary parts, the latter 
being caused by the presence of neutral atoms. When 7-0, O-»1 in agreement 
with the results for a binary gas (6). When wr, »1 and wr, > so that the two 
gases move independently we again have O-»+1 as would be anticipated. In 
this case, although the wave velocity is not appreciably affected by the neutral 
atoms, there may still be considerable absorption. When w7, <1 so that the 
two gases are effectively bound together, O~(1 +7) ' and the velocities are 
I’ cosy and V'=V(i +) '*, being smallert than V because of the additional 
yas density. Finally when w7,~1 and » is not small a type of forced oscillation 
occurs associated with a considerable change in wave velocity and a high degree 
of absorption, 

‘The wave absorption coefficient « is equal in magnitude and opposite in sign 
to the imaginary part of k. When no neutral atoms are present (O=1) this 
absorption is due to collisions between electrons and heavy ions and is given 
by the last term in equation (g). ‘This contribution has been shown negligible (8) 
and will be neglected for the time being. ‘Thus we have: 





* This is seen by substituting for e and V, which were defined in the previous section, and also for 
o,, a, and oy whose values for a binary gas have been given (see, for example ref. (10), p. 658, noting 
that there the term w represents the electron gyro-frequency H,e/m, where m is the electron mass). 
‘The ton gyro-frequency ts £2 //,e/M where M is the ion mass 

+ V’ and not V was the velocity previously found (6); the apparent discrepancy occurs because 
in the present paper the density determining V is that of the ion plasma alone. The real hydro- 
magnetic velocity, provided the gases move together, is V’’ and this value is used below. 
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When w7,, <1 (wave frequency is much less than the relevant collision frequency) 
this reduces to: 


2 2 
WT »7] wT 7) 


«= 275( + m8 = 2S +9): me 


This absorption coefficient, due to friction between the two interpenetrating 
gases, may be compared with that due to Joule losses (ref. (6), equation (20)). 
At the chromospheric level where the hydrogen is half ionized (7 = 1) the electron 
density is about 3 x 10°cm™~*, the conductivity o, about 2 x 10-%e.m.u. and tr, ~1 
(see Appendix 1). Ina steady magnetic field of 100 gauss the absorption given 
by equation (11) is about 10'° times greater than that due to Joule heating (and 
10” times greater than that due to normal viscous effects). In the corona the 
ratio is lower but heating by either mechanism at these levels is unimportant 

5. Hydromagnetic waves from granules.-Granules are convective cells moving 
upwards with velocities of about 1 kmsec ' and having lifetimes of a few minutes. 
The gas density in the granule is about 107‘ gcm * so that the upward flow of kinetic 
energy is about 5 x 10’ergcm™*sec’', perhaps a thousand times or more that 
needed to maintain the corona in non-spot regions. Alfvén (2) has suggested 
that some of this energy is transferred upwards by hydromagnetic waves. 

‘The actual processes involved in the upward transfer of this energy are largely 
speculative. However, there are sufficient data available to allow a semi- 
quantitative discussion which may be of some value. 

Upward transfer of energy by mass motion may be ruled out. Even if the 
granules were moving into a vacuum above the photosphere, its total energy 
(thermal motion of its particles plus mass motion) would only raise it one o1 
two hundred kilometres. If it were being thrust up from below by a larg: 
pressure gradient then the whole granule would be seen to emerge into the 
chromosphere.* ‘The transfer must be by wave motion and in a weak magnetic 
field of say 5 gauss the only effective wave is the S wave (see Section 1). ‘Vhis 


has a velocity of about 1okmsec! 


compared with about o-2kmsec ! for the O 
and E waves, which gives a measure of the relative effectiveness of energy 
transfer. 

As the waves travel upwards the density decreases while the magnetic field 
strength remains nearly unchanged. ‘The elastic forces in the medium become 
increasingly electromagnetic in character and there is a tendency for the O and 
k waves to grow at the expense of the 5S wave. ‘This process is no doubt helped 
by collision between waves from different granules and collision of the waves 
with irregularities in the gas distribution. 

A critical level is reached at about 1500 km above the photosphere where 
(when H,=5 gauss) the gas pressuret and the magnetic pressure //,?/87 attain 
equality. Above this level the O and E waves outstrip the S wave attaining 
velocities of ~130km sec! at 6000 km where the gas density has a value 
~2x10"%gem'!. It is likely that a significant proportion of the energy of 
the original granules is contained in these waves. 


* Additional energy in the form of higher gas temperature or mass velocity, as suggested by 
‘Thomas (11), also seems to be contrary to observational evidence 

+ For this and other such data a model of the lower chromosphere given by van de Hulst (ref. (12) 
p. 238) is used. 





320 J. H. Piddington Vol. 116 


‘The ‘ period” of the granules is a few (say 3) minutes (12), corresponding 
to an angular frequency w ~0-03 sec '. ‘The O and E waves will contain 
harmonics of this frequency, the proportion of which will be increased by 
collisions with irregularities and with other waves. An application of the 
absorption formula indicates that waves with angular frequency 0-3 sec™! are 
absorbed too strongly to be effective at the upper levels. We will assume that 
ten per cent of the original energy (that is 5 x 10°ergcm™~*sec™') is retained in 
the form of O and E waves of frequency 0-1 <w<o-3 and determine the likely 
effect of these waves. 

‘The fate of the 5S waves is not clear. When propagated along or in the absence 
of a magnetic field they are sound waves and the evidence suggests that they would 
be dissipated around the 1500 km level (12). However, when propagated 
obliquely to a magnetic field, they should have greater coherence and may reach 
higher levels and possibly play a part in heating the corona and perhaps in the 
formation of spicules. 

6. Absorption and reflection..-Up to the critical level near 1500 km the 
wave motions are controlled mainly by hydrodynamic forces. Hence the 
absorption discussed above, due to electromagnetic forces, is not appreciable. 
At a level of 2000 km the gas density ~5 x 10° gem * and if H,=5 gauss, 
' ‘The value of 7, will be about one-tenth that at 6000 km (given 


’ 


in Appendix 1), that is r,~o-1 sec. Equation (11) then gives «~10~* so that 


V’ ~20 km sec 


the wave decays by a factor ¢ in a distance 1o'km. As far as absorption is 
concerned most of the wave energy should reach a level of 6000 km. 
Reflection of plane hydromagnetic waves in a horizontal, stratified atmosphere 


with wave normal and //, both vertical has been discussed by Ferraro (13). 


‘The gas density varies as p = pye “** where 2, is the scale-height of the atmosphere ; 
the corresponding variation in the perturbation magnetic field is found to be 
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where J, 1s Bessel’s tunction of order unity and V,’ the wave velocity at the base 
of the atmosphere. Since e 7) =(p/py))'*=V,'/V’ the Bessel function may be 
written J,[(472)/A] where A is the wave-length at level =. When (272,)/A>1 
the perturbation magnetic field is found to vary as p'*. But this is just the 
variation necessary to maintain a constant Poynting flux in the variable-density 
medium so that no reflection occurs but merely an increase in wave-length. 
When the latter increases to equality with the scale height, reflection occurs, 
the last maximum of the standing wave pattern occurring when (272%,)/A~o-g. 
Above that level the field decays by a factor of 10 in about three scale-heights. 
Ferraro’s discussion concerned a medium with z, constant. However, it 
seems safe to assume the same general result if 2) varies (as in the chromosphere) : 
that is, no reflection until (272,)/A=1. ‘The results may be further extended 
to the case of oblique propagation provided we consider only the O wave whose 
perturbation magnetic field remains perpendicular to the steady magnetic field 
(6). ‘The equations remain unchanged except that 2, is replaced by 2, sec y, 
where ys is the angle of inclination of the wave normal to the vertical. It is seen 
that the waves penetrate not only to greater distances but to greater heights. 
In the middle chromosphere 2,=1200km or more so that reflection 
commences above a level where A~7500okm. If w=o-2radsec', then 
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V’=240kmsec™!, corresponding to a density of about 3-5 x 10 “gem and 
a level of about 10000 km. As we are interested in absorption at lower levels 
reflection is not important at the frequency concerned. However, waves of 
much lower frequencies will be reflected before reaching levels where their 
heating effect is important. 

7. Heating at higher levels.—'Vhe above discussion suggests that above solar 
yranules there may be an upward flow of electromagnetic energy of perhaps 
5x 10°ergcm™*sec"'. Averaged over the solar surface there may be a flow of 
5x 10°ergcm™*sec"! which penetrates to levels of 1o'km or more. ‘The 
region where heating is most important is the transition zone from mainly 
neutral to mainly ionized hydrogen centred at a level between 6000 km and 
7000 km. 

In the H1-»>Hit zone we have +r, =1 sec (see Appendix 1) and the 
proton and neutral atom densities are about 3 10%cm* so that »=1 and 
V’=1-4x10’cmsec"!. The absorption coefficient (equation (11)) is about 
10°*cm~! for waves of frequency o-2radsec™!. At levels where the hydrogen 
and helium are mainly ionized the absorption is negligible owing to the small 
value of 7. At levels where the gas is mainly unionized the absorption is much 
the same as in the H1-+-H zone but the heating effect is negligible because 
of the high gas density. ‘lhe total thickness of the heated layer, allowing fot 
some oblique propagation, is likely to be a few thousand km and the total energy 
absorbed about 10° ergem *sec '. ‘This is comparable with the estimate of 
energy needed to replenish the quiet or spot-free solar atmosphere, in which 
only a small proportion of the average coronal energy is released (14). 

The peak gas velocity due to a passing hydromagnetic wave is |’, /1, 
where H, is the maximum of the perturbation magnetic field. In the upper 
part of the transition zone V’ may exceed 2 x 10‘ cmsec ! and occasionally /1, 
may equal H, in magnitude. ‘lhe corresponding kinetic energy of a proton ts 
34x10 “ergs or ~200eV. If the proton collides with a neutral atom which 
does not share the plasma velocity then each particle receives energy of ordes 
100 eV, resulting in ionization and thermal motion corresponding to a temperature 
of order 1o*deg. K. Since the wave quarter-period (w= 0-2 radsec ') is ~8 sec 
and the collision period 7, > 1sec such collisions will be reasonably trequent. 
The process might be likened to the boiling away of a liquid (the dense unionized 
gas) to form an overlaying vapour. 

8. Discussion.—The new mechanism of absorption of hydromagnetic waves 
disposes of the main objection to the theory of solar heating by these waves. 
It seems likely that waves due to granules cause enough heating in the 
H11-» Hu zone (and the corresponding He zone) to account for the quiet 
corona. 

Corresponding effects on a much larger scale are likely to occur in the 
neighbourhood of sunspots where violent hydromagnetic disturbances are often 
observed. In extreme cases the heating mechanism may cause flares. 

A prediction of the theory is that most of the significant heating should occut 
in the surface layer of lightly ionized gas, that is near the surface of the visibli 
chromosphere or of prominence material. In this region there should be a very 
marked departure from conditions of local thermodynamic equilibrium 

The theory of energy transfer by hydromagnetic waves seems to hav 
advantages over corresponding theories involving purely hydrodynamic effects. 


>? 
at 
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A discussion of these theories of solar atmospheric heating, due to Biermann, 
Schwarzschild and others is given by van de Hulst (12). The main difficulty 
seems to be that energy is dissipated at too low a level in the chromosphere. 
‘The most recent suggestion, made by Schwarzschild, is that the pressure waves 
concerned might have initial material velocities of only 10-? or 10-* that of the 
granulation. ‘The energy of these waves would tend to be dissipated in the upper 
chromosphere and corona. ‘The obvious objections to this suggestion, which 
involves throwing away all but a tiny fraction of the initially available energy, 
are avoided when the extra-elastic forces introduced by a magnetic field are taken 
into account. 


APPENDIX 


The collision period of neutral atoms.—-A plasma of density N (NV electrons 
and N protons or other heavy ions) drifting with velocity v, through an atomic 
gas of density N, (= 9N) will experience a force Fem * retarding its motion. 


’ 


‘lo a reasonable approximation we may write 


N 
F Nv, ("4 ") (12) 
T, 


Tj 


where mand M are the electron and heavy ion masses and 7, and 7, their respective 
collision periods with neutral atoms. 

For weak hydromagnetic waves the values of 7, and 7, are found from the 
collision cross-sections at gas kinetic energies. ‘The cross-section of helium 
atoms for slow protons is about 2-6 in units of wa)? (Mott and Massey (15), 
p. 289). ‘The relative mobilities of various heavy ions in He and H, gases (16) 


are in the ratio about 2:1. Allowing for the decrease in mobility because the 
hydrogen was in the molecular form the collision cross-section of H atoms (for 
protons) is about 2-3 x 10 '*cm*. ‘The corresponding cross-section for electrons 
(ref. (15), Fig. 21) is about 3°6x 10cm. ‘The collision period of electrons 
is less than that of heavy ions because of the greater collision cross-section and 


the greater velocity. However, these effects are more than offset by the much 
smaller electron mass so that the term m/r, may be neglected in a reasonable 
approximation. 

From equations (2) and (12) we now have (remembering that v,=(v—v’) 
and NM =p) 

T, = WT; 

as would be expected when electron collisions are neglected. 

At the chromospheric transition level Hi -» Hit reasonable values of 
proton density and temperature are about 3 x 10°cm* and 1o' deg. K. The 
corresponding value of 7, is 0-9 séc. 
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ON ‘THE EQUILIBRIUM OF MAGNETIC STARS 


L. Mestel 


(Received 1956 April 14)’ 


‘The effect of a magnetic field on the equilibrium of a gaseous star is studied 
by a perturbation method. It is shown that the pressure-temperature field re- 
acts on the meridional magnetic field not through the condition of hydro- 
static support, as in a liquid or barytropic star, but through the energ: 
equation. ‘This reaction is in general very small, so that the imposed field 
may be considered as arbitrary, in contrast to Ferraro’s result for a liquid 
star (1). ‘lhe electron-pressure and magnetic force terms in the equation ot 
conduction lead on the average to a slight modification of the law of tse- 


rotation 


1. Introduction.—I\n a recent paper, Ferraro (1) has considered the mechanical 
equilibrium of an incompressible, uniformly rotating magnetic star; he shows 
that the equation of hydrostatic support implies that the magnetic field (assumed 
symmetric about the axis of rotation) must satisfy a certain non-linear condition. 
He derives the explicit solution for a dipole-type field, finding the dependence on 
radial distance to differ considerably from Cowling’s decaying field (2). He 
therefore concludes that Cowling’s field will not satisfy the conditions of equili- 
brium, and so must set up mass motions. 

Ferraro’s result can be generalized to state that if the density of the stellar 
material is strictly a function of pressure alone, then the star can remain in 
equilibrium only if the non-spherical perturbing fields (centrifugal and magnetic) 
satisfy jointly a certain complicated condition. ‘The physical reason is that 
polytropic equilibrium is neutral for purely internal displacements, so that the 
star lacks the extra degree of freedom necessary for it to withstand an arbitrary 


perturbing force; if one component is given, the pressure field adjusts itself in a 


way that severely limits the other component. 

However, it can be seen that this conclusion does not hold if the equation of 
state is that of a perfect gas, as in real stars; for then an arbitrary perturbing field 
can be balanced by non-spherical variations in temperature as well as in pressure. 
In a convectively stable region, the pressure field reacts on the meridional magnetic 
field not through the equation of support, but through the thermal equation. 
Through the bulk of the star, this reaction is small as long as the magnetic force 
is small compared with gravity, and so the imposed magnetic field may be con- 
sidered as essentially arbitrary. 

An important point of method is involved here. Normally in applied mathe- 
matics if a. problem is difficult, the model is simplified and the results treated as a 
first approximation to the truth. Such a procedure is justified only if we know 
on general grounds that the results obtained are qualitatively correct, although 
quantitatively approximate. For this reason, the assumption in the present 
problem of incompressibility or barytropy is misleading. 


Received in original form 1955 January 21. 
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‘The discussion is limited to fields symmetric about the axis of rotation. In 
Section 2 the equation of support is discussed for incompressible, barytropic 
and perfect gas stars respectively. In Sections 3-8 the mutual interaction of the 
magnetic, thermal and rotation fields is studied by a method of successive approxi- 
mation. 

2. The condition of mechanical equilibrium. —Consider a fluid mass ot density 
p (not necessarily uniform), in equilibrium under a magnetic field H_ maintained 
by currents j, and with an arbitrary rotation field Q. ‘The equation of hydro- 
static support is 


(i) pRNra(Qar) + pVd Vp oO, (1) 


( 


where ¢@ is the gravitational potential, p the gas pressure, r the position vector 
trom the centre, and c the velocity of light (Gaussian units being used). ‘The 
sufhx m selects the meridional component. 

On taking the curl of (1), we have 


jaH 
Voave+va(( ) p2n(Qrry ) Oo. 
c m 
If p is constant, (2) reduces to 
vai(/B) RrA(Qar)> =o; (3) 
pe m 


or in cylindrical polar coordinates (w, 4,2) referred to the centre and the rotation 
axis 


(jo) + =< (joe) > - ns 


— (H,)?4 po |Q}P=o, (4) 
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where use has been made of the equation VAH = 47j/c. If the magnetic field ts 
zero, (4) implies that in equilibrium | | must be a function only of w, the distance 
from the rotation axis. In a magnetic star with such a rotation law, condition (4) 
restricts the form of the magnetic field: Ferraro’s result 


(jy )+ <— (jvH,,) O 


follows if the field is assumed purely meridional. 

If p is not a constant, but is strictly a function of pressure only, then the same 
result (3) holds, as is seen by subtracting from (2) the vector product of (1) with 
Vp/p; but now p is determined by the equation of support (1) and Poisson's 
equation V2d 47Gp. In anon-magnetic star, the same condition 6 |Q |/dz =o 
restricts the rotation field: if H does not vanish, however, the constraint exerted 
on the rotation and magnetic fields is in general more complicated than (4), because 
of the dependence of p on the non-spherical perturbing forces. 

‘hus in both the incompressible and barytropic cases, the condition of mech- 
anical equilibrium exerts a constraint on the magnetic and rotation fields; if they 
do not in fact satisfy this condition (3), rapid mass motions are set up. As in a 
highly conducting fluid the field is dragged by the matter, while each element of 
fluid tends to conserve its angular momentum, both the rotation and magnet 
fields are distorted until they jointly satisfy (3). 

By contrast, consider a gaseous star with p ~p7’, 7 being the absolute tem- 
perature. Equations (1) and (2) can now be satisfied, for arbitrary perturbing 
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forces, by allowing non-spherical variations in 7 as well as in p and p. Ina con- 
vectively stable region, mechanical equilibrium is maintained at the expense of 
radiative, leading to slow material currents in meridional planes, which maintain 
thermal balance through the energy they transport (3). ‘The electromagnetic 
effects of flow across the field can then be determined in a higher approximation. 
In a convectively unstable region, strong interaction between the perturbation 
and the convection currents is likely. ‘Turbulent viscosity will tend to reduce 
non-uniformities of rotation; a weak magnetic field will be tangled up by the 
convection, while a stronger field will inhibit flow across the lines of force. But 
whatever modifications the perturbing fields undergo, mechanical equilibrium 
will be preserved by temperature variations, to which the convection currents 
adjust themselves. In the rest of the paper, only stable regions will be treated. 

However, in considering the interaction between a spherical gaseous star and 
4 perturbing rotation or magnetic field, a distinction must be drawn between 
fields that are imposed “rapidly” or “ gradually”. If changes in the fields take 
place in a time comparable with or longer than the Kelvin» Helmholtz time-scale, 
then the temperature field has the time to adjust itself to the slowly changing 
density and pressure fields so as to preserve mechanical equilibrium. ‘This ts the 
case, for example, if the star’s magnetic field is a distorted part of the local galactic 
field dragged in with the condensing proto-star; or if the field has been built up 
against self-induction in a time of the order of Cowling’s decay-time (2). Dis- 
placements of the elements of the star by the perturbing forces will react back on 
the rotation and magnetic fields, but provided the forces are small compared with 
gravity the distortion of the imposed fields will be small. 

On the other hand, suppose the rotation and magnetic fields are changing in 
a time short compared with the Kelvin-Helmholtz time-scale. ‘The resulting 
variations in the density of an element of matter are accompanied by adiabatic varia- 
tions in pressure, as there is no time for the temperature field to adjust itself so as 
to preserve mechanical equilibrium. Mass motions, with distortion: of the 
rotation and magnetic fields must take place, but it is difficult to draw any general 
conclusions about any quasi-equilibrium state that may be reached. In parti- 
cular, it is not obvious that the distorted magnetic field should differ only slightly 
from the imposed field, even if the magnetic pressure is small compared with the 
thermal pressure. However, except briefly in Section 8, we shall deal only with 
slowly varying fields. 

It should be emphasized that Ferraro’s result (5) and its generalization (3) do 


depend on the pressure being strictly independent of temperature. In a real 


liquid or degenerate body there will always be a small thermal contribution to 
pressure, and again it is possible to satisfy the condition of equilibrium (1) by 
allowing non-spherical temperature variations, whatever the shape of the per- 
turbing fields. In fact, equation (3) is just the condition that the perturbations 
to the temperature be of the same order as those in the pressure and density ; if 
(3) does not hold, larger temperature variations are required, with correspon- 
dingly larger thermal currents. ‘hus even for liquid bodies or white dwarfs, 
it is doubtful whether (3) represents a real limitation. 

3. Gaseous magnetic stars._We now assume that the star is in mechanical 
equilibrium under a meridional magnetic field H,, maintained by azimuthal 
currents j,. For the moment rotation is ignored. ‘hat perturbations to the 
field H, must arise can be seen from the equation that replaces Ohm’s law tn a 
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highly ionized gas (4) 
H VY jrH 
a0(es 2 Le, o 
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where o is the unreduced conductivity, E the electric field, p,, a, the electron 
pressure and density, and —e the electron charge. ‘lhe non-spherical pressure 
density field yields an effective e.m.f. in meridional planes, while thermally- 
driven material currents induce electric currents by flow across the field. ‘The 
following scheme of successive approximation therefore suggests itself : 

In the absence of a magnetic field H,, the electron pressure term Vp,/n,e 1s 
given by 

a a. 
(Z+1) e pm  (@+iye™ 

where again the sufhx zero refers to the unperturbed spherical star, and A and Z 
are the mean atomic weight and number respectively. ‘his force on each electron 
is just cancelled out by an electrostatic force 


Amy, 


(Z Fi 1)e Vdy. (3) 


E, 
This is the Pannekoek-Rosseland field, with a corresponding charge density 
V.E,/47e.s.u. The field is built up in the very short time 1/0, the same as the 
decay-time of charge density within a conductor without any non-electric force 
acting on the electrons. 
The currents j, maintaining the field H, are driven by an electric field E,, so 
that 
47, _ 47 
VaH, se = o,E,, 
tf c 
I OH, 


VaE,=--—, 
& ¢ 


V.H,=o. 


As usual in conductors, the displacement current is negligible. 
The magnetic force j,,H,/c in the equation of support involves H,* and so 
the perturbed quantities ¢, p, p, 7 are written 


hb by + by eee P=PotPoat--- 
p=Po+ Pr T=T,+T,+.... 


The perturbations satisfy 
V"4_= — 4nGp, 
from Poisson’s equation ; 
a ine a 
p2= i, (pols + pyTy) 


from the equation of state; and the perturbed equation of support 


. A H e 
(hath) + poV be + poVdby — VPs =0. 
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The electron pressure term in (6) now yields, besides (7), 
A Vp» V / jaH 
Caewe_ Ps a Po | wm lee 4 os | (14) 
(7 +1)e py Po CPy 


‘Thus the last two terms in (6) give electromagnetic effects of order H,?; we there- 


fore write 
‘. Amy ( j,rH, on (Bs 
Jo »(E, (Z 4 1)e Vd, ’ ~ Cy ) Le c CPy )) 


Amy Amy (4 A =) 
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 . 
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and also 


VaE, ! (16) 


V.H,=0. 

‘The term in V¢, and the irrotational part of the term in j,,H,/c merely change 
slightly the electrostatic field and the charge density, but in general the last term 
leads to non-vanishing electric currents in meridional planes, maintaining an 
azimuthal field. ‘The interaction of these currents with the meridional field H, 
and the rotation of the star (neglected so far) will be considered in Section 6. 

The velocity v of matter flow in meridional planes is given by the energy 
equation 


={ 1 OT Top RT TM» . 
yaa) ; *) + = v.¥(In —- ) (17) 


pu r p ot py p 
V.F : i” 

pP po 
where F is the energy flux, « the energy liberation per gram, and y the ratio of the 
two principal specific heats. In the absence of a magnetic field, the Joule heating 
term j,”/po vanishes ; all the quantities are spherically symmetric, and the energy 
equation determines the approach of the star to strict thermal equilibrium, for 
which «—V.F/p=o, in a time of order of the Kelvin-Helmholtz time-scale. 
With a perturbing magnetic field present, the radial velocity will vary over a sphere : 
the spherically symmetric part of the perturbation will expand or contract the 
sphere as a whole, again in a time of Kelvin-Helmholtz order, while the rest drives 
the meridional circulation, as in the rotational case (3). ‘The two components of 


circulation are of order H,*, and are related by the equation of continuity 


V .(pv ~ 2 oO. 8 

(pv) + a °) (1 ) 

‘The terms of order /,* in (6) yield, besides an extra meridional electric current, 
the azimuthal current 

j . (e ' v.sH, a2 (a “ 37s_ (19) 

3 0 3 Cc Zec Po 2 , ZA 1]? } 

the last term arising from the perturbation in the conductivity(a2c 7%*), Again 

E,, j, and Hy, are related by Maxwell’s equations (16) with suffix 3 replacing 2. 
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If terms of order 4 or higher in H, are required in (6), perturbations of the 
equation of support of second order in magnetic force must be computed. 

4. The field A,.—From (9) it follows that E, satisfies Cowling’s equation of 
decay (2) 
470, cE, 


Va(Va€E,) > > (20) 
Ce « 
With H, assumed purely meridional, in spherical polar coordinates (7, 4, 4) the 


eigen-solutions of (20) are azimuthal vectors of magnitude 
7p 


4 
rf Ar) ap ¢ “rng, (21) 


where f,, satisfies a second order linear equation. ‘lhe series of eigenvalues 
r,, are found from the conditions of finiteness at the origin and continuity at the 
stellar surface; however, the order of the decay-time of a field of scale / is seen 
from (20) to be ~470,/*/c?. For a dipole-type field, the slowest decaying com- 
ponent has a time-constant of about 10" years (2). At time ¢=o an arbitrary 
magnetic field can be expanded in terms of these eigen-solutions, thus fixing the 
field H, at all later times. ‘The higher order effects may then be considered term 
by term. 

5. The azimuthal field and the stellar rotation._-F rom (15) and (16), the vector 
potential A, of the field H, satisfies the supplementary relation 

V.A,=0 


and 


~Va(VaA,) = sree (Se + An (Ao). (22) 


c \ ot Z(Z+1)e\ py 

The bar over the last term indicates that the irrotational part has been subtracted. 
We ignore for the moment the decay of the term in //,; then, as shown by 
Cowling (4), equation (22) implies that a field is slowly built up against self- 
induction through the action of the impressed e.m.f. | Am) Z(Z + 1)ec}(j, 4 Ay/ po). 
The asymptotic azimuthal field H, satisfies 


P i 479% ( Amy in rH, . 
V H, _ c2 Z(Z hi S )( Py ) ° (23) 


This field H, can be expanded in terms of the azimuthal ecigen-solutions of the 
decay equation 


70, OA, 
-~Va(VaA,)= ee > (24) 
the time of build-up of each eigen-solution is the same as its time of decay, if it 
were imposed on the star with no quasi-electric forces to drive the maintaining 
currents. As in fact the driving term in(j, AH,) is itself slowly decaying, it is 
clear that the field H, would grow to a maximum, and then decay along with the 
field H,. 
However, in fact the star does not allow this azimuthal field Hi, to be built up. 
So far we have neglected the effect of rotation on the equation of conduction (6). 
If k is the unit vector in the positive azimuthal direction, then a rotation field © 
leads to an induced meridional e.m.f. 
Qdr sind 


-——(kaH,). 
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Let us ignore for the moment the terms in Vp, and j,AH,/c in equation (6). 
Then the assumption that H, be purely meridional implies that the meridional 
electric field should be electrostatic, so that no currents flow, maintaining an 
azimuthal field. ‘lhe condition for this is 


Qrsinf 
Vi ( 


(kH,) O, (26) 
which at once yields Ferraro’s well-known law of isorotation, restricting the 
ingular velocity to be uniform along a line of force. 

Now reintroduce the last two terms in (6), which we have seen tend to build 
up slowly a meridional current field. In the presence of the meridional field H,, 
the magnetic force has in general a non-vanishing component j,aH,/c in the 
azimuthal direction; as there are opposing it only the very small Coriolis and 
viscous forces, this force must affect the rotation field, and so also will modify j, 
through the motional induction term Qrsin@k~H,/c. ‘The change in the rota 
tion law will cease only when j,»H,=0; in general this implies j,=0. ‘The 
time-scale of adjustment is comparable with the period of a hydromagnetic wave, 
which for any sizeable magnetic field is much shorter than the characteristic time 
of growth or decay of the field. We may therefore expect the rotation field to 
suffer a gentle progressive modification with time; as the impressed meridional 
e.m.f. slowly overcomes self-induction and tries to build up a current field, the 
rotation field is immediately altered so as to cancel out the currents. ‘The maxi- 
mum deviation from isorotation occurs after the time in which the field H,, would 
have been built up, if the effect of rotation on (6) could have been ignored ; after 
this time, the impressed e.m.f. will have overcome self-induction, and the rotation 
field will satisfy 


Am lis | | 
Va (o yg to) ‘sin Ok ? 
A(Z + te pyr sin @ dosent , 


in place of (26), leading to 


Amy lint 


Q hs. 
~~  £(Z+1)e per sind 


constant along a line of force. (28) 
However, over the bulk of the star this correction to Ferraro’s law is small; e.g. 
if H, is Cowling’s dipole-type decaying field, (28) reduces to the constancy of 
({2— 10 “H/p,) along a line of force. As even a slowly rotating star like the sun 
has (2~3 x 10 © at the surface, it follows that the departure from isorotation is 
negligible in all cases of interest. 

6. The perturbation to the meridional field. —\n computing the current j, from 
(19), we may ignore the term in j,4H,/c, for we have seen that this azimuthal 
force merely serves so to modify the rotation field as to cancel itself out. ‘The 
vector potential of H, is then an azimuthal vector satisfying 


1% OA eT. 
Va(VaA,)= ye my ~ (v.04, + 35#e,)), (29) 
0 


again implying a slow build up of M1, against self-induction to a maximum 


VaH,- (vik, t : 7h). (30) 


and then a slow decay as the primary field H, decays. 
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‘The velocity v, is calculated as in the rotational case (3). From (11) and (13) 
we have an equation for 0¢,/00; dp,/00, dp,/06, OT,/08, and finally 0(v,),/0@ are 
then found from (11), (13), (12) and (17). ‘The condition of zero outflow from a 
sphere of given mass determines (v,), uniquely, and (v,), is found from the 
equation of continuity (18). 

The orders of magnitude of the different terms in (17) are as follows 


V.F i. j,aH, . 
( =) ~M\ = | / oa bo’) 


ROT, BT yep. BT y|j, rH, R |j,aH, 
= ees : es Pol dy \~ ; , 
je ot fh po ot “a t 
i,” oj, je |2 
oe ao en Pots 
Po%o 47 Lpyay 
where L, M, and R are stellar luminosity, mass and radius, 7 is the time of mag- 
netic decay, and / again the scale of the field. (If H, is a superposition of fields 
of widely varying dimensions, then we select the largest term). 
Consider first a large-scale field, with /~R. ‘Then inspection of (31) shows 
that the last two terms are comparable in magnitude, while the first is larger than 


Pot € 


rf 


the other two by a factor of about 200. Hence from (17) the velocities through 
the bulk of the star are of the order 

La, ai 

M1 gnR | Poh Fo (32) 
The term (}c7,/7,) | €,| in (1g) is small compared with |v,,H, |, and so the 
asymptotic value of | V~H,| is of order (470,/c?) |v, 0H, |, giving 


TT \32 hk 
As wT rT.) H*(4q): (33) 


H, (p p,) M 
(Barred quantities are expressed in solar units, and the sufhx ¢ refers to central 
values. 7°. is taken as 2 ~ 10’ deg. K, and p, from Miss Gardiner’s integration of 
a Cowling Model star (§)).  “Uhus through the bulk of the star only an enormous 
field H, can give rise to velocities large enough to modify seriously the meridional 
field. 

If/- Rthen H,/H, is greater than (33) by a factor (R//); with /- R/200 the 
last two terms in (31) become dominant and give an even larger perturbation. 
However, the decay-time of such a small-scale field is much shorter than the 
stellar lifetime, and unless there is operative some process that continually crinkles 
the large-scale field we may ignore it. 

‘The perturbation technique of this paper may break down near the edge of 
the star, as |j, »H, |= o,)| E,,H,/c| drops tozero less rapidly than po( — dy’). “Vhis 
leads to a large v,, while the equation of continuity (18) gives 7, even larger, by a 
factor rp,'/py. ‘Thus near the surface, with 1~R, 


d ( p ( T\32 
} 
H, dr od 7 ) (= 
aie ) 4 - 5 ' 
a (? ) mi u) 
p 
Putting in photospheric values for p and 7, we have 
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At the base of the Unséld convection zone the numerical factors will be much 
smaller. ‘hus the perturbation technique seems adequate for main-sequence 
stars with large-scale magnetic fields unless the fields are of the intensities found 
by Babcock (6). On the other hand, a locally concentrated strong magnetic field, 
of the sort that may give rise to sunspots cannot be treated as a perturbation, and 
an entirely new approach is necessary. 

If the field H, should be such that H,/H, ~1, then in fact the thermally driven 
matter currents would distort the lines of force so that they follow closely the flow ; 
near the surface the field would become more or less horizontal, because of the 
large horizontal component of flow. Perturbation technique could then be 
applied to this distorted field unless the magnetic force becomes comparable with 
the gravitational force. Further discussion of this case is beyond the scope of 
this paper. 

We may conclude, then, that the meridional magnetic field imposed on a star 
may be considered as essentially arbitrary, provided it is not too strong nor changes 
too rapidly in space. ‘he contrast with the situation in a liquid or barytropic 
star has already been pointed out. ‘lhe results are intuitively satisfactory: as 
long as the magnetic forces are small compared with gravity we expect them to 
have small consequences except when they cannot be balanced by small pertur- 
bations in gravity and gas pressure, as with azimuthal forces, or with general 
meridional forces in a barytropic star. 

Rotation... Yhe effect of rotation on the equation of conduction (6) has 
been dealt with in Section 5. As long as we are considering only first order 
perturbations in the equation of support (1), then the distortions of the pressure- 
density field due to magnetic and centrifugal forces may be computed separately 
and superposed. Again, the rotational distortion of the star will react back on 
the imposed fields, but as long as the ratio of centrifugal to gravitational force is 
small this perturbation will be small; in any case, the field H, introduced in 
Section 3 will be the field present in the magnetically and rotationally distorted 
star, apart from the perturbations Hy, etc. 

We have assumed so far that the field H, is purely meridional and have seen 
that this implies that the angular velocity is initially constant along a line of force, 
but suffers very slight progressively increasing non-uniformities. However, a 
rotation field conforming to this law will still in general yield a centrifugal field 
not derivable from a potential ; as first pointed out by Biermann (7), such a centri- 
fugal field drives meridional electric currents j,, that build up an azimuthal field 
H,, analogous to the j, and Hy, of Section 5. ‘The asymptotic field will be 


Cc 


VaH, = $7%e/( ay )(@r(Beor), (36) 


(7+ 1)e 


where again a bar denotes that the potential part has been subtracted out. How- 
ever, in the presence of the meridional field H,, as soon as currents j,, begin to 
flow, the azimuthal force j,,~ H,/c reacts on the rotation field; as in the case 


discussed in Section 5, we may expect a gentle progressive modification of the 


m 


rotation field sufficient to keep the magnetic field meridional. ‘The equation 


analogous to (27) reduces to 


CQT 2Amy CQL 2Amyc ; 
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Neglect of the first term in each bracket gives Ferraro’s law. ‘lhe deviation from 
isorotation is of the order of the ratio 


singe (8 (F 38 
(Z+1)e i) se A): (39) 


In order that centrifugal force should not be too great compared with gravity, 
{2 must not much exceed 10™* radians/sec. anywhere within the star; thus only 
with fields as small as 10~’ gauss does (38) approach unity. With such small 
fields the coupling between the rotation and magnetic fields is in any case negligible, 
as the period of a hydromagnetic wave along the lines of force is longer than the 
valactic age. 

Thus we see that in the presence of a sizeable meridional field, the electron 
pressure and magnetic force terms in (6) do not build a large azimuthal magnetic 
field, but simply lead to a slight non-uniformity of drift across the lines of force. 
This is a particular case of a result due to Sweet (8). 

Should the initial rotation field be non-isorotational, then the hydromagnetic 
interaction between the rotation and meridional magnetic fields will cause an 
oscillation about the mean state in which the field is meridional. ‘lhe period of 
the oscillation is that of a hydromagnetic wave, which for any sizeable field is 
much shorter than both the time of magnetic growth and decay and the Kelvin 
Helmholtz time-scale. ‘The comparatively rapid changes in the centrifugal force 
and in the meridional magnetic force will therefore lead to adiabatic changes in the 
density and temperature of each element of matter. As pointed out in Section 2, 
this perturbed density-temperature field will not in general satisfy the equation of 
support, and so further mass motion must occur, with rapid distortion of the 
centrifugal and magnetic fields and adiabatic temperature density variations, 
until a state of equilibrium is reached. 

The rotational thermal currents v_ discussed by Sweet (3) will also react on 
H, through the induction term vo, H,/c in (6). ‘This perturbation become: 
comparable with H, if 


rant a — 


r ct 4mol’ (39) 
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/ being the scale of the field. ‘This condition asserts that the time of flow across 
a distance of order / should be comparable with the decay time of the field, From 
Sweet’s results with o~107!’ e.s.u. 


Vg~=10 1% LR® M*)Q2?, 


4ralvg ~ 2( l ( LR ()2 
emo Be) ape): 4) 


Thus quite moderate angular velocities will drive rotational currents strong 
enough to distort the field H,, whereas we saw in Section 6 that the magnetically 
driven thermal currents have quite negligible effects unless H, is enormous. 


and 
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kither the lines of force will be made to follow the streamlines of the rotational 
currents, so that Vo AH, =o, or possibly the field H, may be so twisted as to exert 
much stronger mechanical forces, so that the effect of rotation on the thermal 
equation is cancelled out. ‘This result is relevant to the oblique rotator theory 
of magnetic variables, in that only slowly rotating stars will preserve an initial 
obliquity of the magnetic and rotation axes. 

8. The imposed field non-meridional.—So far we have assumed H, to have only 
meridional components. If now H, is taken as purely azimuthal, then exami- 
nation of the equations of Section 3 shows that all the perturbing terms H,, H,.. . 
are also azimuthal. As the rotation velocities are parallel to the field, there is no 
interaction with the rotation field. ‘The asymptotic value of the ratio j,/j, is of 
order ogf1,my/pyec, which is comparable with unity in a stellar envelope for quite 
moderate fields H,. 

However, it is unlikely that the field H, will have no meridional component, 
and so we can again expect interaction with the rotation field through the azimuthal 
component of magnetic force. In general, the field at each point will possess 
an azimuthal component which changes sign periodically as the local rotation 
field fluctuates. Again, the period of hydromagnetic waves along the lines of 
force will be much shorter, for any sizeable field, than the time of build-up of 
the fields H,, H,...., and so only the meridional field exerts a magnetic force 
having a permanent effect on the electromagnetic equation. A star with an 
initial purely meridional field and a rotation field not satisfying the modified iso- 
rotation law is a special case, in that the magnetic field fluctuates about a state in 
which it is simultaneously meridional at every point of the star. 
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Summary 
‘This paper is concerned with homogeneous and isotropic model-universe 

‘The two main alternative descriptions of such universes, (1) by means of a 
Riemannian space-time map characterised by the Robertson-Walker metric 
and (ii) as a kinematic equivalence, are compared. ‘The generating functior 
i(t) and signal-parameter « of the equivalence are related to the scale-functior 
R(t) and radial coordinate r of the corresponding metric Vlilne’s distance: 
is related to the other distances commonly used in cosmology, for which 
expansions in terms of proper distance are given. ‘lhe empirical determina- 
tion from nebular red-shift data of the functions R(t) and y(t) appropriate 
to the observable universe is discussed with reference to Whitrow’'s 
second-order shift-distance criterion and in the light of recently published 
observations \lso, new general shift-ratio formulae for kinemati 
equivalences are derived 


1. Introduction. —Vhere are two main techniques for the discussion o! 


cosmological problems including the construction of model-universes. One 
which I shall call the Riemannian technique, originated in Minkowski’s analysis 
of Einstein’s “ Special Theory of Relativity’? and depends on the use of a 
Riemannian space-time map. ‘This became the basic technique of Einstein's 
“General Relativity”. ‘The other is the technique first developed by Milne 
and Whitrow for ‘‘ Kinematic Relativity”, which I shall therefore call th: 
kinematic technique. Neither of these techniques, however, is restricted to th: 
particular physical theory to which it was first applied, a fact well illustrated by the 
Riemannian treatment of “ Kinematic Relativity” by A. G. Walker and others. 
The present paper will stress the fact (implicit in the work of Robertson, Walker, 
Heckmann and others, but still, perhaps, not generally appreciated) that both 
techniques are indeed only techniques and not theories. As such, they can be 
adapted to the same cosmological hypotheses, and then when applied to a given 
problem they should furnish equivalent results. 

Both techniques have been developed for the investigation of homogeneous 
and isotropic model-universes, although the Riemannian technique is certainly 
capable of incorporating basic hypotheses other than homogencity and isotropy ; 
the same would appear to be true of the kinematic technique, but this does not 
seem to have been investigated. ‘The present discussion will be restricted to 
model-universes of the former type. 

Thus, on the one hand, it has been shown by H. P. Robertson (1) and, 
independently, by A. G. Walker (2) that any cosmological model which satisfie 
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the hypotheses of homogeneity and isotropy admits ipso facto an invariant 
quadratic form, now generally known as the Robertson—-Walker metric, 


PEF eT all nt nonitichn (1) 
(1 +hkr*/4)* 

which characterizes the model as set out in Section 2 below. Here we need 
only note that R(t) is an arbitrary function of t, and k is either —1, 0, or 1. 
The Riemannian 4-space of which (1) is the metric is called the Riemannian 
map of the model and furnishes a kinematic background against which its 
properties can be discussed. Additional hypotheses are, however, needed 
before a particular model can be isolated, that is, before we can assign a 
particular functional form to R(t) and a particular value to k. It is at this 
stage only that the various cosmological theories postulating homogeneity and 
isotropy diverge. For example: (i) if the field equations of ‘‘ General Relativity” 
and certain other reasonable assumptions are imposed, some fifteen different 
models are isolated (3); (i) if Milne’s dimensional hypothesis is imposed, a 
single form of R(t) is isolated (2); (iii) the hypotheses of the steady-state theory 
of Bondi and Gold give rise to a unique model (4); (iv) so does the cosmological 
theory of Dirac (5). 

On the other hand, it has been shown by Milne and Whitrow (6) that the 
issumptions of homogeneity and isotropy imply certain well-defined relations 
between the operationally defined coordinates of the fundamental particles of 
the model. It will be sufhcient to recall at this stage that these relations are 
characterized by an arbitrary generating function ¢(t). Once a particular form 
of if is specified, the kinematic properties of the model are uniquely determined, 
except for the type of geometry appropriate to it: hyperbolic, spherical or 
Euclidean. (This property, which in the Riemannian map is governed by the 
value of k, does not seem to have been adequately discussed by the exponents 
of the operational technique. It is of importance even in spatially one-dimensional 
problems like that of the shift-distance law when comparison is made between 
theory and observation, but its effect does not appear until terms of the third 
order are considered.) ‘hus, in this technique also, when applied to homogeneous 
and isotropic models, the kinematic background is largely characterized by a 
single unspecified function, 4(t), which can only be determined by additional 
hypotheses, such as, for example, Milne’s dimensional hypothesis, the “ perfect 
cosmological principle”’ of the steady-state theory, etc. 

It would appear, therefore, that any homogeneous and isotropic model 
with given kinematic properties should possess a description in terms of either 
technique and that there should exist a relation between the functions R(t) 
and (t) appropriate to it. In general, formulae appertaining to the one 
technique should be translatable into formulae appertaining to the other and 
results obtained by the two techniques should be comparable. One would not 
expect to find a complete parallelism, for evidently the kinematic technique is 
of narrower scope, if only because in it there are difficulties in ascribing 


operational coordinates to particles beyond a possible “horizon of vision”’ and 
also because in three dimensions the method still lacks complete development. 
But within the narrower field of the kinematic technique it is evidently desirable 
to establish Riemannian equivalents. Undoubtedly, the most important step 
towards correlating the two techniques has been Robertson’s and Walker’s 
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derivations of the metric from the hypotheses of homogeneity and isotropy. 
‘This showed that these hypotheses lead either to the metric (1) or to the concept 
of the kinematic equivalence, depending solely on which technique is used. 
Walker also showed (7) how, by a suitable clock-regraduation, the model 
characterized by the metric (1) can be made formally identical with the 
Milne-Whitrow “ uniform motion” equivalence. Walker’s further investigations, 
however, were concerned with the Riemannian equivalents of Milne’s ‘‘ Kinematic 
Relativity’ theory (free paths, etc.) as distinct from the kinematic technique. 

In Milne’s later work, and in Walker’s, only two particular equivalences are 
considered—the “‘uniform motion” equivalence and the “stationary”’ equi- 
valence. ‘These are regarded as alternative descriptions, in terms of two different 
time-scales, of the unique physical model isolated by the dimensional hypothesis. 
Yet the kinematic technique can be applied to cosmological problems in its own 
right, apart from the dimensional hypothesis, just as the Riemannian technique 
can be applied without the restriction of Einstein’s field equations. For a recent 
example we may point to a paper by Whitrow and Randall (8) in which more 
general equivalences are examined. In connection with such investigations the 
following hitherto apparently neglected questions are of interest, and it is the 
primary object of the present paper to attempt to answer them: 

(i) How can one obtain the scale-function R(t) corresponding to an arbitrary 
yenerating function ¢(f) ? 

(ii) What is the relation of the coordinate r in the metric (1), specifying a 
particle on a given line through the origin, and the signal-parameter « which 
enters into the well-known Milne-Whitrow signal-function @(t)= yal "(t) (9) 
and equally specifies a particle on a given line? 

(11) What is the relation of Milne’s distance (“‘radar’’-distance) to the 
distances commonly used in the Riemannian technique ? 

(iv) Once a relation between R(t) and +/(t) is established and /(t) is found 
empirically by way of R(t), to what clock-graduation of the universe does #/(f) 
correspond ? 

The last is a question of physics, whereas the others are purely mathematical. 
Question (iii) is dealt with in Section 2, the remainder in Section 3, which may 
be read independently of Section 2. 

\ further object of the present paper is an examination of kinematic equi- 
valences and Riemannian models (which can now be undertaken concurrently) 
in the light of the shift-distance law with particular reference to a recent paper 
by G. J. Whitrow (10). 

All expanding models are characterized by a Doppler red-shift which, to a 


first approximation, is proportional to distance. ‘he factor of proportionality 


can usually be made to agree with observation by adjusting a parameter. For 
many well-known models, however, the second-order term in the ¢ xpansion of 
the shift-ratio in the form 


ax + bx? +..., 


where x denotes distance,* cannot be brought into line with observation by the 


* It is a matter of choice which of the several possible cosmological “ distances’”’ is to be used in 
this expansion. ‘The values of the coefficients (with the exception of the first) will, in general, 


I 


depend on this choice, 


24 
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mere adjustment of a parameter. ‘The second-order term could thus be used to 
eliminate such a model. ‘The accurate empirical determination of the shift- 
distance law is, therefore, of great importance in this connection. ‘The only 
results available until very recently were those obtained by Hubble before 1937. 
In view of their unavoidable uncertainties no very definite conclusions could be 
drawn from the second-order term. Nevertheless, cautious use has been made 
of that term by Heckmann (11) im conjunction with data obtained from nebular 
counts, to determine some of the parameters of ‘‘ General Relativity’? models. 
A similar procedure seems to have been adopted by most other writers on the 
subject. 

Recently, however, in the above-mentioned paper, G. J. Whitrow re-examined 
the direct use of the second-order term, which is of particular interest now that 
far more accurate observational results are at our disposal (12). He considers 
not the absolute values of the coefficients but the ratio b/a*, which we shall denote 
by 7, of the second-order term to the square of the first-order term in the 
expansion of oA/A. Whitrow shows, for example, that 7=} for any non- 
recessional theory of the red-shift, provided that the correlation between shift 
and distance is independent of time and provided also that the expansion of 
bA/A is made in powers of an additive distance. ‘Thus the second-order term 
could be used to eliminate this whole class of models. Another advantage of 
using 7, which has not been stressed, is its independence of any scale error made 
in the distance measurements. ‘his is of particular interest in view of the recent 
drastic revision of the cosmological distance-scale, by a factor* of about 2 or 3. 
Several questions concerning 7 are discussed in Whitrow’s paper, but four, in 
particular, are left open and these will be considered in the present paper: 


(v) Which of the distances commonly used in cosmology are additive in a 
static model ? 

(vi) In the general case, how does the particular definition of distance used 
affect the value of » ? 

(vii) What is the value of 4 for the most general kinematic equivalence and 
Riemannian model ? 

(viii) What conclusions can be drawn from a comparison of the theoretical 
values of » with empirical evidence ? 

The answers to questions (v) and (vi) are contained in Section 2, where 
expansions of the shift-ratio and auxiliary distances are derived in terms of 
proper distance, the latter being additive in static models. Various expansions 
relating shift-ratio, luminosity distance, distance from apparent size and the 
coordinate r in (1) are well known (13), but expansions in terms of proper 
distance do not seem to be available. ‘The incorporation of Milne’s distance 
into these relations also seems to be new. Regarding question (vi) we may say 
that the first coefficient in the expansion of 5A/A, usually called Hubble's constant 
after multiplication by the speed of light, is independent of the particular 
definition of distance used, since all current cosmological “ distances’”’ are equal 


to a first approximation. But 7 is more sensitive to the exact definition of distance 
than was made evident in Whitrow’s paper. For example, it will be seen below 
(Section 4) that the value of 7 can be changed by as much as 2 when proper 


* It is not unlikely, however, that the correction for the more distant nebulae will prove to be 
more complicated than a mere scale adjustment. 
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distance (implicit in Whitrow’s paper) is replaced by luminosity distance in the 
expansion. 

Question (vii) can be readily answered for the Riemannian models (see 
Section 2), as general shift-distance formulae are available in this technique. 
By the correlation of techniques the Riemannian shift-distance formulae can 
then at once be made applicable also to kinematic equivalences. For the latter 
no such general formulae have previously been given. In the Appendix it is 
shown how these shift-distance formulae can be obtained, alternatively, by 
purely kinematic methods. ‘This serves as a useful illustration of the parallelism 
between the two techniques. 

It was question (viii) that, in fact, first suggested the correlation of the two 
techniques. For neither of the two main links between theory and observation 

the relations between luminosity and red-shift and between luminosity and 
nebular density——can be satisfactorily treated by the kinematic technique in the 
absence of a geometrical method. Now, by the correlation effected below, all 
the relevant formulae of the Riemannian technique can be made available to 
the models of the kinematic technique also. Section 4 deals with the general 
problem of comparing theory with observation, as far as the relation between 
luminosity and red-shift is concerned. 

2. Distance relations.—\n the present section | obtain expansions in terms 
of proper distance of the shift-ratio, luminosity distance, distance from apparent 
size and Milne’s ‘“‘radar’’- distance. Proper distance, /, as defined below, is 
the only objective or cosmic definition of distance and is therefore suitable as 
the fundamental variable in a theoretical discussion such as this. 

We must now make use of the Robertson- Walker metric (1) which, as already 
mentioned in Section 1, has been shown to be applicable to al/ theories based 
on a homogeneous and isotropic substratum. It has the following significance : 
(i) ¢ is a cosmic time-coordinate; (11) 4, ¢ are the usual angular measurements 
made at the spatial origin r=o (which can be identified with any fundamental 
particle); (iii) the world-lines of the fundamental particles are the geodesics r, 0, 
¢@=constant, whence r is a “‘co-moving”’ radial coordinate; (iv) light-tracks 
correspond to the null-geodesics of the metric, and, in particular, light-tracks 
through the spatial origin have the equations 6, ¢=constant and 

dt dr 

RQ) ~ * a+ kr*]4)’ (2) 
It should be noted that at this stage the constant ¢ occurring in (1) need not be 
identified with the velocity of light. Moreover, in deriving (1) no assumption 
need be made about the invariance of this velocity nor about the measurement 
of length (see ref. (1), p. 298). In any theory, however, which adopts the 
convention of the constancy of the local speed of light (and which consequently 
has related time- and distance-scales) that velocity can be, and usually is, 
identified with the constant ¢ in (1). 

The proper distance, /, between the origin and a fundamental particle with 
coordinate 7,, at time ¢t,, is defined (14) by the formula 


%. 
1=R() | (3) 


o 1 +hr*/4 
In theories of the type just mentioned this distance can be interpreted as the 
sum of the infinitesimal distance measurements made at the world-moment 1, 


24° 
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by a chain of fundamental observers situated along the space-geodesic joining 
the particle to the observer. In other theories / must be regarded as conventional. 
As a preliminary to obtaining the expansions of the auxiliary distances, as 
well as for the later discussion of the coefficient 7, we derive an expansion of 
the shift-ratio of the light received at the origin from a fundamental particle. 
Assuming for the moment that this is a pure Doppler-shift, one easily obtains 
from (2) the familiar formula (15) 
D= A+OA Rt) (4) 
A Rit.) 
where A and A+ 6A are the emitted and observed wave-lengths and ¢, and ¢, are 
the times of emission and reception, respectively, of the light-signal in question. 
Writing ¢,=t,—7 and expanding, we obtain for z= D—1 the series 


oA R’ {Rr 1k’ tn” RR iR 
. » fo + 73 ( — ——— + - —} +.,.., (5) 
A R (Fa 2 R) (7 Rk? 6 z) 
where we have written R, R’, etc., for R(t,), R’(t)), etc. ‘To this notation we 
shall adhere throughout this paper. Integration of (2) along the track of the 
signal and reference to (3) now gives 
& dt ry” ad / 
l. >.) or | 1+kr* 4 ¢R(tg—T)° (0) 
On introducing a new function s(t), defined by the equation s‘(t)=1/R(t), the 
integral on the left becomes s(t,)—s(t)—7), which can at once be expanded in 
powers of 7, and the same is true of the denominator on the right. In this way 


we obtain the series 
l rR’ 7/R” 1R" 
( 2R° “(F sf) te 
12R’ 1F/f R® R 
ro ++ oe = (2 x z) beeey 
which, when substituted in (5), gives the required formula 


+n. oF ( ow? 5. see R’R' RR" 
2 ck t 2¢2 3 Re R) + (10 7 ti + or) ore (9) 


Inverting, we get 


We now turn to the auxiliary distances. ‘The distance from apparent size, S, 
of a given nebula, is defined by the equation 
Ws 55 

S2? 
where A is the supposedly known cross-sectional area of the nebula in question 
and w is the solid angle subtended by it at the observer. ‘Taking the observer 
to be at y=o and the nebula at r=7,, at t=¢,, it has been shown (16) that 
R(t,)ry 
——— « (10) 
1+ kr,*/4 


‘Vhe distance from apparent brightness, or luminosity distance, L, is defined by 
the equation 


B 

L*’ 
where B is the supposedly known absolute brightness of the nebula and + its 
observed brightness. Again, it has been shown (17) that, if ¢,; and ft) are the 


b 
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times of emission and reception of the light by means of which the nebula at 
y=r, is observed, then 
RX(to) iad (11) 
R(t,) 1 + kr,?/4° 
Comparison of (10) and (11) gives, by reference to (4), the well-known 
relation 
L = DS. (12) 


On the other hand, comparison of (10) and (3) gives 
; oe. ' 
S= R(t) sin Ru)’ /, R(t,) sinh 
(according as k=1, 0, or —1) and consequently 
ki 
ORE, ) 


/ 
R(t,) 


S=/ 
This can be written as 
ki 
l- —, +..., I 
6k? ( 4) 
since R= R(t,) = R(t,) + O(/). Substituting (14) and (g) in (12), we also obtain 
the following series expansion for L: 
a a k 
L=l1+P— p(s se ~ 7 eee (15) 
cR CR ck OR? ‘ 
Finally, Milne’s distance, which is the fundamental distance in ‘‘ Kinematic 
Relativity”? and which will here be denoted by M, is defined by the equation 
M = he(t, —t,), 
where ¢t, and ¢, are the instants of emission and reception (after reflection) 
respectively of a light-signal emitted by the observer and reflected (at time ¢,, 
say) at the particle whose distance is required. Writing t,=t)—7, ty=ty-o 
and integrating (2) once more along the track of the signals, we see that 
[" dt , dt 
R(t) ‘ R(t)’ 


whence 
. dt : dt 
I , R(t) . P R(t)° 
ikxpanding both integrals by the method already employed in the derivation 
of (7), we obtain 
o* R o* 2R’*—RR rT? R r® 2R’*— RR 
R“ 2R" 6 R® tents a(j Rt” R’ 


; 
vhich, when solved for a, gives 


oO M m c Re 


n,n » BR 


Substituting for 7 from (8) we get the following series expansion for M 


Ph (RR 
N =| + ( + + 6 
iml+ Cal Re R) ia (10) 


‘This shows that as far as terms of the first and second order are concerned, 
M and /are interchangeable, which is of particular practical interest in connection 
with formula (15). (Compare also formulae (g) and (37).) 
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With reference to the recent paper by Whitrow (10) mentioned in the 
Introduction, we note that cr is clearly a linearly additive “ distance”’ in the 
sense employed by Whitrow (i.e. such that, if A, B, C are events along a 
light-track and in that order, then distance AC = distance AB+distance BC). 
It follows that the other distances, /, S, L, M, being non-linear functions of ez 
(see (7), (14), (15) and (16)), cannot, in general, be linearly additive in this sense. 
Only when the model is static, so that R(t)=constant, do we have cr=/= M. 
Even in a static model, the distances L and S are not necessarily additive. Only 
when the model, in addition, is spatially flat (k=0) do we have cr=l=M=S. 

3. The metric of the b-equivalence.—\t has been shown by A. G. Walker (7) 
that, whatever be the form of the function R(t) in the metric (1), it is always 
possible to find a regraduation from the universal time ¢ to another universal 
time ¢’ relative to which the model characterized by (1) is equivalent to the 
“uniform motion’’ equivalence of ‘‘ Kinematic Relativity’. ‘lhe required 
regraduation is given by 

cdt = 
7} R(t)’ (17) 


and it can easily be verified that this leads to the linear distance-epoch relation 
M/T =ctanho(r,) for a fundamental particle at r=r, relative to the origin r =o, 
where M and T are Milne’s conventional distance- and epoch-coordinates 
assigned in the usual way to the particle in question by means of light-signals 
and using ¢’-time, and where 


dr 


a(r;) : + hr®/4° 


Since the origin can be chosen so as to coincide with any fundamental particle, 
it follows that we have here exactly the situation of the “uniform motion” 
equivalence. 

On the other hand, a well-known formula of ‘ Kinematic Relativity”’ (18), 
sometimes referred to as expressing the ‘‘ main theorem on equivalences’’, gives 
the regraduation necessary to convert the equivalence generated by the function ¢ 
into that generated by any other function #. If’ is the time of the d-equivalence, 
then ¢’, the time of the #-equivalence, is given by 

= p(p (t’))’), 
where p and s are arbitrary constants which, tor our purposes, may be taken as 
unity. ‘The uniform motion equivalence is generated by the function 4(t) = ¢ (19). 
Thus a transformation from uniform motion time ¢ to the time ¢” of the 
#-equivalence is given by 

t” y(t’). 
The last formula in conjunction with (17) shows that the regraduation 
cdt 

OP | R(t) { 
will convert the model characterized by the metric (1) into the -equivalence. 
If (1) is itself the metric of the #-equivalence, we must have ¢=t" and therefore 


us 


1 

R(t) 

‘This formula gives the scale function A(t) corresponding to the generating 
function /(t), and vice versa. 


1d 
rls | s 
—- log ys *(t). (18) 
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An important function in the kinematic technique is the function «(?) given 
by the infinitesimal member of the signal function group, ¢ + ew(t). ‘This function 
is related to the generating function (ft) by the equation {(1/a(t)) dt = logy '(t) 
(20). In terms of this function relation (18) takes the very simple form* 


R(t) c(t). (1g) 


A further point of interest is the relation between the “co-moving” 
coordinate r in (1), which permanently appertains to a given fundamental 
particle, and the signal-parameter x occurring in the Milne-Whitrow signal 
function (9g) O(t)=%aup(t), which also permanently appertains to a given 
fundamental particle. From standard formulae of ‘‘ Kinematic Relativity’ (20) 
we obtain the relauon 


‘to dt 
ot 
t, and ¢, being, as before, the times of emission and reception of a light-signal 
from a particle characterized by the parameter z, to the observer. On the other 
hand, using (1g) and (2), we obtain 
te dt ‘te cdt ' db 
| att}, RO) |, rR 


log %,, 


a(’;), 
whence we get the required relation, 
z=expo(r), 


having discarded the particularizing suffix 1. 

The most important application of the correlation here established is to 
enable us to associate with kinematic equivalences a method for studying 
questions involving space-geometry. [t appears that only with the help of the 
corresponding metric can such topics as, for example, luminosity distance, 
distance from apparent size and nebular counts be satisfactorily discussed for 
a given equivalence. 

Conversely, the correlation throws certain light on the Riemannian technique. 
Milne’s concept of clock-equivalence, his recognition of alternative time-scales 
and of the connection between time-scales and distance-scales (depending on a 
conventional speed of light) are relevant to all cosmology. ‘he fact that an 
empirical determination of R(t) automatically decermines /(1), and our knowledge 
that (t) depends on the clocks that define ¢, force us to examine more closely 
the time-scale implicit in our observational work. ‘Vhis is somewhat hidden, 
since in practice none of the most significant measurements are actually made 
with clocks, but rather with photographic plates, spectroscopes, etc. ‘The 
time-scale enters only when these primary measurements are translated into 
theoretical terms, usually on the assumption that the observed nebular red-shift 
is a pure Doppler effect. 

If we assume that the red-shift can be regarded as a pure Doppler ettect, 
then when we so regard it we are automatically committed to the uniform time 
of “atomic vibrators’’, which I shall call ¢-time (not to be confused with Milne’s 
dual use of this symbol to signify also the time of the uniform motion equivalence). 
For in any other time-scale, in which the frequency of atomic vibrators changes 
secularly, this change constitutes a non-recessional effect contributing to the 
spectral shift. ‘hus it appears that f-time is generally implicit in our work 


| owe this observation to Dr (. J. Whitrow 
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and the generating function 4(t) found empirically by way of R(t) will therefor« 
refer to /-time 

‘The generating function appropriate to the observable universe and a given 
physical time-scale is, therefore, seen to be amenable to observational determination 
or, at least, verification. Milne’s identification of f-time with the time of the 
uniformly expanding equivalence, generated by #/(¢)=¢, can thus be put to th 
test. Since this identification depended on arguments essentially based on the 
dimensional hypothesis, we see that the justification of this hypothesis, which 
is perhaps the cornerstone of Milne’s theory, largely depends on the empirical 
confirmation, subject to the assumptions of ‘‘ Kinematic Relativity”, of this 
particular form of i. 

We conclude this section with a table showing the scale-functions R(t) 
appropriate to some of the more familiar forms of kinematic equivalences. We 
find that several of these correspond to well-known Riemannian models, provided 
the curvature index k is suitably chosen. ‘The last column of entries shows the 
value of Whitrow’s coefficient 7 =3/2— R’R/2R®. ‘This will be of significance 
in the following section. 

"Tasie | 


| 
Kinematic equivalence R(t) | Ruemannmian model n 





stationary equivalence Einstein’s static model | (4)(10) 


uniform motion equiy (a) one of the ‘ empty 

alence models of General 
Relativity 

(b) Milne’s 4 hydro- 


dynamical model 


Page's uniformly accel 


erated equivalence (22) 


Whitrow and Randall's for n=% : Einstein-de 
equivalence (8) Sitter model 
for? } Dirac’s model 
(5) 
Whitrow equivalence (a) de Sitter’s model 
e exp 
(23) é 





1 | (b) steady state model 


(4) 

















4. The second-order term in the expansion of the shift-ratio, and its comparison 
with the empirical data.Formula (g), for the shift-ratio interpreted as a pure 
Doppler effect, may be re-written in the form 


/ \2 


l l 
i. + 7) (;) Tees (20) 
cR/R' and 


3 RR 
=~ sR?’ (21) 


where, provided R’ + 0, | 


rT) 


‘The coetticient 7 in the expansion (20) is, in some sense, characteristic of the 


Various cosmological models. In ‘Table | of the preceding section the last 
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column gives this coefficient for some of the better known kinematic equivalences 
and Riemannian models evaluated according to (21) except for the first entry. 

If, on the other hand, we consider a static model and a red-shift consequently 
due to some non-Doppler phenomenon, we shall still have a relation of the 
form (20), although 7 will not now be given by (21). However, each theory 
postulated for the explanation of the observed red-shift must predict its own 
value for 7. Whitrow (10) has proved the general. result that if such a theory 
vives a time-independent shift-ratio then » = 4.* 

We cannot, of course, rule out the possibility that both the expansion of 
the universe and some other phenomenon contribute to the red-shift, but we 
shall ignore it here. 

We now turn to the observational data concerning the shift-ratio. ‘This 
ratio, 2, 1s correlated by observers with the photographic magnitude of the 
nebulae, and from these we can calculate the corresponding luminosity distances 
L. If we assume that 

s=aL+BL?+..., (22 
we can attempt to find the coefficients « and £ by “fitting the best quadratic 
curve”’ to the data. Implicit in this procedure is the assumption that the infinite 
series is fairly represented by its first two terms, at least for the distances so far 
surveyed. We can obtain some estimate of the validity of this assumption by 
examining one or two theoretical models. From calculations depending on 
formula (6) we can find the theoretical value of ///, corresponding to a given 
shift-ratio z in a given model, which may then be substituted in (20) to exhibit 
the pattern of the series. ‘l’o take two widely different cases, consider the models 
characterized by (a) R(t)=e"", k =o (e.g. the steady-state model) and (h) R(t) =¢, 
k=o (a “uniformly expanding” model). For (a) we find //l,=2/(1 + 2) and 
for (6) U/lj=In (1 +2)/(1+2). ‘Taking z=0°25 we find that the series (20) 
takes the following forms in these cases: 

(a) 0°25 =0°200 + 0'040 + 0'008 +..., 

(b) 0-25 =0°178 + 0°048 + 0°010+.... 
From (12) and (13) we have that L =(1 + 2)*/ in our cases (k=0), whence we 
find that the series (22) takes the following forms: 

(a) 0°25 =0°312— 0°097 + 0-06! 

(b) 0°25 = 0'278 ~0'039+0°014 

We therefore see that for these models and a shift-ratio of o-25 the second-order 
term in (22) 1s about two or three times as large as all the remaining terms put 
together. For smaller shift-ratios the convergence of course improve and 
the value z=0-2>5 may be regarded as about the frontier of present observations 
‘The maximum error in the empirical value of 2 in (22), obtained by “ fitting 
the best quadratic” to the data, may, therefore, be of the order of 50 per cent 
so that we can have some confidence in at least the sign of 2 so obtained, provided 
that additional sources of error are not excessive. 

Let us now rewrite (22) in the form 


LL‘ 
7 +0'(3 ) ; (23) 
where Lg=1/« and 7'=f/a*. Then if we know ¢ and f we can calculate 7’. 


* Whitrow’s result here applies, since it refers to the expansion of © in powers of an “‘ additive’ 
distance, which / is in the case of a static model See the remark at the end of Section 2 
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In order to relate 7’ to the coefficient 4 of (20) and (21) let us set 

L= DS. (24) 
It was seen in Section 2 that p=2 on a purely Doppler interpretation of the 
red-shift and on the assumption of constant nebular brightness. However, in 
a static model in which some non-Doppler agency causes the red-shift, the 
absence of the so-called number effect (or recession effect) in the calculations 
for L can be shown to lead to the value p=} (24). In “ Kinematic Relativity’, 
on the other hand, the variation of the absolute brightness of nebulae (Bat ') 
implied by the dimensional hypothesis (21) leads to p = 3, as can easily be verified. 
Only if we allow for the possibility, mentioned earlier, of a red-shift caused by 
more than one agency, would the relation between L and S be of an altogether 
more complicated type than (24). 

Referring to (14) we can rewrite (24) as 


L= Del+ O(F), 
whence, using (20), 
“a 
l4 pj { O(1*), 
0 
which is a general relation analogous to (15). When this is substituted in (23) 
we obtain an expression which must be identical with (20). Comparing 
coefhicients we then obtain /, = /,, and 
n=’ +p. (25) 
‘This relation allows us to compare our theoretical results with observation. 
‘There are two published sets of figures available for analysis. ‘The first set 
was obtained by Hubble before 1937 and analysed by Heckmann (25). It 
comprised data concerning ten nebulae with shift-ratios ranging from 0-004! 


to 01403. With due stress on the uncertainties both of the data and of the 


series expansion method as such, Heckmann gave the following values for the 
cocthcients in (22) (the unit of distance being 10° parsec): 
(1) z=O°178g+0°0160, f 00071 + 00160 
(26) 
(11) %=0°1782 + 0°0035, } 0'0150 + 0°0035 
the second set being obtained by omitting from the analysis the three most 
doubtful nebulae. ‘he following are the corresponding values of 7’ in (23 
(1) 0°22+0°50, (iI) 7’ O747+0°11. 
\lthough the distance scale has been drastically revised since these results were 
published, that portion of the correction which is merely multiplicative will 
affect only the values of and 8 but not of »’. I have, therefore, included the 
above figures for purposes of comparison with the new ones, 

Hubble’s data are now superseded by the publication, by Humason, Mayall 
and Sandage (12), of new data obtained with the most modern apparatus and 
corrected in accordance with the most refined theoretical techniques. ‘hese 
authors, moreover, carefully analysed their data-—amongst other things—for much 
the same purpose as ours, namely to determine second-order terms. ‘lhe 
publication of the present alternative analysis would, perhaps, be unjustified 
if its results completely agreed with those of Humason, Mayall and Sandage. 
‘The fact, however, that there is a difference (significant, though not large enough 
to affect the main conclusion) seems to indicate some uncertainty in both results. 
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‘The discrepancy can be traced to the use of different variables in the fitting of 


theoretical curves to the empirical data, and the sensitivity of the ‘‘ least squares ”’ 
method to such differences. Whereas these authors seek to determine the 
coefficients in the relation m= 5logcez + Bz +D, | seek x and f in the relation 
e=aL+ BAL. I shall comment on our differences below. 

The relevant data are contained in ‘Table XIII of the above-mentioned paper 
(12). We begin by converting the (corrected) photographic magnitudes, P., 
into luminosity distances, L, by the usual formula log L =0-2(P,— My) +1, 
where for the absolute magnitude MV, we take the authors’ value of 19°78. 
This leads to the following correlation between the shift-ratio, z, and L (in 
10° parsec). 


"Tasee Il 





I 
1°O14 


‘Then, fitting the best quadratic curve to these data by the method of least squares, 
| tind that 
z=0°0472L, + 0-0042L*. (28) 


‘The first term in this quadratic would imply a Hubble coefficient of 142 km/sec 
per mepaparsec, 

\lthough we have here neglected the tricky question of the magnitude of 
the possible errors involved, the convex trend of the z L graph and therefore 
the positive sign of the quadratic term in (28) seem beyond doubt. From (28) 
and (23), we obtain, in fact, 

7) E*@, (29) 

Some months ago Professor G. C. McVittie kindly informed me that he had 
computed the values of R’/R and R"/R from data on red-shifts available in 1954 
July and that these implied a value of 7’ lying between 1 and 3, which is 
compatible with (29). 

The values (27) of 7’ which result from Heckmann’s analysis of Hubble’s 
data would present fewer theoretical difficulties than the new value (29). By (25) 
the former permit the ranges 

(1) 0°22 <9 <0°75, (11) O'OS 4) O'l4 
for static models with time-independent red-shift (p iF 

(1) O78 <y- 1°78, (11) O-GQ2> y II4g 
tor the model of ** Kinematic Relativity ”’ (p -) and 

(1) 1°28 <y < 2°28, (11) 1°42 104 
for other expanding models (p= 2). 

hus the stricter sets of figures (11) would rule out static models with time- 
independent red-shift (y= 4) and the model of ‘* Kinematic Relativity” (4 = 4). 
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Both (i) and (ii) would rule out, amongst others, the model of the steady-state 
theory (y= 1). A uniformly expanding model (R= t, 7 = 5) without the particular 
assumptions of “‘ Kinematic Relativity’’ regarding intrinsic brightness would 
appear to be well favoured by Heckmann’s figures. 

‘The new data indicate that the value of 7’ implicit in Hubble’s data was too 
low. ‘The former lead to a convex 2 L graph (viewed from the axis of L), whereas 
Hubble’s data led to a slightly concave graph. If, for the moment, we accept 
no other consequence of the new data than this, 1.e. that 4’ >o, we find from (25) 
that 7 2 for the ordinary expanding model (p=2). Reference to (21) shows 
that this implies R” ~ 0, and so we are led to a class of models which, until recently, 
were rejected because of the well-known time-scale difficulty inherent in them 
(26). (Owing to a more rapid expansion in the past they are of even more recent 
origin than models with uniform expansion at the present rate, and the latter 
led to a low age of the universe not consistent with geological and other data.) 
The newly revised scale resolves this difficulty to some extent. However, before 
we discard, on the basis of criterion (2g), practically all the models that have 
been put forward in the last two decades* we must carefully consider the 
following possibilities: (1) the shift-ratios of the absorption-line spectra of the 
most distant nebulae may have been systematically increased by some non- 


Doppler agency; (ii) the distances of the remote nebulae may be systematically 


underestimated, owing perhaps to (iii) the possible non-constancy of absolute 
nebular brightness in ¢-time. ‘Iwo further possible sources of systematic distanc: 
errors are (iv) the Stebbins Whitford effect and (v) internebular absorption. 
Both these, however, make for overestimates rather than underestimates of 
distance. 

Humason, Mayall and Sandage find that the coefficient B in their relation 
m~slogez + Bz + Dis —1:180+0°875. From the theoretical expression 


RR’ 
B 1-086 (1 —, n), 
R’? 


where y is a correction term depending on the factors (i), (iv) and (v) above 
and which we shall neglect for the moment, this is found to imply 

RR 

Re? 


whence 4’ =0'5 + 0-4 for an ordinary expanding model. When an “aperture 


+ or8 


correction” is applied, these authors find 

Ri 

Ie 
whence 7 1'o+0°4. ‘This is considerably lower than our value of 1-g. Even 
so, the authors consider that on present theories of nebular evolution and with 
the present data concerning the Stebbins Whitford effect, the term jp cannot be 
made sufficiently large for the tentative conclusion R” —o to be avoided. 
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APPENDIX 


Shift-distance formulae obtained by the kinematic technique 


‘The two results of this Appendix, viz. formulae (33) and (37), could be 
written down at once by translating the Riemannian formulae (4) and (g) into 
their kinematic equivalents according to Sections 2 and 3. I here obtain these 
formulae by purely kinematic methods as an example of the parallelism between 
the two techniques. Such general formulae for the shift-distance relation in 
the kinematic technique do not appear to have been obtained previously. 

Consider two fundamental particles, A and B, belonging to a kinematic 
equivalence and let 6(t) be the signal function appertaining to these two particles. 
Let an observer attached to A emit a light-signal at time ¢,, which is received 
at B at time ¢,, whereupon the observer attached to B instantly emits a signal 
which arrives at A at time fy, all times being measured by the local clock, so that 

ty Ht»), ly Wt,). (30) 
A second signal is then emitted at B at time ¢, + dt, and received at A at time 
ty+dty. If dt, is the duration at # of a single pulse of the emitted light, then 
the Doppler-shift in this light when received at A is given by 
A+ 0A d 


lo ; 
D 5 di, 0’(t,), (31) 


where A and A+ 6A are the wave-lengths as measured at B and A respectively, 
the last equation resulting from differentiation of the second member of (30). 
If w(t) is the function occurring in formula (1g) we have the relation (20) 
w(O(t)) 
v(t A 32 
(t) ot) (32) 
On substituting this in (31) and using (30), we obtain 
Ty) 
Dw ; ( 
w(t,) 33) 
If we write t,;=t)—7, we can expand the right-hand side of this equation in 
powers of 7, thus: 


D=14 : .- + 7? <n ene apa (34) 


where we have written w, w’, etc., for w(t,), w'(ty), ete. now remains to 
relate 7 to the distance AB. 
‘rom equation (32) we obtain 
dt “O'(t) dl 
| ct) , w(H(t)) cay 8) 
and therefore 
" dt ; mi dt 
J 4, w(t) J, A(t) 
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If we introduce a new function s(t), defined by the equation s‘(t) = 1/a(t), and 
write t, =t, —o, then equation (35) reduces to 
8(ty) — (ty — 7) = 2[s(ty) — s(t, — 7)], 


both sides of which can now be expanded, thus 


; 9 ” a ) / 3 "9 , 
a” 2u) * ww T 7” @ Tv” 25 * WwW Ww 
——+...=2 + ——, + > —- ti... |. (36) 
hy w” w 2 w* 6 w 


Now, by the usual definition (in “‘ Kinematic Relativity’’) of the distance M 
between A and B, we have o=2M/c, so that equation (36) can be solved for + 
in terms of M, giving 

M I VM? wy’ I x (= 


2 > 


( at" @ 2 ¢ Ww w 


which, when substituted in equation (34), leads to the following shift-distance 
relation : 


I VWs ne wo" 


\ | ( a , 
D he 6a > eves (37) 
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Summary 


A modified form of Boyle’s Law is derived for a spherical mass of gas at 
uniform temperature in equilibrium under its own gravitation and an external 
pressure, and a typical p-V’ curve is shown. Any such mass of gas is gravi- 
tationally unstable if sufficiently compressed. It is shown also that a sufh 
ciently large isothermal gas sphere is unstable, and the condition for this ts 
compared with Jeans’ formula for gravitational instability 

The relevance of the phenomenon to the formation of the nebulae is briefly 


considered 


1. Introduction.—\t has recently been suggested by ‘Terletsky (7) that for a 
large mass M of gas, of volume V and temperature 7, containing NV molecules 
under boundary pressure p, the equation of state should be not 

pV =NkT (1.1) 
but 
pV = NRT -2GM?V-"3, (1.2) 
where & is Boltzmann’s constant, G is Newton's constant of gravitation, and & is 
a constant depending on the shape of the mass. ‘lhe proposed correction otf 
Boyle’s Law arises because, for a large mass, one has to take account of the gravi 
tational interactions between the molecules. 

‘Terletsky is interested in the statistical fluctuations in density of the gas, and 

he uses the ordinary formula for the volume fluctuation 
(V—V/y kT 
pe V*p/aV)’ wm 


where V is the mean volume. If this is applied to (1.1) the right-hand side 
becomes 1/N and the volume (or density) fluctuation diminishes as the mass 
increases ; if, however, it is applied to (1.2), one finds 


7_ Vy (Gm? N'|-! 
(F ak » en | , (1.4) 


~— 3 kT" 
where m=M/N. From (1.4) it follows that the fluctuation can become large 
for large masses as well as for small ones. For example, supposing that the gas 
is hydrogen at 10° deg. K, and is a uniform sphere (in which case « = }(47)!*, as will 
be shown later), one finds that the right-hand side of (1.4) becomes large when 
the ratio mass/radius is about 2 x 10”, 

This result is especially interesting because of the present situation in 
cosmology. One of the difficulties of the point-source cosmological models of 
general relativity concerns the formation of the nebulae. If one supposes that the 


V2 oN 
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primordial matter was uniformly distributed but subject to ordinary fluctuations, it 
seems from the work of Lifshitz (6) that the time of 4 x 10° years which has elapsed 
since the start of the universal expansion is not long enough to allow the develop- 
ment of the present degree of heterogeneity. Some mechanism which allows 
large fluctuations in the cosmic medium might help to explain the formation of 
the nebulae. 


In discussing formula (1.2) ‘Terletsky seems to have in mind a static mass of 
gas of uniform density. ‘This is unsatisfactory because if the mass is finite the 
equation of hydrostatic equilibrium is not satisfied, and if it is infinite the gravi- 
tational field strength at a point is not properly defined and the use of Newtonian 
gravitation theory is questionable. ‘here is also the objection that in a proper 


treatment of fluctuations one should consider, as (¢p/0V), tends to zero, the 
higher derivatives such as (0*p/0V*),, as is done in the theory of gases near the 
critical point. Further, if for a certain value of the ratio mass/radius (ép/oV ), 
changes sign (as it does for (1.2)) this suggests that such a mass of gas would be 
unstable, and not that it would be in a critical state as ‘Verletsky proposes. 

In this paper I obtain the equation of state for a large gas sphere of constant 
mass and temperature, but non-uniform density, and | construct its p-V curve. 
This amounts to an investigation of the equilibrium under its own gravitation 
and an applied external pressure of an isothermal gas sphere containing a fixed 
mass. Boyle’s Law (1.1) 1s obeyed only as V+ #. If the mass/radius ratio 1s 
large enough the sphere is unstable and the gas tends to collapse towards the 
centre; there is, however, no evidence that gravitation, though responsible for 
the instability, could confer on the gas the properties of a critical state. 

The p-V curve is obtained in Section 2, and this is followed in Section 3 by a 
discussion of the stability and a comparison of the results with those of Jeans on 
yravitational stability. ‘The paper ends with Section 4, in which the main 
conclusions are briefly summarized, 

The p-V curve. Consider a spherical mass of perfect gas (of uniform 
chemical composition) in isothermal equilibrium under its own gravitation, and 
let O be the centre. \t eve ry point P the gas obeys the local equation of state 


p=kpT/m, (2.1) 


where p, p and 7 are the pressure, density and temperature at ?, and mis the mass 
of a molecule. ‘he gas satisfies also the equation of hydrostatic equilibrium 
dp 4nGp(r 


. »y"dy 
dr r* },? alg 


1d (ce CO 
red pP a) a ” 


eqns. (2.1) and (2.2) together give the equation for the density distribution ot a 


1d ‘i: 47Gmp 

> ra ° beat 

r-dr\p tr) k1 
his equation was studied in detail by Emden (2), and is discussed in standard 
textbooks. (See, for example, Chandrasekhar (1).) 


isothermal gas sphere: 


eqn, (2.3) may be put into a form analogous with ‘Terletsky’s equation (1.2) 
in the following way. Consider a spherical volume V, with centre O and radius 
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r, containing N molecules and under pressure p at the boundary. We have 


an ee oe 
N=" | erry, (2.4) 
and we may write 

4nkT (7 


py*dy 
m 0 


pV -NkT= "rp 


3 


using (2.1). Using (2.2) we now find 


47 (" s% 1, 
\9 dy ry, 


6 r V 
pV =NkT- : 7G pady | px*dx., (2.5) 


This is the corrected form of (1.1) when gravitation is taken into account, p 
being given by (2.3). 

Although we have obtained it for an isothermal gas sphere by using the local 
equation of state (2.1) and the equation of hydrostatic equilibrium (2.2), 
equation (2.5) may be derived without assuming either (2.1) or (2.2). If one uses 
the virial of Clausius, (2.5) follows on the assumptions, first, that the spherical mass 
of gas is in a steady state, secondly, that its mean kinetic energy is equal to 
(3/2) NRT, and thirdly that one may in computing the virial use ordinary Newton- 
ian gravitation theory. Included in the last assumption is that the gravitational 
force at a point distant 7 from the centre arises only from matter within radius r 
Strictly, therefore, p does not need to satisfy (2.3) for (2.5) to be true, and 
Terletsky’s expression (1.2) (with «= 1/5(47/3)'*) follows from (2.5) by taking 
p=const. ‘The objection to this, as already stated, is that if the sphere is finite, 
p=const. does not satisfy the equation of equilibrium, and if it is infinite the 
gravitational force is not defined, and the third assumption above is not justified. 

Let us now transform (2.3) by making the following substitutions : 


p=Ae", r= Bh2)-N2E, (2.6) 
where A is an arbitrary constant, and 


B=kT/4nGm., (2.7) 


a4 (2M) _ ov 
3 (eg) <<. (2.8) 


The boundary conditions to be satisfied at the centre r=o are 


‘Then (2.3) becomes 


P= Pes 


where p, is the central density, supposed given. If we choose A=p, these con- 
ditions become, in the new variables 


. dip 
a =0, dé 


=0 at =0, 


which completely determine the solution of (2.8). ‘The solution, so determined, 

may be applied to gas spheres of different central densities by choosing different 

values of A. ‘Thus we may regard A as the parameter of the family of isothermal 

gas spheres. Any complete isothermal gas sphere extends to infinity: this means 

that in (2.8) e~’ (and therefore p and p) has no zero for finite positive values of €. 
25 
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To construct the p-V curve, let us imagine a Boyle’s Law experiment con- 
ducted in a spherical container of variable radius which contains a fixed number 
N, of molecules of gas, and which is kept at constant temperature 7. If the 
pressure p, on the outside is successively varied, and after each variation the gas 
is allowed to come to equilibrium, it will assume the configurations of the central 
cores of a series of isothermal gas spheres. ‘To each point (py, Vo) (Vo being the 
volume of the container) of the Boyle’s Law curve will correspond a different 
value of A. 

Although it does not seem possible to find a closed formula for py in terms of 
V,, at constant N, and 7, we can find the derivative (dp)/0V5)y, 7 in the following 
way. From (2.2), (2.4), (2.6) and (2.7), we have 

yw 
N= . pe2r-Wege - (2.9) 
We are interested in variations of py and V, at constant N,, so we take 


4” case | y-102 3 (\ | 
ON 8/25 < A-N2E 2 oO, 2.10 
0 m perc 7) (Z “a ( ) 
where the suffix o means that the value of a quantity on the boundary is to be 


taken. [quation (2.10) gives 


\3 Lf . difs 12 d ‘ dbs 
), Bhthe 7% OA +A -"* >» O€ 2. 
os di), dé dé | Jy S09» sandal d 


| d ge | fe" 
2 _ oe °. 
\déL> déjf, °° 


0 0€& 
Now éf + ( x) oA; 
} OA _ 


substituting this and (2.12) into (2.11), and using (2.6), we find 


. difs nee or 2 
te) an} sA-* 6(F) + ASE 8 © "| + B-**E,? Org 
77 0 


which gives, on using (2.6) and (2.9) 


and from (2.8) 


2Ap,0r 0 
OA 


m N,/477r," VoPo 
Krom 2.1 we have 
, 


( Opy ) k ( Opo ) AT (2) 
OVo) N,.1 m \oV, N,, T 4nmr,,” \ or, Ne, r 


dpy = O(Ae”*) = e "0A — Ae (Ze), Bo 


apt J ‘ dfs’ dy dip\ 0. | 

; io | = ( dé), on : dé} 4 Or, ror 
al mN, mN,A 
| - ( E) . 477,*B ros : a 5) 


where (2.6) and (2.9) have been used. Reintroducing p from (2.6) into (2.15) 
and using (2.13), we find that (2.14) gives after simplification 


3V," «NRT 
3PoVo 


_ (**) ‘Gm*N« 
( Opy 2Po 3/ — Opy Vo" . 
).. 1 


V4 
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If in (2.16) one puts G =o and supposes that the density is uniform one finds the 
usual expression for op/eV for an ideal gas. 

Although I have not been able to integrate (2.16) in finite form, it is not 
difficult to sketch the solution curves. It is rather easier, however, to obtain 
the solution of physical interest from Emden’s data (2) in the following way. 
‘This procedure amounts to choosing that integral of (2.16) in which py, is of order 
V, for large V4. 

Emden gave the solutions of (2.8) in the form of a table of values of ys and 
djs/d€ corresponding to certain values of the independent variable €. If one 
assumes a definite constant value of A one can find from the table the density 
distribution for an isothermal gas sphere of given central density A by using (2.6). 
Instead of doing this, let us allow A to vary whilst keeping N (number of mole- 
cules inside radius r) constant; in this case we have from (2.9) 
dp 

dé" 
Then points on the p-V curve are obtained by expressing p(=kTp/m) and 
V( =‘nr’) in terms of ys and djs/dé, and using Emden’s data to calculate them. 
From (2.6), (2.7) and (2.17) we find 


MT (4) 
47G*m®N* ii) 


47 GmN °}: 
3 kT Edds/dé 


Using this method we find p-V curves such as that shown in Fig. 1. which 


2.17) 


12 4" 3/2 
— mN poss 


refers to a mass of 10" gm. of hydrogen at o C. For large volumes the curve 
approaches the ordinary curve p= Nk7/V, but for larger values of p it spirals 
round the singular point 


2 af Gm? Nn) 
P= Gm N®’ ot as 


of the differential equation (2.16). 


ormesere| 


80 
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3. Stability.—For points on the curve between V = » and A, the highest 
maximum of p, op/dV is negative and the equilibrium is stable. At A the equi- 
librium is unstable because a small increase in the external pressure causes a 
decrease in the volume which results in a reduction of the internal boundary 
pressure as the state of the gas moves along the part AB of the curve. ‘Thus an 
isothermal sphere of gas, if sufficiently compressed, is unstable, because a small 
increase in external pressure will make it collapse. 

It follows also that any sufficiently large isothermal gas sphere is unstable. 
For consider a sphere of radius R in equilibrium under the correct boundary 
pressure, and with definite central density. Let p(r) be the pressure at radius 
r(<R) and V(r) and N(r) the volume and number of molecules inside r. 
Throughout the sphere there will be slight pressure and volume fluctuations, 
and for spherically symmetric fluctuations (0p/0V), 7 at radius r is given by 
(2.16). Now the numerator and denominator of the second fraction on the right- 
hand side of (2.16) may be written 


/2)"83Gm2N2 dis\2 
_— (47/3)'8GmEN? 4er(F) (3-1) 


pV 
kTN (a3) 
—-—, =! 3.2 
3pV 
in which form the two expressions may readily be calculated from Emden’s data. 
For smaller values of € (or r) both these functions are positive and dp/dV is nega- 
tive as in Boyle’s Law. For larger values of r both (3.1) and (3.2) oscillate in 
sign, (3.1) becoming zero before (3.2). Let V =v be the smallest volume for 
which (3.1) (and therefore 0p/0V) becomes zero; then one easily finds 


O*p 
sp2 oO. 
ab N, T(V =v) 


‘Thus for a small fluctuation dV in the volume of the sphere of mean volume 2, 


I 


we have 


O*p 
dp ( ) dV24+ O(dV*) <o. 


dV?) y 
‘Therefore if a fluctuation occurs which results in a slight decrease in the volume 
of this sphere, this will produce a decrease in the pressure inside its boundary 
which will lead to a further reduction in its volume. ‘Thus an isothermal gas 


sphere of volume greater than v will be unstable because small fluctuations will 
make the volume wv collapse towards the centre. 

‘The lowest value of € at which 0p/0V becomes zero is found from (3.1) and 
I:mden’s table to be about 6:5. Hence from (2.6) the critical radius is 


kT 1/2 
6-612) -12 8( , 
-) an) (3-3) 


where, it will be recalled, A is the central density. Formula (3.3) may be written 
in various forms. For example, we find 


kT \}2 
‘ 0°-76( . (3.4) 
mGp : 


where p is the mean density inside the sphere of radius r, ; 


kT \¥2 
or r= 0°49 (<a (3-5) 
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where p, is the density at the critical radius. Finally the critical ratio of 
mass radius is 
M(r,) kT 
= 2°4 


r, mG” (3-6) 

Reverting to the p-V curve for a fixed mass of gas we can now see that any 
point on the spiral part of the curve below A refers to a state of unstable equilib- 
rium. (It is necessary to point this out because along BC, for example, op/eV is 
positive, so it might seem that this part of the curve refers to stable states.) ‘The 
point A refers to the critical value of € (€,= 6-5), and on the spiral part below A, 
for example at D, € €,. ‘Vhus for an isothermal gas sphere represented by D, 
there must be some interior value r, of r at which €=€,, and the gas inside r, is 
unstable to small fluctuations. It therefore follows that the only stable part of 
the p-V curve is the part from V = o to A. 

Finally, | shall compare the criterion of gravitational instability derived here 
with that discovered by Jeans (4). Jeans considered a mass of gas, uniform 
except for small condensations and rarefactions which are supposed to propa- 
gate through it as waves of length /. He found that the velocity of propagation 


dp Gpl? 2 
dp 7 


7 dp 
Gp dp E 


was 


This becomes imaginary if 

(3-7) 
and Jeans concluded that if the linear dimensions of the mass of gas were greater 
than (7/Gp dp/dp)'* the motion would be unstable and the mass would break up. 


If one considers a sphere and uses p=k7'p/m, one finds that to the three criteria 
for instability (3.3), (3.4) and (3.5) there corresponds the single equation 


kT \'2 
4 (3.3) 


and to (3.6) corresponds 


kT 
3°3—7- 


(3.9) 


It seems worth while to emphasize that the Jeans’ condition for instability of 
a gaseous mass is (3.7) and not, as some writers seem to think, that the gravitational 
energy (taken as positive) shall exceed twice the thermal energy. If a gas is under 
no forces except its own gravitation, then, from the virial theorem, a necessary 
condition for it to be in a steady state is that the gravitational energy shall be equal 
to twice the thermal energy. It happens, as Jeans pointed out, that the equilib- 
rium dimensions, as given by the virial theorem, are about the same as the maxi- 
mum dimensions for a stable mass as given by (3.7). Hence Jeans concluded 
that an isolated mass in equilibrium would have no tendency to break up. ‘The 
situation is quite different if the mass under consideration is not isolated: if, for 
example, it consists of a local concentration in a gaseous medium which exerts a 
pressure on it. In this case the virial theorem does not apply, but the mass might 
still be unstable. 

From a comparison of formulae (3.3) to (3.6) with (3.8) and (3.9), it seems 
that the instability investigated in this work occurs in masses with sizes and 
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densities roughly the same as those predicted by Jeans’ formula. However, the 
question is not so straightforward as it appears, for the following reasons. 

In the first place, Jeans’ formula is open to the criticism, already mentioned, 
that one should not assume a constant density for a gravitating mass of gas, whether 
it is finite or infinite. Secondly, it does not follow that Jeans’ method will give 
similar results when applied to a non-uniform mass: in fact, Jeans (§) came to 
the conclusion that similar perturbations applied to an isothermal gas sphere did 
not make it unstable. ‘Thirdly, the numerical factors in equations (3.4) and (3.6) 
are important because, owing to the non-uniform density, r6'* and M(r)/r are 
not monotonically increasing functions of r; they have in fact maxima and 
minima, and tend to limits as r -- #%, so that they each lie within a certain numerical 
range. ‘lhe upper limit of the range is in each case less than the value from (3.8) 
or (3.9), in the direction of which p is assumed to be constant. ‘hus, for example, 
the critical ratio M(r,)/r, in (3.9) is not achieved for any isothermal gas 
sphere. 

4. Conclusion.—¥rom Sections 2 and 3 it appears that, at any rate for an iso- 
thermal gas sphere, the corrections to Boyle’s Law proposed by ‘lerletsky result 
in gravitational instability, and not in the gas having the properties of a critical 
state. Although one condition for a critical state, i.e. 0p/AV =o, is satisfied at an 
appropriate volume, the second condition, 0°p/0V* =o, is not, and the fact that 
the second derivative is non-zero implies instability. 

‘The instability is interesting because the condition for it resembles the well- 
known Jeans’ condition, but it is here obtained without assuming that the gas 
has a uniform density. ‘The difference between the two conditions suggests that 
Jeans’ formula should be used with caution, and with due regard to the fact that 
no large cloud of gas actually existing in nature satisfies the conditions assumed 
in deriving it. 

‘The bearing of this work on cosmology is as follows. According to the point- 
source models of relativistic cosmology, there was, during the early history of the 
universe, a finite or infinite expanding mass of gas which at any given time was 
of uniform density. In order to find out whether such a model can explain the 
formation of the nebulae two additional theories are required: first, one which 
gives the probability that density fluctuations of a given magnitude shall occur ; 
and secondly, a theory which describes the motion and growth of such pertur- 
bations once they have occurred. 

A rigorous formulation of the first theory, which would involve an adaptation 
of statistical mechanics to general relativity, is not available, and for want of 
something more appropriate, one naturally falls back on classical statistical mech- 
anics, and obtains the likely fluctuations by the formula (1.3). One then uses 
these in the second theory, available in the work of Lifshitz, and finds that the 
probable fluctuations could not have given rise to the nebulae in the time available. 
‘This implies that in the point-source models the formation of the nebulae is vastly 
improbable. 

‘To save the point-source models one might make some special hypothesis about 
the early history of the universe ; for example, that the primordial cosmic medium 
was turbulent, as has been suggested by Gamow (though it is not in fact clear 
that primordial turbulence would explain the formation of the nebulae); or one 
might even suppose that the universe started its expansion complete with in- 
homogeneities of the required size, at some time slightly later than the singular 
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state of the point-source models. Alternatively, one may examine more care- 
fully the application to cosmology of classical fluctuation theory. 

It was as a contribution to this latter examination that the present work was 
undertaken. If it had turned out that large density fluctuations could take place 
without upsetting the stability of a system such as an isothermal gas sphere, this 
would have been some indication that a rigorous treatment of fluctuations in the 


expanding universe might yield perturbations of the magnitude required to 


explain the formation of the nebulae. In fact, however, the investigation lends 
no support to this idea because it reveals no large fluctuations, but shows that the 
isothermal gas sphere is unstable. 

The instability could be of interest only for the later nebular evolution. An 
isothermal gas sphere of nebular mass represents a considerable inhomogeneity, 
and for its formation within the appropriate time a large initial fluctuation in the 
homogeneous cosmic medium would be required. ‘Uhus the instability by itself 
does not solve the problem because one first needs a mechanism to produce 
density concentrations corresponding to the isothermal gas spheres. 

Granted some mechanism to produce the isothermal gas spheres, the criterion 
for instability fits quite well the hypothesis that the nebulae were formed from 
them. For example, if one takes p = 10-** g/cm® one finds that the biggest stable 
sphere of hydrogen has the mass of a nebula (10"g) if the temperature is about 
10° deg. K. ‘This agrees well with the temperature of a condensing nebula 
deduced by Hoyle (3) on other grounds. 
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